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Page  2.  (Ho  Typing). 

AERODYNAMIC  CHARACTERISTICS  OF  BODIES  IN  THE  TRANSITION  REGION  AT 
HYPERSONIC  SPEEDS  OF  FLOW. 

V.  N.  Gusev,  T.  V.  Klimova,  A.  V.  Lipin. 

SUMMARY. 


Is  conducted  the  analysis  of  the  special  features/peculiarities 
of  the  hypersonic  flow  around  bodies  in  the  transition  region,  which 
lies  between  the  region  of  free  molecular  flows  and  the  zone  of  flow 
of  gas  as  continuous  medium.  Are  discussed  questions  of  the 
simulation  of  actual  conditions  in  the  transition  region  on  the  basis 
of  the  vast  experimental  material,  obtained  in  the  low-pressure  wind 
tunnel  of  TsAGI  [lIArH  -  Central  Institute  of  Aerohydrodynamics  im.  N 
Ye  Zhukovskiy].  Work  gives  the  aerodynamic  characteristics  of  the 
broad  class  of  bodies. 

ADOPTED  DESIGNATIONS. 

cj  -  aerodynamic  coefficient; 

<,-—73 —  -  drag  coefficient; 

fee  uoz  _ 


DOC  =  83067301 


PAGE  3 


d  -  diameter 


D  -  diameter 


Ji  £y/tfjp  — 


M, 


*“p.«2 


n  -  exponent 


Pr  -  Prandtl 


lift  coefficient; 


of  blunting; 


of  maximum  cross  section; 


lift-drag  ratio; 


coefficient  of  pitching  moment; 


in  law  f i—r"; 


number; 


Re,  -  Reynolds  number,  in  whom  the  coefficient  of 

viscosity/ductility/toughness  is  calculated  from  the  temperature  of 
stagnation. 


Reoi  — 


Poo  ^oo  L 


»*0 


I 


Poo  ^oo  £ 
1*0 


Re, 


■oo  ' 


P,t/ccP 

f*0 


S  -  area; 


r#*  TV  -  temperature  of  stagnation  and  the  temperature  of  body 
surface; 


-  speed  in  the  undisturbed  flow; 


DOC  »  83067301 


PAGE  4 


a  -  angle  of  attack; 

«BI,  ~  angle  of  attack  when  K  =  Kmu; 

y  -  specific  heat  ratio; 

9  -  half-apex  angle  of  the  cone  and  semicone; 

5  -  thickness  of  plate, 

»t  -  »//; 

^  -  elongation  of  plate,  equal  to  the  ratio  of  its  width  to  the  chord 
in  the  root  section; 

x  -  sweep  angle; 

Pao  -  density  in  the  undisturbed  flow; 

u  -  coefficient  of  viscosity/ductility/toughness; 

(i,  -  coefficient  of  viscosity/ductility/toughness  with  T=T, . 

Indices  "t"  and  "n"  designate  respectively  wind  tunnel  tests  and 


full-scale. 
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Introduction. 

The  technical  progress  in  aviation  and  rocket  engineering  led  to 
the  intense  development  of  theoretical  and  experimental  studies  in 
the  region  of  aerodynamics  of  hypersonic  flows.  The  greatest  number 
of  these  investigations  relates  to  two  sufficiently  to  well  studied 
regions  of  general  gas  dynamics.  One  of  them  -  this  is  usual  gas 
fluid  dynamics  in  which  the  characteristic  mean  free  path  of 
molecules  is  much  lower  than  the  significant  dimension  of  body.  The 
detailed  presentation  of  the  theory  of  hypersonic  flows  in  this 
region  with  its  numerous  applications/appendices  is  contained  in 
monographs  [1,  2].  Another  region  -  dynamics  of  the  free  molecular 
and  adjacent  it  medium  where  the  course  of  gas  proves  to  be  such 
rarefied  that  the  characteristic  mean  free  path  of  molecules  becomes 
much  more  than  the  significant  dimension  of  body.  Latest  achievements 
in  this  region  are  presented  in  monograph  [3].  Are  least  investigated 
flows  of  rarefied  gas  in  the  intermediate  transition  region.  A  strict 
theoretical  studies  of  such  courses  can  be  carried  out  only  on  the 
basis  of  kinetic  theory,  which  uses  an  equation  of  Boltzmann.  The 
solution  of  this  equation  is  at  present  connected  with  the  great 
mathematical  difficulties. 
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Not  easy  here  proves  to  be  the  experimental  path.  Creation  of 
hypersonic  low-density  flows,  their  diagnostics  and  conducting  in 
them  experimental  investigations  -  most  complex  task  of  contemporary 
experimental  aerodynamics  of  hypersonic  flows.  At  present  on  one 
experimental  installation  it  is  not  possible  to  carry  out  complete 
simulation  during  model  tests  in  hypersonic  flow,  since  besides  Mach 
numbers  and  Re  in  the  wind  tunnel  it  is  necessary  to  reproduce  the 
high  value  of  the  enthalpy  of  hypersonic  flows.  The  coincidence  of 
these  conditions  in  one  installation/setting  up  is  extremely 
difficult.  In  connection  with  this  experimental  research  of 
hypersonic  flows  divide  into  the  investigations  of  the  effects  of 
hydrodynamic  character,  caused  by  change  Math  numbers  and  Re,  and 
investigation  of  the  effects  of  imperfect  gas,  caused  by  high  energy 
of  flow.  To  the  discussion  of  hydrodynamic  special 

features/peculiarities  indicated  above  of  the  hypersonic  flow  around 
bodies  in  the  transition  region,  which  lies  between  the  region  of 
free  molecular  flows  and  the  zone  of  flow  of  gas  as  continuous 
medium,  and  is  dedicated  this  work.  In  it  is  systematized  the  vast 
experimental  material,  obtained  in  the  low-pressure  wind  tunnel  of 


TsAGI . 
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Simulation  of  actual  conditions  in  a  transition  region  1 . 


FOOTNOTE  l.  As  the  basis  of  this  part  of  the  work  is  assumed  V.  N. 
Gusev's  report  "Problems  of  simulation  in  the  dynamics  of  the 
rarefied  gases"  at  the  III  All-Union  conference  on  the  dynamics  of 
the  rarefied  gases  (Novosibirsk,  1969).  ENDFOOTNOTE. 


Questions  of  the  similarity  of  hypersonic  flows  in  the 

transition  region  were  examined  in  work  [4],  On  the  basis  of  the 

equation  of  Boltzmann  in  it  it  is  shown  that  under  the  power  law  of 

interaction  of  molecules  for  observing  the  similarity  in  the 

mode/ conditions  of  hypersonic  stabilization  besides  the  geometric 

similarity  it  is  necessary  to  maintain/withstand  the  equality  of  the 

following  parameters:  Reynolds  number  Re,,  in  whom  the  coefficient  of 

viscosity/ductility/toughness  is  calculated  from  the  temperature  of 

stagnation,  exponent  n  in  the  dependence  of  the  coefficient  of 

viscosity/ductility/toughness  on  temperature  temperature 

factor  TjT0,  where  T,  and  T9  -  respectively  the  temperature  of 

stagnation  and  body  surface,  specific  heat  ratio  y ,  number  of  Prandtl 
of  accommodation.  Of.  the  parameters  mentioned,  the  Refl  number 
Pr  and  coef f icientAit  is  basic  similarity  criterion:  with  its  change 

the  values  of  the  aerodynamic  characteristics  of  the  streamlined  body 
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can  change  to  the  orders. 

When  7'.T0=1  basic  laws  governing  the  supersonic  flow  around 
characteristic  bodies  over  a  wide  range  of  a  change  in  criterion  Re, 
were  revealed  in  work  [4].  Data  for  the  transient  region  were 
obtained  in  the  low-pressure  wind  tunnel  employing  the  presented  in 
work  [4]  procedure.  For  the  zone  of  flow,  close  to  free-molecular, 
when  is  applicable  the  theory  of  the  first  intermolecular  collisions, 
the  aerodynamic  characteristics  of  bodies  were  obtained  by 
calculation  [5,  6].  Experimental  data  for  the  case  of  viscous  flows 
of  continuous  medium  were  borrowed  from  works  [7,  8],  theoretical  - 
from  works  [9,  10].  For  the  cone  with  the  half-angle  of 
solution/opening  9  and  plate  whose  relative  thickness  5/L=0.03,  these 
data  are  cited  in  Fig.  1-6.  During  the  calculation  of  aerodynamic 
coefficients  as  the  characteristic  area  was  selected  the  area  of  body 
S  in  the  plan/layout;  pitching  moment  was  calculated  relatively 
"nose/leading  edge";  through  a  was  designated  the  angle  of  attack. 

As  it  follows  from  given  data,  the  results  of  experiment  and 
calculations  according  to  the  theory  of  the  first  intermolecular 
collisions  and  according  to  the  theory  of  the  viscous  flows  of 
continuous  medium  make  it  possible  to  obtain  information  in  the 
entire  transient  zone  of  flow.  In  this  case  experimental  and 
theoretical  data,  obtained  taking  into  account  the  parameter  of 

similarity  Re0,  will  agree  sufficiently  well  with  each  other. 
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Fig.  5.  Key:  (1).  Plate. 


Should  be  focused  attention  on  the  nonmonotony  of  a  change  of 
some  characteristics  in  the  transition  region  where  their  values 
prove  to  be  more  than  the  values,  obtained  in  the  free  molecular 
flow,  and  in  certain  cases  except  the  maximum  indicated  they  have 
even  and  a  minimum  (for  example,  the  coefficient  of  lift  force  of 
cone).  All  th as  speaks,  that  the  laws  of  the  behavior  of  aerodynam 
characteristics  in  the  transition  region  are  very  complicated  and 
they  cannot  be  obtained  by  simple  interpolation  of  given  for 
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continuous  medium  and  free  molecular  courses. 

During  the  simulation  of  actual  conditions  in  the  wind  tunnels 
besides  the  basic  criterion  of  similarity  Re,  it  is  necessary  to 
maintain/withstand  a  whole  series  of  other  important  similarity 
criteria.  However,  the  installations/settings  up  existing  at  present 
in  which  the  temperature  of  stagnation  is  close  to  the  room,  by  no 
means  provide  this  simulation.  In  connection  with  this  arises  the 
question  about  the  procedure  of  the  recalculation  of  the  results  of 
tube  experiment  for  the  actual  conditions  with  the  noncoincidence  of 
some  similarity  parameters.  In  the  case  of  the  thermodynamically 
ideal  gas  these  parameters  include  7\,/T*  by  n,  y,  a  number  of 
Prandtl,  accommodation  coefficients.  As  are  shown  available  in  works 
[11,  12]  data,  under  the  conditions  of  the  viscous  hypersonic  flows 
at  the  speeds  of  less  or  the  order  of  the  first  space,  effects  of 
imperfect  gas  to  the  aerodynamic  characteristics  are  insignificant, 
and  air  can  be  considered  as  the  thermodynamically  ideal  gas.- 

It  is  known  that  for  air  at  low  temperatures  (T<150°K)  jx~T, 
i.e.,  n«l,  and  with  T>400°K  value  n  is  close  to  0.67.  In  other  words, 
at  low  temperatures  of  the  molecule  of  air  they  behave  as  Maxwellian, 
i.e.,  with  large  T  they  are  nearer  to  the  elastic  spheres. 
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Therefore  obtained  in  the  low-temperature  wind  tunnels  data  cannot  be 
transferred  to  nature,  even  if  we  ensure  a  sufficient  cooling  of 
models  for  obtaining  the  full-scale  values  TJT0. 

For  the  courses,  close  to  the  free  molecular  ones,  a  question 
about  the  recalculation  of  tube  data  to  the  full-scale  ones  was 
examined  in  monograph  [3].  In  the  region  of  viscous  hypersonic  flows 
for  this  purpose  was  applied  the  approximate  law  of  similarity  of 
Cheng  [11,  13].  The  detailed  analysis  of  the  effect  of  different 
similarity  criteria  on  the  aerodynamic  characteristics  of  simple 
bodies  in  this  region  is  given  in  works  [14-16].  The  results 
represented  in  them  of  the  parametric  analyses  of  the  flow  around  the 
simplest  bodies  make  it  possible  to  explain  the  degree  of  the 
nearness  of  the  aerodynamic  characteristics,  obtained  under  the  wind 
tunnel  and  the  actual  conditions,  and  in  certain  cases  [14]  to 
introduce  the  appropriate  corrections  during  the  complete  simulation 
on  criteria  TjTt,  rt ,  7. 

It  is  possible  to  attempt  to  find  the  approximate  laws  of 
simulation  with  the  aid  of  the  experimental  data  (for  example,  see 
[17,  18]).  In  this  case  one  should  only  have  in  mind  that  the 
boundary  of  free  molecular  hypersonic  flows  and  the  aerodynamic 
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characteristics  of  bodies  in  such  flows  depend  on  the  form  of  body, 
temperature  of  its  surface  and  laws  of  interaction  of  the  molecules 
between  themselves  and  body  surface  [3], 

Let  us  present  arbitrary  aerodynamic  value  r,  in  the  form 

r/  =  ria#  +  (rjo  —  ci  «)/£  Re*  (~7^ )  ’  ~TL '  n*  •  > 
where  c<«  and  cit  -  value  ct  with  Re0=®  and  to  zero,  and  /3  is  the 
function  of  the  similarity  parameters  and  in  the  general  case  depends 
on  the  shape  of  body.  Relying  on  experimental  data,  it  is  possible  to 
attempt  to  find  such  value  0,  at  which  the  dependence  of  functions  f 
on  the  criteria  of  similarity  TJT0  and  n  will  be  weak.  In  this  case 
the  data,  obtained  in  the  low-temperature  wind  tunnel,  can  be 
converted  to  the  full-scale  ones.  Actually/really,  let  c[  -  obtained 
in  the  wind  tunnel  value  of  value  ct  when  TJT0=\  and  n«l,  and  c" 
corresponding  to  it  full-scale  value  ci  at  fixed  value  TjTt<t  1  and 
n=0 . 67 . 


Then 


C‘  3»  “f“  (fii  0 -  Ci  a)  f 


T. 


where 
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It  is  obvious  that  with  this  recalculation  it  is  necessary  to 
know  values  cl<t  and  •  obtained  under  the  wind  tunnel  and  the  actual 
conditions.  The  latter  can  be  determined  by  calculation. 

For  the  drag  coefficient  of  sphere  and  cone  the  identification 
of  parameter  0  was  produced  on  the  basis  of  experimental  data,  given 
in  the  works  [  ^  17].  Analysis  showed  that  with  /3=- 0.1  the 
experimental  values,  obtained  at  the  essentially  different  values  of 
temperature  factor,  will  agree  sufficiently  well  with  each  other  over 
a  wide  range  of  a  change  in  Mach  numbers  and  Re. 

Page  9. 

The  results  of  this  comparison  are  given  in  Fig.  7  and  8.  There  solid 
lines  gave  the  results  of  calculation  according  to  the  theory  of 
first  intermolecular  collisions  [5]. 

It  is  obvious  that  approximation  method  examined  above  of  the 
recalculation  data,  obtained  in  the  low-temperature  wind  tunnels,  by 
the  full-scale  ones  is  not  strict  and  cannot  encompass  entire 
diversity  of  bodies  and  ranges  of  a  change  in  the  similarity 
criteria.  However,  the  absence  at  present  of  a  sufficient  quantity  of 
experimental  data  at  substantially  the  higher  values  of  gas  enthalpy 
makes  it  possible  to  propose  nothing  the  best.  In  proportion  to  the 
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accumulat ion  of  such  data  this  gap/spacing  will  be  eliminated. 


It  is  necessary  to  keep  in  mind  that  in  certain  cases  with  the 
accuracy  acceptable  for  practical  purposes  there  is  no  need  for 
increasing  the  temperature  of  stagnation  of  flow  in  the  wind  tunnel 
to  its  full-scale  values.  For  example,  the  given  in  work  [14] 
estimations  show  that  in  the  region  of  viscous  interaction  heating 
gas  in  the  precombustion  chamber  to  T„=2000°  K  provides  the 
simulation  of  the  full-scale  values  of  resistance  and  lift  of  plate 
with  an  error  less  than  4%,  when  TJT0>Q,05.  The  latter  fact  in  many 
respects  facilitates  the  way  of  the  straight/direct  simulation  of 
actual  conditions  in  the  wind  tunnels. 
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[23]).  Are  represented  below  systematized  experimental  data,  obtained 
in  the  low-pressure  wind  tunnel  of  TsAGI.  The  coefficients  of  drag 
cx,  lift  c,,  pitching  moment  and  lift-drag  ratio  K  of  the  broad 
class  of  bodies  are  given  in  the  dependence  on  the  angle  of  attack  o 
when  r./r0-l  and  n*l.  Aerodynamic  coefficients  are  everywhere  related 
to  velocity  head  q  and  characteristic  area  S,  pitching  moment  is 
related  to  the  significant  dimension  of  L  and  is  calculated  relative 
to  the  end  section  of  model. 

Tests  were  conducted  in  the  low-pressure  wind  tunnel  of  TsAGI. 
For  the  formation  of  the  supersonic  flow  served  calculated  conical 
nozzle  with  the  aperture  angle  of  25°.  As  the  working  gas  was 
utilized  air  at  a  room  temperature.  A  maximum  change  In  Mach  number 
along  the  length  of  model  did  not  exceed  2%. 

During  the  determination  of  aerodynamic  coefficients  were 
utilized  the  values  of  the  parameters  of  the  undisturbed  flow  at  the 
point,  which  corresponds  to  the  middle  of  model.  As  applied  to  the 
low-pressure  wind  tunnels  the  procedure  of  the  recalculation  of 
aerodynamic  characteristics  in  the  nonuniform  flow  the  not 
corresponding  characteristics  in  the  uniform  for  some  simplest  bodies 
in  the  maximum  modes/conditions  of  flow  is  presented  in  the  works 
[24,  25j  For  the  cone  with  a  half-angle  of  solution/opening  of  0=30° 
at  the  angle  of  attack  a*0  the  following  from  work  [24]  correction  to 
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cx  for  the  heterogeneity  of  flow  in  given  conditions  of  experiments 
did  not  exceed  4%.  For  the  cones  with  the  smaller  aperture  angles  it 
was  less. 

For  measuring  of  the  aerodynamic  forces  and  moments  of  forces 
were  utilized  three-component  magnetoelectric  weights.  The 
sensitivity  of  balance  comprises  the  tenths  of  milligram.  A  relative 
error  of  measurement  of  forces  does  not  exceed,  as  a  rule,  ±2%,  the 
moments  of  forces  -  ±5%. 

Taking  into  account  errors  of  measurement  of  the  flow  parameters 
a  relative  error  in  the  determination  of  aerodynamic  coefficients 
comprised  in  the  majority  of  cases  of  5-6%.  Somewhat  greater  errors 
occurred  for  the  models  of  small  sizes/dimensions  -  they  reached 
7-10%. 


The  installation/setting  up  of  symmetrical  models  relative  to 
the  axis  of  flow  was  realized  with  the  aid  of  the  mechanism  a  of 
weights.  Angle  of  attack  at  was  considered  equal  to  zero  when 
simultaneously  they  turned  into  zero  normal  to  the  model  force  and 
pitching  moment.  An  error  of  measurement  a  composed  10'.  With  the 
same  accuracy  the  "angular  dimensions"  of  models  accepted  correspond 
to  real  ones. 


DOC  *  83067301 


PAGE 


Cone.  The  aerodynamic  characteristics  of  the  thin  pointed  cones 
with  different  aperture  angles  (7°s20S4O°)  at  the  fixed  value  of 
number  Re<>£=162,  cone  designed  along  the  length,  are  represented  in 
Fig.  9-12.  During  the  calculation  of  aerodynamic  coefficients  as  the 
characteristic  area  S  here,  as  subsequently  is  selected  the  area  of 
the  basis/base  of  cone,  and  as  the  significant  dimension  of  L  -  its 
length. 

The  analogous  characteristics  of  the  pointed  cones  with 
half-angles  of  solution/opening  of  6=30°,  45°  and  60°  at  the  fixed 
value  of  number  Reoo  =  65.  basis/base  calculated  according  to  the 
diameter,  are  given  in  Fig.  13-15. 

At  large  angles  of  attack  the  aerodynamic  characteristics  of  the 
pointed  cones  with  half-angles  of  solution/opening  of  0=15°,  20°  and 
30°  at  the  fixed  value  of  number  Reoo  =  97  are  given  in  Fig.  16-19. 

Page  11. 

The  effect  of  the  flat/plane  blunting  whose  relative 
size/dimension  d*d/D,  where  d  and  D  -  respectively  the  diameters  of 
blunting  and  basis/base  of  cone,  to  the  aerodynamic  characteristics 
of  slender  cones  with  different  half-angles  of  solution/opening  at 
the  fixed  value  of  number  ReJi  *162  for  0=2°. 5  and  3°. 75; 
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ReJi— 154  for  9-6°  and  ReSz  — 92  for  5  =  10°  is  illustrated  in  Fig. 
20-35. 


For  the  low  angles  of  attack,  i.e.,  when  the  dependences  of 
aerodynamic  coefficients  on  the  angle  of  attack  are  linear,  given 
higher  experimental  data,  which  correspond  to  number  Reo£  =  162,  are 
represented  in  Fig.  36-38  in  the  form  of  dependences  cx0,  c'y0  and 
on  the  half-apex  angle  of  the  cone  9  at  different  values  of  d  of 
flat/plane  blunting  {cx0,  c\o,  m,0  -  corresponding  values  cx,  dcyjdn ,  dmjdn 

with  a=o),  At  the  fixed  value  of  number  Reot  the  aerodynamic 
characteristics  of  cone  in  proportion  to  decrease  9  considerably 
exceed  the  values,  valid  under  the  conditions  for  inviscid  ideal 
flow.  The  effect  of  a  small  blunting  of  cone  on  its  aerodynamic 
characteristics  under  the  same  conditions  is  unessential,  at  least  in 
comparison  with  the  analogous  effect  during  the  ideal  flow. 

The  effect  of  spherical  blunting  on  the  aerodynamic 
characteristics  of  cone  over  a  wide  range  changes  in  the  parameters  6 
and  d  are  given  in  Fig.  39-70.  Criterion  of  similarity  ReSz>  for  each 

9  was  here  by  constant,  its  values  were  given  below: 

»  . 3°, 5  5°  r.S  108  12°  ,5  15*  17*, 5  20° 

ReJ0  . .  20  29  43  57  72  87  102  118 

Cylinder.  The  aerodynamic  characteristics  of  cylinders  (9= 0)  at 
different  values  of  L/D,  where  L  and  D  -  length  and  the  diameter  of 
cylinder,  are  given  in  Fig.  71-78.  As  the  characteristic  area  S  is 
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selected  the  area  of  basis/base,  as  the  significant  dimension  of  L  - 
length  of  cylinder.  At  the  fixed/recorded  length  of  cylinder 
(Reoi=162)  the  experimental  data  are  given  in  Fig.  71-74,  with 
fixed/recorded  diameter  (Re0To=13)  -  in  Fig.  75-78. 

Body  of  revolution  with  the  generatrix  Re~x/4.  The  aerodynamic 
characteristics  of  the  bodies  of  revolution,  which  have  the  equation 
of  generatrix  in  the  form  of  exponential  monomial  with  the  exponent 
of  are  given  in  Fig.  79-82  at  the  different  values  of  the 

relative  thickness  D/L,  where  D  and  L  -  diameter  of  basis/base  and 
the  length  of  body.  As  the  characteristic  area  is  accepted  the  area 
of  basis/base  S,  as  the  significant  dimension  -  length  of  body  of 
revolution  L.  Criterion  Reot,  designed  along  the  length  of  model,  for 

each  value  of  D/L  is  given  below: 

. 0.115  0,16  0,32  0,80 

Re0i .  228  162  162  98 

for  the  comparison  in  Fig.  79  and  80  dotted  lines 
plotted/applied  the  drag  coefficients  and  lift  of  acute  cone  at  the 
same  values  of  D/L  and  Reoz . 

Sphere.  Resistance  of  sphere  was  studied  both  theoretically  and 
by  experimentally  many  authors.  3asic  experimental  data  are  obtained 
in  works  [4,  17,  18]  and  are  given  in  Fig.  7.  The  values  of  the  drag 
coefficient  cx  of  sphere  with  the  needle  of  variable  of  length, 
pertaining  to  drag  coefficient  of  the  sphere,  in  hypersonic  flow  of 
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rarefied  gas  are  given  in  Fig.  83.  On  the  same  figure  for  the 
comparison  are  cited  analogous  data,  obtained  in  work  [26]  with  the 
large  Re  numbers.  It  is  evident  that  with  the  small  Re  numbers  a 
change  in  the  length  of  needle  weakly  affects  resistance  of  sphere. 

Page  12. 

Plate.  Are  given  below  the  dependences  of  the  aerodynamic 
characteristics  of  plate  on  its  geometric  parameters.  As  the 
characteristic  area  is  accepted  the  area  of  plate  S  in  the 
plan/layout,  as  the  significant  dimension  -  chord  length  in  the  root 
section  L. 

In  the  case  of  rectangular  plate  the  effect  of  its  relative 
thickness  8t  =  8  L  at  several  fixed  values  of  the  elongation  of  a  plate 
X  and  number  Reot,  designed  along  the  length  L,  are  illustrated  Fig. 
84-95,  and  the  effect  of  the  elongation  of  plate  X  at  constant  values 
H  and  Reoi  -  Fig.  96-104.  Fig.  105-107  presents  dependences  cx,  cy 
and  K  on  the  angle  of  attack  a  for  the  rectangular  plates  with 
constant  values  of  S  and  !,  =  8/1=0,028,  where  l  =  yis, over  a  wide  range  of 
change^.  Number  Reoi  here  was  calculated  from  the  reference  length  l 
and  was  equal  to  39. 

The  effect  of  sweep  angle  *  on  the  aerodynamic  characteristics 
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of  triangular  plate  at  different  values  Reor  and  is  shown  in  Fig. 

108-115.  The  values  of  the  aerodynamic  characteristics  of  triangular 
plates  over  a  wide  range  of  change  X  at  constant  values  of  S  and 

—  0,028,  where  /  =  are  given  in  Fig.  116-119.  Number  Reo/  here 
again  was  calculated  from  the  reference  length  l  and  was  equal  to  39. 

As  it  follows  from  given  experimental  data,  viscosity  effect 
significantly  changes  the  character  of  the  flow  around  body,  as  a 
result  of  which  the  dependence  of  the  aerodynamic  characteristics  of 
bodies  on  their  geometric  parameters  becomes  different  from  that 
which  is  observed  with  the  large  Reynolds  numbers.  For  example,  the 
given  in  works  [19,  20]  results  of  the  systematic  studies  of  the 
aerodynamic  characteristics  of  the  plates  of  the  fixed/recorded  area 
showed  that  with  a  change  in  the  elongation  X  the  dependence  of 
maximum  lift-drag  ratio  Kmn,  has  a  maximum  at  the  finite  value  X.  The 
results  of  these  investigations  are  given  in  Fig.  120-123  in  the  form 
of  the  dependences  of  maximum  lift-drag  ratio  and  angle  of 

attack  am„,  which  corresponds  Kmtx  on  the  elongation  X  and  relative 
thickness  5,  for  the  plates  of  different  planform.  As  it  follows 
from  given  data,  in  the  hypersonic  low-density  flow  the  wings  with 
the  elongations  X=0.4-0.8  have  advantages  from  the  point  of  view  of 
lift  effectiveness  in  comparison  with  the  wings  of  other  elongations, 
moreover  with  the  sufficiently  small  wing  chord  ratio  8, <0,03 
independent  of  its  form  optimum  is  wing  with  the  elongation  X=0.6. 
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The  aerodynamic  coefficients  of  the  plates  of  various  forms  and 
fixed/recorded  area  S,  streamlined  at  the  high  angles  of  attack 
(a*"90°),  are  represented  in"£able  1. 
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Key:  (1).  Coefficients.  (2).  Form  of  plates.  (3). 

Circle/circumf erence .  (4).  Ellipse.  (5).  Rhomb.  (6).  Rectangle.  (7) 
Equilateral  triangle. 
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Number  ReJ,  was  calculated  in  this  case  over  the  reference 
length  t and  was  equal  to  65.  As  with  the  large  Re  numbers  [27] 
the  form  of  flat/plane  blunting  virtually  does  not  affect  the  value 
of  aerodynamic  coefficients. 


Three-dimensional  bodies.  The  aerodynamic  characteristics  of 
three-dimensional  bodies  are  analyzed  based  on  the  examples  of 
semicone  and  combination  of  wing  with  the  semicone.  As  the 
characteristic  area  is  accepted  the  area  of  model  in  plan/layout  S, 


as  the  significant  dimension  -  chord  length  in  the  root  section  L. 
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The  effect  of  blind  sectors  of  body  on  its  aerodynamic 
characteristics  was  investigated  based  on  the  example  of  semicone.  In 
the  case  when  the  conical  part  of  the  semicone  is  turned  towards  the 
flow,  aerodynamic  coefficients  for  the  semicone  are  compared  with  the 
appropriate  characteristics  of  complete  cone.  When  Reoi  =  130 
corresponding  data  for  the  cone  with  a  half-angle  of  solution/opening 
of  0“15°  are  cited  in  Fig.  124-127. 

In  the  case  when  the  flat  surface  of  semicone  is  turned  towards 
the  flow,  the  aerodynamic  characteristics  of  semicone  with  5=15°  with 
number  ReJi*=130  are  compared  with  the  appropriate  characteristics  of 
triangular  plate  with  the  sweep  angle  of  x*75°  (Sts=0,025)  in  Fig. 
128-131. 


As  follows  from  given  data,  in  the  mode/conditions  hypersonic 
flows  of  rarefied  gas  the  effect  of  blind  sectors  of  body  on  its 
aerodynamic  characteristics  it  becomes  essential. 

Fig.  132-139  gives  the  values  of  the  aerodynamic  characteristics 
of  the  model,  which  is  the  combination  of  the  blunted  semicone  with 
the  tapered  wing.  The  sweep  angle  of  wing  of  x=75°,  chord  length  in 
the  root  section  is  equal  to  the  length  of  semicone  L,  and  the  length 
of  leading  edge  -  to  diameter  of  the  flat/plane  blunting  d  of 
semicone.  Number  Reoi=130,  angles  of  attack  a  are  considered 
positive,  when  the  flat  surface  of  model  is  turned  towards  the  flow. 
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Fig.  42. 
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Fig.  57. 


Fig.  58 
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Fig.  96 
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Fig.  107. 
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Fig.  120. 


Fig.  121. 


Fig.  120.  Key:  (1).  rectangular.  (2).  triangular.  (3).  elliptical 
(4).  rhombiform. 


Fig.  121.  Key:  (1).  rectangular.  (2).  triangular.  (3).  elliptical 
(4).  rhombiform. 


Fig.  122. 


Fig.  123. 
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ig.  122.  Key:  (1).  triangular  plate. 


Fig.  123.  Key:  (1).  triangular  plate. 


Fig.  124.  Key:  (1).  cone.  (2).  semicone 
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APPLICATION  OF  THE  MONTE  CARLO  METHOD  IN  DYNAMICS  OF  STRONGLY 
RAREFIED  GAS. 

V.  A.  Perepukhov. 

Summary. 

Are  described  setting,  method  and  some  results  of  solving  the 
specific  problems  by  the  Monte  Carlo  method  for  two  types  of  flows  - 
free  molecular  and  almost  free  molecular  (the  slight  disturbance  of 
free  molecular  flow). 

CALCULATION  OF  FREE  MOLECULAR  FLOWS. 

In  the  case  of  convex  body  the  solution  of  the  problem  of  free 
molecular  1  flow  does  not  represent  labor/work;  however,  during  the 
study  of  flow  by  the  free  molecular  flow  of  concave  bodies  and  bodes 
of  complex  form  it  is  necessary  to  consider  the  possibility  of 
contact  with  the  element/cell  of  the  body  surface  of  the  molecules  in 
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question,  reflected  from  other  elements/cells  of  surface  or  from 
other  bodies,  i.e.,  to  consider  interference. 


FOOTNOTE  1 .  Flow  is  free  molecular,  if  the  number  of  Knudsen 
Kn=A/5->®,  where  d  -  significant  dimension  of  body,  X  -  minimum  local 
mean  free  path  of  the  molecules  of  gas.  ENDFOOTNOTE. 


Analogous  situation  appears  also  during  the  calculation  of  internal 
free  molecular  flows  (for  example,  flows  in  the  channels),  when 
molecule,  "stray"  within  the  channel,  can  many  times  clash  with  its 
walls.  In  the  most  general/most  common/most  total  setting  such  tasks 
for  arbitrary  law  of  reflection  are  examined  in  works  [1-3].  The 
mathematical  formulation  of  the  problem  for  the  concave  bodies  is 
given  in  work  [1].  Fundamental  equation  is  here  the  integral  equation 
of  Fredholm  of  the  second  order  with  the  symmetrical  kernel  for  the 
particle  flux  to  the  element/cell  of  the  surface 


N  (dSJ  =  (dS,)  -f- 


J  N(dSJ 


G  (</5,  dS2)  dSt, 


(1) 


where  G(dSxdS,)  -  probability  that  the  molecule  reflected  from 
element/cell  dS,  falls  on  the  element/cell  of  surface  dSx.  The 
solution  of  this  equation  in  the  general  case  presents  great 
difficulties.  Moreover,  and  the  subsequent  calculations  of  the 
aerodynamic  characteristics  of  body  require  bulky  calculations  on 
computer (s) . 
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In  this  setting  in  works  [4,  5]  were  calculated  the  aerodynamic 
characteristics  of  the  circular  cylindrical  and  spherical  surfaces, 
converted  by  concave  side  to  the  flow,  which  encounters  at  the  angle 
of  attack  a  (Fig.  1).  Since  in  these  works  it  was  disregarded  with 
the  thermal  velocities  of  molecules  in  comparison  with  the 
macroscopic  gas  velocity  and  were  not  examined  the  cases,  when  some 
elements/cells  of  surface  were  shaded  by  others,  then  the  obtained 
results  were  valid  for  the  angles  of  attack  a,  which  satisfy 
relationship/ratio  a +u*rr/Z,  where  u  -  half-angle  of  the 
solution/opening  of  segment.  The  reflection  of  molecules  from  the 
surface  was  assumed  to  be  diffuse.  In  work  [4]  was  allowed  the 

a  (/  QQ  . .  • 

miscalculation,  which  then  was  corrected  in  [5],  If  ioe ~  ~YWTai  ^  ' ’ 
and  the  accommodation  coefficient  of  energy  ag=t.l,  then  the 
temperature  of  the  molecules  T7  =  T9ss;tob  reflected  and  the  effect  of 
concavity  can  be  disregarded/neglected.  In  this  case  the  interference 
of  bodies  in  the  free  molecular  flow  will  be  determined  only  by  the 
shading  action  of  bodies  on  each  other.  In  work  [ 6 j  was  undertaken 
the  attempt  compute  the  aerodynamic  characteristics  of  the  internal 
surface  of  hemisphere  with  arbitrary  Sa> ,  however  in  it  was  allowed 
error  in  mathematical  calculations  and  its  results  were  not  accurate. 
The  aerodynamic  characteristics  of  concave  cylindrical  surface  when 
•$»>!  were  designed  in  work  [7],  moreover  taking  into  account 
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shading/blanketing,  i.e.,  with  the  arbitrary  a.  To  the  study  of  the 

interference  of  the  groups  of  bodies  in  the  free  molecular  flow  is 

devoted  a  small  number  of  works.  Let  us  note,  two  of  them.  Work  [8] 

examines  the  task  of  flow  of  the  free  molecular  flow  about  two 

identical  plates  of  the  finite  dimensions,  perpendicular  to  each 

other  and  a. ranged/located  at  the  zero  angle  of  attack.  Work  [9] 

examines  the  task  about  the  free  molecular  flow  of  gas  in  the  flat 

ducts  and  the  lattices. 

% 

During  the  calculation  of  internal  flows  in  the  majority  of 

works  is  used  the  method  of  Clausing.  Fundamental  equation  in  this 

case  is  also  the  integral  equation,  analogous  to  that  led  above,  for 

the  flows  of  a  number  of  particles  to  the  element/cell  of  surface. 

For  the  tube  of  round  cross  section  this  equation  takes  the  following 

form:  i 

N  (/)  =  vV,  (/)  -f  f  G  (x)  dNt  ( x ).  (2) 

Here  l  -  length  of  tube;  N^/)  -  a  number  of  molecules,  which 
fly  without  the  collisions  with  the  surface  from  the  entrance  to  the 
output  of  tube;  N*  -  number  of  molecules,  which  fly  into  the  tube  and 
encountering  its  surface  at  least  one  time;  G(x)dN2(x)  -  a  number  of 
molecules  which  enter  the  tube  and  encounter  wall  in  interval 
x-(x+dx)  and  then  either  immediately  or  after  further  collisions  they 
leave  it.  Basic  difficulty  during  the  use  of  this  method  consists  in 
the  determination  of  the  form  of  the  function  G(x).  Examples  of  the 
calculation  of  internal  courses  by  the  method  of  Clausing  can  be 
found  in  a  number  of  works  (see  for  example,  [10-12]). 


DYNAMICS  OF  RAREFIED  GAS  AND  MOLECULAR  GAS  DYNAMICS(U) 
FOREIGN  TECHNOLOGY  DIV  WRIGHT-PATTERSON  AFB  OH 
25  AUG  83  FTD-ID<RS)T-0673-82 

F/Q  20/4 
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For  calculating  the  internal  flows  in  works  [13,  14]  is  used  the 
Monte  Carlo  method,  with  the  help  of  which  were  calculated  the  flows 
in  the  tubes  of  different  configuration  under  the  law  of  the  diffuse 
reflection  of  molecules  from  the  surface.  The  essence  of  the  Monte 
Carlo  method  in  connection  with  internal  flows  is  reduced  to  the 
following.  Let  there  be  the  tube  with  the  section  at  entrance  and 
at  output  S2,  then  the  gas  flow  through  tube  N*N2P2 2-N2P2 2 ,  where  N2 
-  number  of  molecules,  which  fall  into  the  tube  through  section  S2, 
and  N2  -  a  number  of  molecules,  which  fall  into  the  tube  through 
section  S2,  ?12  -  probability  that  the  molecule,  which  passed  through 
section  S2,  will  reach  section  S2.  The  procedure  of  determination  Pl2 
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and  P,,  is  such.  In  section  Sx  are  developed  with  certain  probability 
density  the  coordinates  of  the  point  from  which  begins  its  motion  the 
molecule.  Then  with  the  probability  density,  which  corresponds  to  the 
function  of  the  distribution  of  molecules  in  the  undisturbed  flow,  is 
developed  the  velocity  vector  of  molecule.  If  particle  immediately 
flies  through  the  tube  before  section  S,,  then  from  section  Sx  is 
started  new  particle,  in  this  case  in  the  memory  of  electronic 
computer  is  memorized  one.  But  if  particle  falls  on  the  surface  of 
tube,  in  accordance  with  law  of  reflection  is  developed  the  velocity 
vector  of  particle,  which  is  reflected  from  the  surface.  When  as  a 
result  of  many  reflections  particle  nevertheless  achieves  section  S2, 
then  in  th^memory  of  machine  is  written/recorded  one,  but  if 
particle  achieves  again  section  Sw  '.nen  -  zero.  The  ratio  of  a 
number  of  particles  ,V„  of  those  flown  from  section  Sx  to  section  Sa, 
to  total  number  of  particles  N,  launched  in  section  Sx,  approaches 
Px  a  with 


11 


(3) 


Analogously  is  determined  value  PJt. 


The  use/application  of  the  Monte  Carlo  method  proved  to  be 
efficient  and  during  the  solution  of  the  problems  of  external  flow. 
For  the  first  time  this  method  was  proposed  and  used  in  work  [15],  in 
which  were  determined  the  aerodynamic  characteristics  of  the  internal 
surface  of  hemisphere,  which  moves  at  a  velocity,  arbitrary  according 
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to  the  value  and  the  direction.  During  the  calculation  by  the  Monte 
Carlo  method  automatically  drop  out  all  difficulties,  which  appear 
during  the  determination  of  the  aerodynamic  characteristics  of  the 
isolated/insulated  concave  bodies,  which  have  points  of  inflection, 
and  especially  the  groups  of  bodies. 

Formulation  of  the  problem.  Let  us  consider  the  overall  diagram 

of  the  solution  of  the  problem  of  free  molecular  flow  of  the  Monte 

Carlo  method  about.  Let  us  assume  is  assigned  the  body  of  arbitrary 

form  (or  the  group  of  bodies).  Let  us  assume  that  the  flow  of  gas 

about  the  body  everywhere  free  molecular,  are  assigned  the  function 

of  the  distribution  of  the  molecules  of  the  non-traveling  flow  and 

the  law  of  interaction  of  molecules  with  the  surface  in  any  form, 

even  in  the  form  of  table.  It  is  first  of  all  necessary  to  select  the 

control  surface  whose  form  depends  on  the  form  of  the  function  of  the 

distribution  of  the  molecules  of  the  incident  flow  and  on  the  shape 

of  body.  The  only  condition  for  its  selecting  is  the  following:  not 

one  molecule  of  the  incident  flow,  which  reaches  body,  must  omit 

control  surface.  It  is  natural  that  it  is  better  to  choose  this 

surface  in  such  a  way  that  it  would  take  as  less  than  the  molecules, 

which  fly  past  the  body.  Let  us  register  the  flow  value  of  the 

— »  — * 

molecules,  which  possess  the  speed  in  the  range  from  V  to  v+dv, 
through  the  element/cell  of  control  surface  d<t>: 

dN  —  fa>  V*dV  sin  {»  cos  tydydvdfi,  (4) 
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where  w  -  angle  between  the  normal  to  element/cell  d$ 

(standard/ normal  inside  the  control  surface)  and  the  direction  of 
speed  V; 

<p  -  orbital  angle. 

Page  57. 

We  choose  with  certain  probability  density  of  the  coordinate  of 
element/cell  of  d$  on  the  control  surface,  velocity  V  and  its 
direction,  given  by  angles  <//  and  <p.  We  determine,  did  fall  the  flow 
of  molecules  dN  on  body  surface  (or  one  of  the  group  of  bodies).  If 
the  flow  of  molecules  achieved  body  surface,  then  it  brought  into  the 
point  of  entry/incidence  on  the  surface  the  following  molecular 
signs/criteria:  the  flow  of  molecules  dN,  pulse  streams 
dNV dNVy,  dNVt  in  direction  x,  y,  z,  energy  f low  (1/2)1NV J  and  the 
"flow"  of  moment/torque  dN(i?x\r),  where”^  -  radius-vector  of  the  point 
indicated.  All  these  values  are  memorized  in  storage  cells  of 
computer(s).  If  on  the  body  surface  there  is  no  accumulation  of 
particles,  then  the  arrived  flow  of  molecules  dN  completely  will  be 
reflected  from  the  surface,  but  if  there  is  an  accumulation,  then 
will  be  reflected  only  the  part  of  flow  dN.  The  law  of  interaction  of 
molecules  with  the  body  surface  is  assigned,  i.e.,  is  assigned  the 
probability  of  reflection  in  this  direction  with  the  given  speed; 
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therefore  we  develop  with  certain  probability  density  the  speed  of 
the  flow  reflected  and  its  direction.  With  the  reflection  the  flow  of 
molecules  takes  away  from  the  surface  the  following  molecular 
signs/criteria:  the  flow  of  molecules  dN01,  pulse  streams 
<riVoi  Kto  <Wo.  V0yU  dN0jV ^  in  directions  x,  y,  z,  energy  flow  dN0  x  V20l/2 
and  the  "flow"  of  moment/torque  dN0 x  (RxVai).  All  these  values  are 
memorized  in  storage  cells  of  computer(s). 

Further  we  determine,  crossed  the  flow  reflected  body  surface  or 
not.  If  it  crossed,  then  we  find  the  coordinates  of  point  of 
intersection  and  send  into  storage  cells  of  computer(s)  the 
appropriate  molecular  signs/criteria  and  again  develop  reflection, 
etc.  This  process  is  continued  until  the  flow  of  molecules  » 
reflected  after  the  n  reflection  flies  out  from  the  control  volume. 
This  entire  process  is  called  one  testing.  Furthermore,  one  testing 
is  called  such  process  when  the  flow  of  molecules  dN  after 
"starting/launching"  from  the  control  surface  not  at  all  intersects 
body  surfaces  and  flies  past  it. 

After  conducting  k  of  tests  we  compute  local  aerodynamic 
characteristics,  summarizing  all  values  of  the  corresponding 
molecular  signs/criteria,  memorized  upon  contact  of  the  flow  of 
molecules  with  the  element/cell  of  body  surface  dFt,  and  subtracting 
the  sum  of  the  corresponding  molecular  signs/criteria,  memorized  at 
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the  moment  of  the  reflection  of  the  flow  of  molecules  from 
element/cell  dF} •  Difference  we  divide  into  a  number  of  drawings  k  and 
a  value  of  the  area  of  element/cell  dFt. 

During  the  calculation  of  total  aerodynamic  characteristics  it 
is  necessary  to  memorize  molecular  signs/criteria  at  the  moment  of 
the  first  entry/incidence  of  the  flow  of  molecules  dN  to  the  surface 
and  at  the  moment  of  its  latter/last  reflection,  after  which  the 
particles  will  fly  out  from  the  control  volume.  Differences  in  the 
sums  of  these  two  values  for  the  appropriate  molecular  signs/criteria 
are  divided  into  a  number  of  tests  k. 

Scanning/sweep  of  process  on  the  time,  i.e.,  successive  tracking 
the  particle  fluxes  dN,  is  possible  because  the  molecules  do  not 
interact  with  each  other  and  flow  pattern  can  be  then  represented  as 
the  imposition  of  many  flow  patterns  of  body  of  the  separate  flows  of 
molecules  dN.  It  is  possible  to  give  another  interpretation  to  the 
method  presented.  If  we  integrate  expression  (4)  with  respect  to  V, 

\l>,  <p  and  over  the  entire  control  surface,  then  this  integral  N  is 
equal  to  total  number  of  molecules,  which  fly  through  the  control 
surface.  Ratio  dN/N  is  nothing  else  but  the  probability  of  the 
flight/span  of  molecules  from  this  by  speed  V  through  this  point 
control  surface. 
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Page  58. 

Then,  after  making  the  procedure  described  above  not  for  the  flow  of 
molecules  dN,  but  for  one  molecule,  it  is  possible  to  compute  the 
average/mean  values  of  the  aerodynamic  character ist ics ,  which 
correspond  to  one  molecule,  which  flies  through  the  control  surface. 
For  the  determination  of  the  aerodynamic  characteristics  of  body, 
determined  by  the  total  flux  of  its  reaching  molecules,  it  is 
necessary  the  mentioned  above  aerodynamic  characteristics,  calculated 
for  one  molecule,  to  multiply  by  the  flow  value  of  molecules  N 
through  the  control  surface.  These  two  approaches  correspond  in  a 
sense  by  known  in  the  method  Monte  Carlo  to  the  diagrams  of  uniform 
and  essential  selection.  In  the  diagram  of  essential  selection  the 
accuracy  with  the  given  number  of  tests  is  above,  but  in  certain 
cases  the  use/application  of  this  diagram  causes  difficulty  due  to 
the  sharp  increase  in  the  number  of  computer  operations,  connected 
with  obtaining  of  random  numbers  with  the  assigned  law  of 
distribution  from  the  uniform  sequence. 

The  calculations  of  aerodynamic  characteristics  were  performed 
for  the  cases  of  diffuse  and  mirror  laws  of  reflection;  therefore 
will  be  given  below  basic  formulas  for  the  drawing  of  the  "start"  of 
particle  from  the  control  surface  and  the  formula  for  the  drawing  of 
the  report/event  of  the  reflection  of  particle  from  the  body 
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(analogous  formulas  can  be  obtained  also  for  the  arbitrary  law).  Let 

to  the  quiescent  body  attack  the  flow  of  molecules  the  function  of 

distribution  of  which  takes  the  following  form: 

/.  =  «- (2 .*7,)-**  exp  {-  [( VM -  ,)’  +(V,-Ua „)*+( V, - (/* ,)’l|  - 

Then  flow  dN  (y  axis  is  directed  along  the  normal  to  d$)  can  be 
registered  thus: 


dN-  Jbzgz.  b  + 1 ,  +  ®  (V  2  b  S.)]  j  <6  X 

X /( v,)/( l/x)/( Vt) dVtdVydV,i£  ; 

/W—J— 

V* 

/(^)-^  Vye-(yy 

/(I/j - ‘  e-<^-^cc)* ; 

V  « 

+(1  +  «(y2AS»)|J*. 

In  these  formulas 


IZ  ™  Oe  fl!  £4oy  - Uxb\  £/»  J  ™  £/oc  C 5 


a*  4-  £2  +  c2=  1; 

i/.  — _ El  .  5  _  ^  °  . 

K2/?roo  ’  "  K2/?r»  ’ 

'4° "  T  {*~  +  V*' *S».ll  +  ♦  (««  /2")i  J . 


Page  59. 


The  probability  of  the  "start"  of  molecule  from  the  point  of 
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control  surface  with  coordinates  d$  at  a  rate  of  v  is  equal  to 

1-  d-jr -nV,)f(V,)f(V,)dV,dV,dV,  . 

Let  us  take  the  known  procedure  of  obtaining  random  numbers  with 
the  assigned  law  of  distribution  from  pseudorandom  numbers  [16], 
evenly  distributed  in  the  interval  [0;  1].  For  Vx,  Vy,  V,  it  consists 
of  the  following.  Since  when  Vt  y  Vf  mtx  +  3  (where  Vym»i  —  value  V„ 
with  which  f(Vy)=fmn(vy))  function  f  is  in  effect  equal  to  zero,  let 
us  bound  the  range  of  change  Vy  by  interval  (0,  V0¥]: 

“V  V„  =  J-  (S.  4  +  Vi  +  sy  +  3. 

*-«<•  v'>~Tr-  that 

V°y  /«.x(^y)  _ 

f  A  S»*  +  v'2  +  625l  \ 

As  a  result  of  simple  conversions  we  obtain 

2fi(«-+VTFF3»)  +  3l  v; 


W«  +  V2  +  A25^ 

X  exp  {  (V2+4«S*  -  Woe)]*  -  +  /2  +  625|b)+3]  X 

xv;-bsj  ). 


Further  from  the  sequence  of  the  random  numbers  evenly 
distributed  in  interval  [0,  1]  we  choose  two  numbers  (£2  and  $2)  and 
we  check  inequality  £x<f*  ($2).  If  inequality  is  fulfilled,  then 

V,  -  *2  V9yVWZ  =  l,  ^  U9. 

For  Vx  and  Vt  the  procedure  is  analogous.  For  their  drawing  are  valid 


the  following  formulas 


5,  <  exp  [—  (6$4  —  3)2]; 
<  exp  [ —  (65«  —  3)2]. 
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With  executing  of  first  inequality  we  have 

Vx  =  .(gk±|af_z.3>.  Uoo, 

woo 

with  executing  of  second  inequality  we  have 

. (6E»  +  SoeC  —  3),, 

V  M - ‘  ?  v/ 00* 

woo 

Let  us  register  the  flow  of  the  diffuse  reflected  molecules  with 
the  function  of  velocity  distribution 

/=«.(*2^r.)-Mexp[-(2/?rj-1  V*j 

from  the  element/cell  of  body  surface  dF: 

dNm  -  nm  (* 2  RTJr**  exp  I  -  (2  RTJr'  V*]  V9  sin  <{»  cos  W  d<?  d  V  dF. 

Page  60. 


Total  number  of  molecules,  reflected  from  element/cell  dF,  is 


equal 


The  probability  of  reflection  with  this  speed  V  in  this 
direction  (9,  u)  is  equal  to 


where 


f(V)-2V3e~v'dV ;  /i8)- 2  sin  8  cos  8; 

f(u.)  —s  JL  •  v  *  v 

yuv  2*  ’  v  (2 RTJ1*  '  • 


Formulas  for  the  drawing  9,  v: 

sin  8  —  |*  = 

For  drawing  V  is  applicable  this  procedure.  We  choose  * j .  and 
5/+i  and  we  require  executing  'f  the  inequality 
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If  it  is  carried  out,  then  V*-ljku  where  kx  =  3.  For  the  law  of 
mirror  reflection  Vyatl=-— Vy0T9,  Vx  =  Vx0rp,  Vxatx  =  VxaTf. 

By  this  method  were  carried  out  the  calculations  of  the 
aerodynamic  characteristics  of  the  following  bodies:  the  hemisphere, 
converted  by  concave  side  to  the  flow  at  arbitrary  value  S«  [15]; 
cylinder  with  the  spherical  blunting  and  the  blades/vanes, 
arranged/located  perpendicularly  to  the  axis  of  cylinder;  the  body, 
which  consists  of  the  hemisphere  and  the  circle  and  arranged/located 
under  the  angle  of  attack;  the  body  of  that  consisting  of  the 
cylinder,  the  blades/vanes  and  the  cone;  the  body,  which  consists  of 
three  cones  with  the  spherical  blunting;  the  body,  which  consists  of 
the  cylinder  with  the  spherical  blunting  and  disk  [17].  For  an 
example  we  analyze  the  results  of  calculating  the  aerodynamic 
characteristics  of  hemisphere  and  body,  which  consists  of  the 
cylinder,  blades/vanes  and  cone. 

Hemisphere.  The  aerodynamic  characteristics  of  hemisphere  were 
determined  for  the  laws  of  diffuse  and  mirror  reflection.  Any 
aerodynamic  characteristic  can  be  registered  in  the  form 

n,  —  n  oc  -f-  ^n,0. 

Here  0^—  aerodynamic  characteristic  of  body  in  the  free 
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molecular  flow,  encompassing  the  aerodynamic  characteristic  of  the 
internal  part  of  the  surface  without  taking  into  account  of 
reflection  and  exterior  of  the  surface  taking  into  account 
reflection;  A'n<0  js  caused  by  interference. 

Fig.  1  gives  dependence  ^cy  and  r-,,j  (2*  +  sin  2a) 

(taking  into  account  only  the  first  reflection  of  molecules  from  the 
internal  surface)  on  a  when  Scc  =  oo.  Index  d  corresponds  to  diffuse 
reflection,  index  s  -  mirror;  where  7'»~  temperature  of 

the  molecules  reflected.  Fig.  2  depicts  the  dependence  of  addition  to 
the  value  of  the  "flow  of  moment/torque"  ^S^from 

.  ( '**"£*# ), 

Page  61. 

Fig.  3  depicts  dependence  on  »  ■  Let  us  note 

that  here  in  the  case  of  completely  diffuse  reflection  the 
interference  of  molecules  plays  no  role  and  value  a£*, where 

&Ed—  correction  only  due  to  the  first  reflection  of  molecules  from 
the  internal  surface  of  hemisphere.  This  fact  can  serve  as  the 
further  testing  of  the  accuracy  of  method. 


In  the  free  molecular  flow  for  the  computed  aerodynamic 
characteristics  in  the  case  of  the  laws  of  diffuse  and  mirror 
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reflection  taking  into  account  the  effect  of  the  exterior  of  the 
surface  are  valid  the  following  formulas: 

with  the  diffuse  reflection 

c*  —  +  ^  —  M,  =  M,t  +  2  £lMz4; 

cx  i  =  2  sin  a;  cs2  —  2  sin*  - :  -  =-1J  (*  -2a  -f  sin  2a)  (3S.)-’  ; 

2  ,  .  4 

Cyi —  2~r’  2  COS'  t;  .Mlt— - *-  1  COS*  a; 

ce=a‘cc  i  *  f  e  2  ~  cc  >  ™  sin  a; 

Ciri=(1~^) ""’(t  ~t)  = 

Lce  =  Ace - j^stna; 

with  the  mirror  reflection 

c x  “  t  4"  CM  2  ■+■  xj>  Cy  2  “f"  ^Cys* 
cx ,  =  2  sin  a;  cx ,  =  2  (2  +  sin  a)  sin*  ( -J-  — ;  r, ,  =  -  -5-  cos*  a. 

Here  coefficients  and  c>  are  related  to  0,5 */?* /»«„,  £/£,, 
moment/torque  is  calculated  relative  to  the  center  of  sphere  and  is 
related  to  0,5 */?o mn^Ul,,  coefficient  cE  is  related  to  0,5 itRlmn^LP^. 
Index  1  designates  the  internal  part  of  the  surface  of  hemisphere  in 
the  free  molecular  flow  without  taking  into  account  reflection,  inde 
2  -  the  exterior  of  the  surface  in  the  free  molecular  flow  taking 
into  account  reflection.  From  the  calculations  conducted  it  follows 
that  the  interference  has  an  effect  on  different  aerodynamic 
characteristics  differently.  The  most  sensitive  value  is  cr 
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Page  62. 

In  the  case  of  the  law  of  diffuse  reflection  the  interference  plays 
important  role  for  all  aerodynamic  characteristics,  if  only  body 
surface  has  sufficiently  high  temperature,  i.e., 

Body,  which  consists  of  the  cylinder,  blades/vanes  and  cone 
(Fig.  4).  Calculation  was  carried  out  at  5=45°,  RX=LX»L,.  All 
geometric  dimensions  are  related  to  R=Lx+L,/2  [1+(RX+L,/L, ) * ] . 
Formulas  for  any  aerodynamic  characteristic  are  written/recorded  in 
the  form  H,  —  Hl00-f AHJt  where  n(cc  is  determined  by  first  contact  of 
molecules  with  body;  Ant  —  correction,  caused  by  interference  and 
first  reflection.  At  a>90°  it  is  necessary  to  consider  the  effect  of 
the  end/face  of  cylinder  on  the  aerodynamic  characteristics.  In  this 
case  any  aerodynamic  characteristic  is  written/recorded  in  the  form 
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=  IT/  jo  —  An,  4-  niao  T  -j_  aii,  t> 

where  n*°°T  is  determined  by  the  molecules,  falling  on  the  end/face 
from  infinity,  AHJt-  by  the  molecules  reflected  from  the  end/face. 

In  the  case  of  diffuse  reflection 

Pcdix=  —  mrtoaU^  r.$  cos  a  sin  a;  APTJ  =  0; 

P<*>i  y  —  —  mrUa  Ulo  *Ri  cos2  a;  A PT  y  —  —  kRj  mri*,  Uh,  cos  a. 

In  the  case  of  the  mirror  reflection 

Ptaiy  +  bPi  ,  —  —  2ffIflao  Uh  nR  iCOS2*. 

Fig.  5-8  gives  the  dependences  of  dimensionless  aerodynamic 
characteristics  P^»tP//^,  E «,  Mtat  on  the  angle  of  attack  a. 
Corrections  to  the  aerodynamic  characteristics  in  the  case  of  diffuse 
reflection  are  given  for  Fig.  9-11,  in  the  case  of  mirror  reflection 
-  for  Fig.  12-14.  From  the  given  results  follows  the  conclusion  that 
the  interference  exerts  a  substantial  influence  on  cy  and  M ,  in  the 
case  of  mirror  and  diffuse  reflection  and  "hot"  wall. 


mm 
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Fig.  <?• 

Key:  (1).  taking  into  account  the  first  reflection.  (2).  taking  into 
account  interference. 

Page  64. 

From  the  numerous  calculations  conducted  it  is  possible  to  draw 
the  following  general/common/total  conclusions: 

-  if  molecules  are  reflected  from  the  body,  which  has  one 
significant  dimension,  with  speed,  there  is  much  lower  sjjeed  of  the 
molecules  of  the  incident  flow,  then  the  basic  contribution  to  the 
interference  is  determined  by  shading/blanketing  the  some  parts  of 
the  body  by  others; 

-  if  body  has  the  large  surfaces,  situated  in  parallel  to  the 
incident  flow,  then  interference  can  exert  a  substantial  influence  on 
the  value  of  some  aerodynamic  characteristics  even  in  the  case  of 
"cold"  body; 

-  interference  can  very  substantially  affect  aerodynamic 
characteristics  in  the  case  of  hot  wall,  i.e.,  when  molecules  are 
reflected  at  a  velocity,  equal  in  order  of  velocity  of  incident  flow. 
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Key:  (1).  without  talcing  into  account  interference.  (2).  taking  into 
account  interference. 


.  it —  (O 

,-C  MW>  £».) 


Fig.  14. 

Key:  (1).  without  taking  into  account  interference.  (2).  taking  into 
account  interference. 

Page  65. 

Calculation  of  flows,  close  to  free  molecular. 

With  the  decrease  of  the  local  Knudsen  number  the  collisions  of 
molecules  begin  to  play  the  significant  role  in  the  flow  pattern  of 
body  and  must  be  considered  them.  During  the  determination  of  the 
distribution  function  for  the  flow,  different  from  the  free 
molecular,  it  is  possible  to  use  the  formal  resolution  of  the 
solution  of  the  equation  of  Boltzmann  in  the  series/row  for  1/Kn. 
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Attempt  to  find  solution  in  this  form  taking  into  account  in  essence 
only  of  the  terms  of  order  1/Kn  they  were  undertaken  in  works 
[18-20].  However,  for  this  case  of  lining/calculation  they  were  very 
bulky  and  the  obtained  results  did  not  have  practical  interest.  In 
work  [21]  was  for  the  first  time  proposed  the  theory,  which  was 
subsequently  called  the  theory  of  the  first  intermolecular 
collisions.  Was  solved  the  problem  about  resistance  of  the  circular 
disk,  arrang*»d/located  normal  to  impinging  hyperthermal  flow  of 
molecules,  in  this  case  was  considered  the  correction  to  resistance 
in  the  free  molecular  flow  due  to  the  single  collisions  between  the 
molecules  reflected  and  the  molecules  of  the  incident  flow.  Within 
the  framework  of  the  theory  of  the  first  intermolecular  collisions 
were  evaluated  the  aerodynamic  characteristics  of  infinite  band  [22], 
infinite  cylinder  [23]  and  sphere  [24]  (in  the  linings/calculations 
of  latter/last  work  was  allowed  the  error,  which  led  to  the  fact  that 
the  obtained  allowance  of  the  value  of  drag  coefficient  in  the  free 
flow  was  overstated  approximately/exemplarily  doubly). 


Several  on  another  path  went  in  his  studies  of  flows,  close  to 
the  free  molecular  ones,  Willis  [25].  He  used  instead  of  the  equation 
of  Boltzmann  the  modified  model  equation  of  Crooke  and  he  computed 
his  first  approximation,  after  making  for  the  zero  approximation  free 
molecular  decision.  The  formulation  of  the  problems  in  this  form  is 
qualitatively  close  to  the  formulation  of  the  problem  in  the  theory 
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of  the  first  intermolecular  collisions.  To  a  deficiency/lack  in  this 
method  it  is  necessary  to  relate  the  use  of  model  equation  of  Crooke 
whose  suitability  is  doubtful  for  the  flows,  distant  from  the 
equilibrium  ones.  Furthermore,  and  in  the  case  of  model  equation 
linings/calculations  are  bulky,  so  that  calculations  were  carried  out 
only  for  the  sphere,  cylinder  [26]  and  flow  of  Couette  [25].  In  work 
[27]  for  calculating  the  aerodynamic  characteristics  of  sphere  was 
used  the  method  of  resolving  the  solution  of  the  model  equation  of 
Crooke  in  the  series/row  for  1/Kn  and  was  found  the  first  term  of 
series/row;  the  obtained  results  qualitatively  coincide  with  results 
[26].  In  works  [28,  29]  was  used  the  method  of  successive 
approximations  for  the  integral  kinetic  equation,  for  whict^was 
designed  the  first  approximation.  In  particular,  in  [28]  was 
determined  stagnation  pressure  of  the  mirror  reflecting  sphere. 

The  most  complete  analysis  of  the  flows,  close  to  the  free 
molecular  ones,  was  carried  out  in  the  work  of  M.  N.  Kogan  [30,  31]. 
In  particular,  in  them  it  was  shown  that  in  the  theory  of  the  first 
intermolecular  collisions  for  the  convex  finite  bodies  in  the  cases, 
which  are  of  practical  interest,  it  is  necessary  to  seek  corrections 
and  to  the  free  molecular  value  of  aerodynamic  characteristics  only 
due  to  the  single  collisions  between  the  molecules,  reflected  from 
the  body,  and  by  the  molecules  of  the  incident  flow.  These  works 
examine  also  a  question  about  the  limits  of  the  applicability  of  the 
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theory  of  the  first  collisions,  are  qualitatively  analyzed  the  basic 
types  of  the  flows,  which  have  practical  interest  both  for  the  "cold 
ones"  and  for  the  "hot"  bodies  (i.e.  for  the  actual  conditions  and 
the  test  conditions  in  the  indraft  wind  tunnels),  are  indicated  the 
similarity  parameters  and  the  methods  of  the  recalculation  of 
obtained  theoretical  data  to  the  actual  conditions  and  the  conditions 
for  experiment  in  the  wind  tunnels. 

Page  66. 

M.  N.  Kogan  also  showed  that  if  we  use  the  formal  method  of  solving 
the  equation  of  Boltzmann,  expanding  the  distribution  function  in 
series/row  according  to  1/Kn,  then  the  solution,  found  with  an 
accuracy  down  to  the  terms  1/Kn,  corresponds  to  the  theory  of  the 
first  intermolecular  collisions.  Use/application  of  the  Monte  Carlo 
method  in  the  theory  of  the  first  intermolecular  collisions  made 
possible  to  create  the  universal  calculation  method,  which  makes  it 
possible  to  determine  the  aerodynamic  characteristics  of  convex 
bodies  under  arbitrary  law  of  reflection  of  molecules  from  the  body 
and  arbitrary  law  of  interaction  between  the  molecules. 

Formulation  of  the  problem.  Since  during  the  motion  of  flight 
vehicle  at  the  high  altitudes  its  speed,  as  a  rule,  is  much  more  than 
the  average/mean  thermal  velocity  of  molecules,  then  we  during  the 
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solution  of  tr.e  proc.ems  zi 
of  the  molecules  o:  tre  . '  ze- 
macroscopic  speed  and  -“a- 

hyperthermal  flow  v.:*.  a  sp- 
theory  of  the  first  i-.rermu.- 
examine  only  the  collis.cn* 
body,  and  by  the  molecules  cf  *  re 


c.sreqara  the  tnermal  velocity 
.  '  comparison  with  its 
•  re  oody  attacks  the 
Ac  •  .*  *as  said,  that  in  the 

.  .  ns  . t  .s  necessary  to 

-*.-les,  reflected  from  the 

flow. 


Let  us  register  the  number  of  co*lis.o  which  occur  in  certain 

volume  element  between  the  molecules,  reflected  from  the  element/cell 

-rr  -y 

of  surface  dF  with  the  speeds  in  the  range  from  v,  to  v,+dv,  and  the 
molecules  of  the  incident  flow  per  unit  time  in  the  element  of  volume 


dr: 


dN  =/» ( Vt )/. ( Vt) g^-'  *(2A)~‘  dbd»dVtdV \  dr. 


If  this  expression  is  integrated  over  entire  physical  and 
high-speed/high-velocity  space,  then  we  will  obtain  total  number  of 
collisions  N  per  unit  of  time.  Ratio  dN/N  is  nothing  else  but  the 
probability  of  this  collision.  The  collisions  of  molecules  in  a 
two-fold  manner  affect  the  aerodynamic  characteristics  of  the  body: 
on  one  hand,  due  to  the  collisions  to  the  body  comes  the  further  flow 
of  molecular  signs/criteria,  on  the  other  hand,  due  to  the  collisions 
to  the  body  does  not  come  certain  flow  of  the  molecular 
signs/criteria  which  carried  from  infinity  of  molecule,  which  clashed 
with  the  molecules,  reflected  from  the  body.  Consequently,  for  any 


aerodynamic  characteristic  of  body  taking  into  account  the  first 
collisions  is  valid  the  recording:  n,  —  n<e.1.4-IIi+  —  ni- .  where  Dic.m-L 
value  of  aerodynamic  characteristic  in  the  free  molecular  flow;  ni+  — 
addition,  caused  by  collisions;  n<-—  loss  due  to  the  collisions. 
Developing  randomly  the  sufficiently  large  number  of  collisions  K  in 
the  space,  it  is  possible  to  determine  the  local  and  total  values  of 
the  aerodynamic  characteristics  of  body,  averaging  additions  ELI+  and 
Hi-  according  to  the  number  of  collisions  K.  If  molecules  are 
reflected  from  the  body  with  the  average  speed  which  is  much  lower 
than  the  speed  of  the  molecules  of  the  incident  flow,  then  at  the 
moment  of  colliding  the  molecules  it  is  possible  to  consider  that  the 
speed  of  the  molecule  reflected  is  equal  to  zero;  this  assumption 
makes  it  possible  to  find  in  the  form  of  quadratures  the  function  of 
the  distribution  of  molecules  on  the  body.  In  this  case  the  Monte 
Carlo  method  is  applied  to  the  calculation  of  the  repeated  integrals, 
through  which  are  written/recorded  the  aerodynamic  characteristics  of 
body . 


Page  67. 


Let  us  illustrate  based  on  the  example  of  the  flow  around 
sphere,  how  occurs  the  simulation  of  flow  pattern  with  the  help  of 
the  Monte  Carlo  method.  Let  us  introduce  the  following  assumptions: 
1)  the  reflection  of  molecules  from  the  body  occurs  according  to  the 
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diffuse  law,  2)  molecules  are  elastic  spheres  by  diameter  a.  Let  us 


introduce  dimensionless  parameters  5, 


U  00  J  J,  1 

r r\ — 3  nd  Klloo  — .  —  T-  — 

V  2  /JT,,  2  \r2  tea *  rtx.Ro 


where  T»— -  temperature  of  the  reflected  according  to  the  Maxwellian 


law  molecules 


R  -  gas  constant;  rr a1  -  collision  cross  section  for  the  solid 
spheres;  density  at  infinity. 

Let  us  register  the  number  of  collisions,  which  occur  at  point 
A( x ,  y,  2 )  (Fig.  15)  between  the  molecules,  reflected  from 
element/cell  dF  of  the  surface  of  sphere  with  the  speeds  in  the  range 

-5»  ^ 

v,  v+dvi,  and  by  the  molecules  of  the  incident  flow  per  unit  time  in 
the  element  of  volume  pJdp  sin  ^  d^  d<p: 

dN  —tuianm  (2*  RTJr3'1  exp  ^  —  2 cos  ^  sln  *  ^  °2  850  v  cos  v  x 

X  rfp  exp  ^  j  dp  dF. 

Molecules  were  reflected  from  the  element/cell  of  surface  dF=RJ, 
sin  0d  0d  0  in  the  direction  ,  <p.  At  the  moment  of  collision  the 
direction  of  center  line  is  assigned  by  angles  v ,  u  relative  to  the 
direction  of  relative  speed  <3?j  —  y  —  £/» .  Collision  occurred  at  a 
distance  p  from  element/cell  dF.  The  molecules  reflected  possessed 
the  speeds  in  the  interval  (v;  v+dv) ,  and  their  quantity  is 
proportional  to  exp(  -p/X  u),  where  X  ^  the  mean  free 

path  of  the  molecules  reflected  on  the  molecules  of  the  incident 
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flow.  From  the  law  of  conservation  of  a  number  of  particles  on  the 
surface  of  sphere  it  follows  that  n.  =  2/t,0  cosg 

The  collision  frequency  let  us  register  in  the  form 

dN=  itina,  l/./(9)/(«/W/(p)/(?)/(P)  X 

.  X/(V)/(fl)^0  dVdtydvdpdidpdp, 

where 

/(8)  =  2  sin  9  cos  8; 

/(V0  =  2(2/?T.)-*exp(-2^j  V3; 

/(<b)  —  2  sin  •{»  cos  iji; 

/(*)  —  2  Sin  v COS r,  /( p)  —  jL  exp  J^—  ; 

/(?)=/(fO-/(P)  =  jj; 

moreover  for  all  variable/alternating 
[f(xi)dxt  =  1. 
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Fig.  15. 

Page  68.  • 

Then  the  probability  of  this  collision  exists 
-nSr  -/(•)/(?) /W/(p)/(?)/(P)/(j»)/( VO <*8 
One  drawing  let  us  name  the  collision,  characterized  by  the 
random  sampling  of  eight  parameters,  moreover  each  parameter  was 
chosen  randomly  with  an  appropriate  probability  density  of  /(*/)•  After 
the  collision  was  played,  were  determined  the  trajectories  of 
molecules  after  collision,  was  fixed/recorded  the  fact  of  falling  of 
molecule  into  this  point  of  surface  and  were  memorized  the 
corresponding  values  of  the  molecular  signs/criteria,  yielded  by 
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molecule  to  the  body.  Simultaneously  it  was  checked,  did  not  occur 
collisions  at  the  point  of  space  which  is  projected/designed  along 
the  direction  of  speed  to  the  body  surface;  if  this  is  so,  then  at 
this  point  of  body  surface  we  fix/record  the  report/event  of  the  loss 
of  the  molecular  signs/criteria  which  carried  the  molecule  from 
infinity.  The  parameters  of  each  collision  are  developed  according  to 


V  -  VWT9  V  —  In  Mi+i;  sin  *  = 
sinv  — P=  (  —  2  /2~  KncoJ^-'l  InSi-n; 


the  following  formulas: 


9  =  2wSi+s;  n  =  2i3i+6;  ^  —  2s5i+7; 

sin  8  = 

where  &<+«  —  evenly  distributed  random  numbers  [0-1]. 

After  playing  collision,  we  determine  the  speeds  of  molecules 
after  collisions  Vt  and  Km-i-  Each  molecule  after  collision  has  the 
"weight",  equal  to 

mn*, £4o  ~Rl  —  for  the  flow  of  molecules; 

mn^Uoa-zRo  —  for  the  flow  of  normal  impulse/momentum/pulse; 

mn^U„-RoUt/,/+i~  for  the  pulse  stream  in  direction  £/»; 


mn* 0  2  for  the  energy  flow. 
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If  collision  occurred  above  any  point  of  body  surface,  then  at 
this  point  of  surface  occurred  the  loss  of  the  following  molecular 
signs/criteria: 

nuiaUaxfto  —  for  the  flow  of  molecules; 

'n«oof/oo*/?o(f/ao«)  —  for  the  flow  of  normal  impulse/momentum/pulse; 

/nn« «/?o —  for  the  pulse  stream  in  direction 

(ji 

mrua t/oo *Rq  —  for  the  energy  flow. 

If  it  is  necessary  to  know  the  distribution  of  the  computed 
flows  from  the  body  surface,  it  is  necessary  body  surface  to 
decompose  into  the  elements/cells  (value  of  which  it  is  determined 
from  the  considerations  of  accuracy),  to  count  the  sums  of  the  flows 
of  the  molecular  signs/criteria  interesting  us  and  to  relate  them  to 
the  area  of  the  element/cell  of  surface  and  to  a  number  of  drawings 
K.  For  the  total  characteristics  necessary  to  count  the  sums  of  the 
corresponding  molecular  signs/criteria,  which  arrive  at  entire  body 
surface  and  lost  as  a  result  of  collisions  above  the  entire  body 
surface  relate  them  to  K. 

Thus,  we  compute  H+*  and  H— where  index  k  indicates  the 
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appropriate  molecular  sign/criterion. 

Page  69. 

A  number  of  drawings  K  is  determined  by  the  assigned  accuracy  of 
calculation  according  to  known  formulas  [16],  moreover  dispersion  is 
determined  in  the  process  of  calculation.  By  the  method  proposed  were 
determined  the  aerodynamic  characteristics  of  the  whole  class  of 
bodies. 


1.  Sphere.  (Fig.  16).  Aerodynamic  characteristics  of  the  sphere 

■5 

were  computed  for  arbitrary  value  S.  when  [32].  It  was  assumed 

that  the  reflection  of  molecules  from  the  surface  diffuse  and  the 
molecules  interact:  a)  as  elastic  spheres;  b)  as  pseudo-Maxwell ian 
molecules.  Calculations  showed  that  value  c*  in  the  mode/conditions 
of  the  first  intermolecular  collisions  is  reduced  in  comparison  with 
its  free  molecular  value  approximately/exemplarily  to  15%  within  the 
framework  of  the  applicability  of  theory.  The  results  of  detailed 
calculations  are  given  in  works  [32-34  ]. 

Cone  (Fig.  16  and  17).  In  this  case  the  effect  of  the  first 
intermolecular  collisions  on  different  aerodynamic  characteristics  is 
different.  Allowance  to  resistance,  energy  flow  and  to  the  value  of 
the  pitching  moment  of  "cold"  cone  when  Sco$<l  (with  exception  of  the 
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thick  cones  when  8>60°,  where  6  -  half-angle  of  solution/opening)  is 
negative  and  does  not  exceed  10%  of  value  of  these  aerodynamic 
characteristics  in  the  free  molecular  flow  in  the  region  of  the 
applicability  of  the  theory  of  the  first  collisions.  For  example, 
correction  due  to  the  collisions  for  the  cone  with  the  significant 
dimension  of  h=l  m  when  6=  15°,  *=»  72°,  —  10  for  the  height/altitude  of 

120  km  is  equal  to  approximately/exemplar ily  7%.  Allowance  as  a 
result  of  the  collisions  to  the  free  molecular  value  of  the  lift  of 
"cold"  cone  in  the  dependence  on  the  combination  of  the  half-apex 
angle  of  the  cone  and  angle  of  attack  can  be  both  positive  and  by 
negative  and  on  the  average  are  approximately  50%  of  c,t».  In  the 
"resonance"  case  when  a*jr/2-0,  value  can  be  several  times  more 
than  (in  the  given  higher  example  at  the  height/altitude  of  120 

km  relative  value  of  correction  Ac,  it  is  equal  to 

approximately/exemplar  ily  150%  e*<-»  and  only  at  the  height/altitude 
of  230  km  it  composes  10%). 
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Fig.  16.  1  -  circle;  2  -  square;  3  -  rectangle;  4  -  cone,  0=60°;  5  - 
sphere;  6  -  cone,  0»45°;  7  -  cone,  5*30°;  8  -  circle  (along  the 
flow) . 

Page  70. 

For  the  thick  "cold"  cones  the  allowance  due  to  the  collisions  to 
free  molecular  to  the  value  of  pitching  moment  can  be  both  the 
positive  and  negative,  and  its  value  can  be  much  more  than  the  free 
molecular  value  of  pitching  moment.  In  the  case  of  very  slender  cones 
when  Sas9>l  the  value  of  resistance  taking  into  account  the  first 
collisions  is  more  than  free  molecular  value.  From  the  given 
calculations  it  is  possible  to  draw  the  following  basic  conclusions; 
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-  the  boundary  of  free  molecular  courses  for  different 
aerodynamic  characteristics  of  cone  lies/rests  on  different 
height/altitude  and  there  cannot  be  determined  only  on  the  basis  of 
the  Knudsen  number;  it  in  many  respects  depends  on  the  geometry  of 
task,  laws  of  interaction  of  molecules  with  each  other  and  with  the 
surface  of  cone; 

-  ignorance  of  the  laws  of  interaction  of  molecules  with  the 
surface  (accommodation  coefficients)  to  a  high  degree  affects  the 
determination  of  lift  and  pitching  moment  of  cone  and  to  the 
considerably  smaller  degree  -  the  determination  of  the  value  of  ^ 
resistance  and  heat  flux. 

Plate.  The  results  of  detailed  calculations  are  published  in 
works  [33,  35-40].  Plate  is  most  "convenient"  body  for  obtaining  the 
qualitative  laws  governing  the  flow  around  bodies  by  the  strongly 
rarefied  gas.  In  work  [30]  based  on  the  example  of  plate  were 
determined  different  types  of  courses,  realized  in  the  case  of  the 
flow  around  the  plate,  arranged/located  at  different  angles  of 
attack,  under  different  laws  of  interaction  of  molecules  with  each 
other  and  with  the  surface  of  plate. 


The  calculations  conducted  completely  confirmed  qualitative 
results  [30].  If  plate  is  arranged/located  normal  to  flow  [36,  39, 
40],  then  its  resistance  and  energy  flow  to  it  in  the 
mode/conditions,  close  to  the  free  molecular,  are  reduced  within  the 
framework  of  the  applicability  of  theory  approximately/exemplar ily  to 
10%.  In  the  case  of  the  plate,  arranged/located  in  parallel  to  flow 
[33,  35,  37,  38],  possibly  several  types  of  flows  depending  on  the 
"type"  of  molecules,  law  of  their  interaction  with  the  surface  and 
values  of  numbers  S«,,  S9  and  Km*. Thus,  in  the  "cold"  case  when 
5w<Kn<»  resistance  of  plate  and  energy  flow  in  the  flow  conditions, 
close  to  the  free  molecular,  are  more  than  the  free  molecular  value 
of  maximum  by  ~10%. 


Page  71. 


But  if  are  realized  the  conditions  of  the  so-called  molecular 
boundary  layer  when  I  <  Kni,  then  within  the  framework  of  the 

applicability  of  the  theory  of  the  first  intermolecular  collisions 
value  c%  and  energy  flow  can  be  many  times  of  more  than  their  value 
in  the  free  molecular  flow.  Even  to  the  larger  degree  grows  pressure 
on  the  surface  of  plate.  For  the  plate,  just  as  for  the  thin  pointed 
bodies,  high  value  has  a  precise  knowledge  of  law  of  reflection  of 
molecules  from  the  surface,  since  depending  on  the  law  of  interaction 
of  molecules  with  the  surface  of  the  value  of  aerodynamic 
characteristics  they  can  increase  or  be  reduced  by  an  order.  Fig.  16 
gives  some  results  of  calculations  for  the  plates  of  various  forms. 
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Skimmer  (Fig.  18).  Calculations  were  carried  out  for  the  purpose 

of  the  determination  of  flow  disturbances  at  the  entrance  into  the 

skimmer,  the  caused  by  collisions  molecules  of  hypersonic 

incident  flow  with  the  molecules,  reflected  from  surface  [34],  In  the 

investigated  modes/conditions  the  mean  free  path  of  molecules  in  the 

incident  flow  was  much  more  than  the  diameter  of  the  entrance, 

molecule  they  were  assumed  elastic  spheres,  S.=  10  and  1.57.  From  the 

obtained  results  it  is  evident  that  even  when  Kn,=40  and  S  =10 

a  w 

density  at  the  entrance  increases  by  10%  in  comparison  with  its  free 
molecular  value.  Axial  inlet  velocity  falls,  average/mean 
longitudinal  quadratic  speed  is  reduced,  all  this  together  undertaken 
speaks,  that  at  the  entrance  begin  to  predominate  the  oblique 
collisions.  With  an  increase  in  the  wall  temperature  the  effect  of 
collisions  is  reduced.  Consequently,  cooling  skimmer  can  contribute 
to  the  decrease  of  the  distortion  of  flow  only  in  the  case  of  the 
complete  freezing  of  gas  on  its  surface.  The  calculations  conducted 
show  that  even  in  the  case  when  Kn^^l  and  knD>l,  intermolecular 
collisions  significantly  distort  the  picture  of  flow  after  the 
skimmer. 
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Fig.  18. 
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STUDY  OF  INTERNAL  AND  EXTERNAL  FREE  MOLECULAR  FLOWS  ABOUT  AN 
ARBITRARY  GROUP  OF  COMPLEX  BODIES. 

M.  A.  Zakirov. 

Are  given  formulas  and  universal  program  comprised  on  them  for 
calculating  of  aerodynamic  coefficients,  local  flows  through  the 
surfaces  and  local  parameters  of  gas  in  the  field  of  flow  at  the  flow 
of  the  free  molecular  flow  about  the  group  of  complex  bodies. 
Calculations  can  be  carried  out  under  different  laws  of  interaction 
of  molecules  with  the  wall.  Body  surface  is  assigned  analytically. 

The  velocity  vector  of  the  incident  flow  is  arbitrary  in  the  value 
and  in  the  direction.  Program  makes  it  possible  to  carry  out  the 
flow-field  analyses  of  bodies  by  the  flow  of  light/world  and  by 
hypersonic  inviscid  flow  (according  to  Newton's  theory). 
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ADOPTED  DESIGNATIONS. 

v,  Vhui,  h~v2 ,T,n,  m—  average/mean  and  relative  gas  velocity,  the  most 
probable  speed,  temperature,  density  and  the  mass  of  particles; 

Xxh  Yyj,  Vw„—  coordinate  systems,  connected  with  entire  body,  with  the 
separate  surface  and  with  the  surface  element  on  surface  of (i ,  j, 
n-1,  2,  3); 

unit  matrix  (j,  i=l,  2,  3); 

e±i.  crl,  m±l.  mri— coef f  icients  of  total  aerodynamic  forces  and 
moments/torques ; 

£± ,  Er—  coef  f  icients  of  total  energy  flows  to  the  body; 

E*~  flux  coefficient  of  the  energy  through  the  control  surface; 

iP  *  .  Pr\  | 

,  \  t.~  »  J  —  coefficients  of  the  average/mean  local  particle 

'T±  .  VJ.  >  er)  ) 

fluxes,  normal  and  tangential  impulses/momenta/pulses  and  energy 


through  the  surface; 
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v_l—  probability  of  the  flight/span  of  molecule; 


P\—  static  pressure; 


«x»  VW  T»  •  1  particle  density,  the  average  speed,  temperature, 
Pij  EJ  <M  •  <M  J 

the  tensors  of  the  flux  of  momentum  and  stresses/voltages,  the  flows 
of  complete  and  thermal  energy  along  the  field  of  flow  of  gas; 


a*—  particle  flux; 

M[X\,  mathematical  expectation  and  dispersion; 

N9(Nrk)—  number  of  particles,  played  through  entire  (through  the 
k-th)  face  of  control  surface. 
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Indices. 


<*»  -  undisturbed  flow; 


w  -  condition  on  the  wall; 


k  -  number  of  the  face  of  control  surface; 
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v  -  v- th  surface; 

u,  7  -  local  flows  through  the  surface; 

X  -  local  parameters  on  the  field  of  flow  in  the  vicinity  of  formal 
surface  sx-, 

n  -  coefficients  in  Newton's  theory; 


f  -  flow  of  light/world  (photons). 
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Introduct ion. 

The  tasks  of  the  study  of  free  molecular  flows  appear  in 
connection  with  the  flights  of  bodies  in  the  free  molecular  flow  and 
the  calculations  of  flows  in  the  vacuum  devices. 

The  calculations  of  free  molecular  flows  are  performed  through 
two  directions:  in  the  first  case  solve  the  complicated  integral 
equations  of  the  type  of  Fredholm  the  second  kind  for  the  function  of 
the  distribution  of  particles  [1-6]  reflected;  in  the  second  case 
free  molecular  flow,  they  simulate  by  the  method  of  the  Monte  Carlo 
[7-12]. 

The  solution  of  integral  equation  is  obtained  for  the  case  of 
flow  of  the  hypersonic  free  molecular  flow  about  of  the  concave 
sphere  and  cylinder  with  the  diffuse  reflection  of  particles  from 
surface  [2,  13,  14],  Numerical  calculations  are  carried  out  during 
the  free  molecular  flow  around  the  infinite  cylinder  [15]  and  wedge 
[16];  cone,  hemisphere  and  cylinder  with  wings  [5,  8]?  flat  ducts  and 
lattices  [17]. 

In  the  tasks  about  the  internal  flows,  the  first  settings  and 
solutions  of  which  were  given  in  the  works  of  Knudsen,  Smoluxovsky 
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and  Clausing  (see  vast  bibliography  in  the  works  [1,  5,  6]),  are 
determined  the  probability  of  the  flight/span  of  molecules  and  the 
density  of  gas  during  the  flow  through  the  rectangular  cylindrical, 
the  conical  and  wedge  channels  [17-20].  Is  determined  also  gas 
density  in  the  vicinity  of  sphere  and  cone  [21-23]. 

In  the  present  work  is  given  the  calculation  procedure  by  the 
method  of  the  Monte  Carlo  of  the  aerodynamic  coefficients  of  separate 
surfaces  and  entire  body,  and  also  coefficients  of  local 
stresses/voltages  in  the  body  surface  and  local  parameters  of  gas 
(density,  temperature,  speed,  heat  flux,  etc.)  on  the  field  of 
internal  and  external  free  molecular  flows.  Calculations  can  be 
carried  out  for  arbitrary  combination  of  the  complex  bodies  whose 
surfaces  are  assigned  analytically.  The  velocity  vector  of  the 
incident  flow  of  gas  can  be  arbitrary  in  the  value  and  the  direction. 
For  gaming  the  reflection  of  molecules  it  is  possible  to  apply 
different  models.  Employing  this  procedure  it  is  possible  to  lead 
calculations  for  the  bodies  in  the  flow  of  light/world  and  according 
to  Newton's  theory. 
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Chapter  1. 

UNIVERSAL  ALGORITHM  OF  THE  STUDY  OF  FREE  MOLECULAR  FLOWS. 

The  overall  diagram  of  the  solution  of  problem  is  the  following. 
The  geometry  of  the  arbitrarily  arranged/located  bodies  in  the  space 
is  assigned  with  the  help  of  set  of  parameters  of  surfaces.  In 
accordance  with  the  function  of  particle  distribution  of  the 
undisturbed  flow  is  produced  the  drawing  of  the  random  parameters  of 
particles  at  the  "start"  from  the  control  surface,  which  surrounds 
body.  Are  developed  consecutive  collisions  and  reflections  (in 
accordance  with  the  assigned  density)  of  particles  with  the  body 
surfaces  up  to  their  escape  from  the  control  volume.  At  the  moment  of 
the  intersection  with  the  molecule  of  body  surface  are  summarized  all 
necessary  flows,  throughout  which  then  are  calculated  the 
mathematical  expectations  of  the  corresponding  flows  per  unit  time. 

Page  74. 

1.  On  the  analytical  method  of  the  assignment  to  surface  of  complex 
body. 


Let  us  consider  the  complex  body  whose  entire  surface  can  be 
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represented  by  the  set  of  the  arbitrary  parts  of  the  surfaces  of  the 
second  order.  The  numbered  parts  of  the  surfaces  we  will  call  thus: 
the  first  surface,  the  second  surface,  etc.  Let  us  introduce  the 
system  of  coordinates  Xxt  (/=  1.  2,  3)  connected  with  the  body. 

All  surfaces  are  determined  by  assignment  for  each  of  them  the 

ordered  array  from  the  coefficients  of  the  equations  of  the  second 

order  .  ,  ~ .  ,  .  -  , 

bJmyjym  +  2bJAyj  +  bu=Q  (/,  m  -  1,  2,  3),  (1.1) 

written  for  the  convenience  relative  to  the  arbitrary  system  of 
coordinates  YyJt  oriented  relative  to  system  Xxi  by  direction 
cosines  between  axes  x,  and  yt  and  by  coordinates  xot  of  point 
Y.  In  formula  (1.1)  bJm—  tensor  (affine  tensor);  —  vector;  b44 

-  scalar. 

Analytical  method  includes  also  the  assignment  to  region  D  of 
changing  the  variable/alternating  During  the  assignment  to  D 
region  use  Cartesian,  polar,  cylindrical  and  spherical  coordinates. 
Then  almost  for  all  surfaces,  presenting  interest  from  the  point  of 
view  for  the  theory  and  practice,  the  form  D  region  becomes  simplest 

-  rectangular.  For  the  assignment  to  D  region  are  used:  in  the 
Cartesian  system  of  coordinates  (yj)  the  variable/alternating  yt,  y2, 
in  polar  ( <p ,  r,  y})  -  <p,  r,  in  the  cylindrical  (<p,  r,  y3)  -  <p,  y3  and 
into  the  spherical  («,  9,  r)  -  <p,  9 .  Let  us  agree  on  subsequently  to 
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call  and  cp  external  variable/alternating,  and  y2,  r,  y3  and  Q  - 
internal  variable/alternating. 

Let  us  do  some  notes.  In  the  Cartesian  and  polar  coordinate 
systems  D  region  does  not  set  limitations  on  values  of  y3;  therefore 
for  achievement  of  uniqueness  is  considered  only  the  part  of  the 
surface,  which  falls  into  half-space  y3>0.  By  the  appropriate 
selection  of  the  coefficients  of  equation  (1.1)  it  is  possible  to 
achieve  entry/incidence  into  half-space  y3>0  of  different  parts  or 
entire  surface.  In  the  case  of  conical  surface  it  is  necessary  to 
throw/reject  a  small  vicinity  about  the  apex/vertex,  since  the 
direction  of  standard/normal  at  this  point  is  not  determined. 

Thus,  each  surface  is  assigned  by  set  of  parameters  bjm<  -bn, 
b  *  * '  an *  Pi,  Pir  7 1 ,  7  j  and  Here  (plf  £,),  (ylr  y2)  -  the 

limits  of  the  rectangular  region  D  on  by  external  and  internal 
variable/alternating;  —  sign/criterion,  which  designates  the 
coordinate  system,  in  which  is  assigned  D  region. 

The  following  space  is  the  introduction  of  the  typical  classes 
of  surfaces  and  complex  bodies.  For  example,  if  are  examined  only 
axisymmetr ic  surfaces  (circle,  cylinder,  cone,  spherical  segment), 
then  a  number  of  assigned  parameters  about  the  surfaces  substantially 
is  reduced.  Complicated  surfaces  also  can  be  assigned  by  a  smaller 
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number  of  parameters;  for  example,  the  surface  of  parallelepiped  is 
determined  by  the  assignment  of  the  coordinates  of  three 
apexes/vertexes  and  height/altitude. 

Particle  trajectories  are  calculated  relative  to  axes  Xxit 
therefore  the  equations  of  all  surfaces  are  converted  from  one  axes 
Yyx  to  the  next  Xxt.  This  is  realized  into  two  stages.  First  is 
produced  the  rotation  of  axes  Yy,  so  that  new  axes  Yy\  \\  Xxt  and 
equation  (1.1)  takes  the  form 

I'rt/iJ'*  +  2 =  0;  =  <*ji  amk  bfm\  (i.2> 

ti  *  -  aji  bj «;  yj  -  a}lyx\  ym  -  amky„ ;  i,  k  -  1,  2,  3. 

Page  76. 

Then  is  realized  the  parallel  shift/shear  of  axes  Yy\  before 
the  coincidence  with  axes  Xxt 

yi—xt  x9l,  (1-3) 

and  the  equation  of  surface  takes  the  form 

aikx -i- 2aMx, -f- a44  =  Ch  alk  =  '!»«;  1  (14) 

4  “  4  —  '}’/*  xo  »’>  *“  X0  i  (1 XOk  2  4)  + 

2.  Selection  of  the  form  of  control  surface. 

On  control  surface  (S0)  we  will  produce  gaming  the  random  values 


DOC  =  83067303 


PAGE 


E 


of  the  components  of  vector  of  the  particle  speed  and  coordinates  of 
the  start  of  particle.  In  this  case  the  value  of  the  probability  of 
contact  with  the  body  surface  of  the  particle,  which  starts  with  S„, 

„  A/»  (2.1) 

Pa~  Np 

must  be  closest  possible  to  one.  Here  —  number  of  particles 
fallen  on  body  of  all  that  played  on  S». 

If  we  as  Sa  take  sphere,  then  value  pn  will  be  small  for  the 
bodies,  extended  on  one  or  two  measurements,  but  if  elliptical 
cylinder  or  ellipsoid,  then  are  obtained  complicated  formulas  for  the 
drawing  of  random  variables. 


In  the  present  work  as  the  control  surface  is  selected 
parallelepiped  (5*),  since  for  it  value  P*  will  be  sufficient  high, 
and  the  density  function  of  random  variables  for  faces  are 
simplest. 


In  the  examination  of  flows  within  the  channels  it  is  profitable 
to  develop  flight  of  particle  only  through  the  part  of  control 
surface,  for  example  only  through  one  face  in  plane  of  which  is 
arranged/located  the  entrance.  After  the  drawing  of  the  random 
coordinates  of  particle  on  rectangular  face  St  is  checked  their 
entry/incidence  in  the  section  of  inlet  whose  form  is  arbitrary. 
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3.  On  the  laws  of  the  distribution  of  the  random  parameters  of 
particles. 

We  will  examine  Knudsen's  gas  (Kn-»®) .  If  about  body  (final) 
flows  uniform  equilibrium  (limitless)  flow,  then  the  function  of  the 
distribution  of  the  incoming  particles  is  equal  to  Maxwellian 

1=3 

/o»  =  Soo  (Aoo /w)3fl  exp  [ -  Aob]£b(S| -  Aon  ~  (3.1) 

Here  designation  conventional  (for  example,  see  [1]). 

Let  us  consider  now  formulas  for  the  drawing  of  the  random 
components  of  speed  and  coordinates  of  particles  on  surface  S*. 

Preliminary  analysis  showed  that  the  random  components  of  the 
particle  speed  on  faces  S*  to  more  simply  develop  in  the  Cartesian 
coordinates.  In  the  cylindrical  and  spherical  coordinates  the  random 
variables,  which  characterize  value  and  direction  of  the  particle 
speed,  are  dependent,  which  strongly  complicates  drawing.  Only  when 
V*,  =  0  are  obtained  simple  formulas  for  the  drawing  of  the  velocity 
vector  of  the  particle  isee  formula  (4.10)]. 
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Six  faces  s*  which  satisfy  the  equations 

x,  —  Xu  — 0,  k  —  i;  x(  —  Xji  =  0,  k  =  i-\-  3;  xii<x2i,  (3.2) 

let  us  number  by  index  k=l,  2,  3,  . ..,  6.  Let  us  designate  through 
nm  the  internal  normals  to  faces  Sk. 

Expressing  the  standardized/normalized  particle  fluxes  through 
faces  S»,  we  will  obtain  expressions  for  the  densities  of  random 
components  of  particle  speed: 

components  EM,  parallel  to  standards/normals 

/(£*i)  —  2x_1  (Vr«l)exp  [ — —  V**)*]  ,  .  (3-3) 

i=*k  (Jfe  =  1 ,  2,  3),  i-k- 3(A-4,  5,  6),  8W€(0,  oo); 

components  parallel  to  faces  S„, 

/(Q= («)-«■’  exp  {  -  (U,- W).  (3-4) 

i±k  (ft-  1,  2,  3),  i^fe  —  3  (ft  =  4,  5,  0); 

here 

sin »o cos fc,  —  sin fc);  xW -  exp (- **)•+  V** (1  +  erf x); 

So s  ”  V^oo  A**,  «*/  “  AJi* ; 

afl  and  /3»  -  angles  of  attack  and  slip. 
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Let  us  note  that  function  (3.3)  reaches  maximum  at  the  point 

5o  *  =  0,5  l4  *  +  (0,25  V2k  +  0,5P,  (3.5) 

but  function  (3.4)  -  at  point  Scci- 

Integral  laws  take  the  form: 

distribution  (3.3) 

__  J  /U*<)  &u  =  1  +  Vn* ler* (*«  V'oc *)  —  1]  — 

-exp (-  &  -  V^m-1  (VC*);  <3-6> 

the  distribution  (3.4) 

# 

*fci 

Rkl  M  -  j  fd'ki)  dtki  = -j  [ : I  +  erf  (U<  -  ,)].  (3.7) 

— oo 

For  the  drawing  of  the  random  values  of  components  &*'  it  is 

necessary  to  find  intervals  (5*“,  iu).  the  probability  of  the 

$  * 

entry/incidence  into  which  of  random  variables  ««  is  close  to  one. 

On  the  basis  of  the  preliminary  estimations  for  the  components 
of  speed  parallel  to  standards/normals  nu,  let  us  consider  the 

following  cases: 
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A.  1^08*  >  3,25,  Ini  =V’atk  —  Z,25,  -ki  =  v'mk  -f  3,25; 
6-  0^l/**<3,25,  tf-0.  £^=^*  +  3^5; 

C.  — 4,55<  V»*<0,  &  -  0. 


Page  78. 


For  case  of  A,  by  considering  the  probability 

^[?«<(Vr;*-3,25)l  =  /?w(V';*-3,25)-l- 
-  {V*  I/.'*  [1  +  erf  (3,25)1  + exp  (—3,2s1)}  jr'WLt), 

derivative  of  which 

dRkl  |/~ic exp (—3,25*)  (1  +  erf(3,25))  [  f +  erf  V*k  exp(-V;»t) 
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Accurately  also  for  case  of  B  we  find 

p  10  <  iu  <  ( V‘m  *  +  3,25)j  ( 1  -2,47  •  10~s).  (3. 1 4) 

For  case  of  £  after  numerical  calculations  according  to  formula 
(3.6)  was  obtained  the  empirical  formula 

;«*- 1,523(1/;*  + 4, 55)'«:  (3.15) 

/>(o<£<s;?)^(i-io-*). 

When  l/;*<-4,55  the  quantity  of  molecules,  which  fall  to 
appropriate  face  S*  and  which  constitute  from  all  those  falling  the 
portion,  approximately  equal  to 

Mat 

s,  X  ( *4  [2 X  (  V*  *  )p  =€  i Q-\  (3. 16) 

*=1 

can  be  disregarded/neglected.  Here  Sk—  area  of  face  k. 

Is  analogous,  for  parallel  to  faces  5»,  let  us  find  that 

p  [($« ,  -  3,25)  <  (u  <  (SOT  ,  +  3,25) j  >  ( 1  - 10"»).  (3. 1 7) 

Since  (3.6)  and  (3.7)  are  not  solved  relative  to  &«»  then  the 
drawing  of  random  values  ««  is  produced  according  to  Baird's  method 
[12].  In  this  method  in  comparison  with  Neumann's  method  [25]  the 
condition  for  the  selection. of  random  numbers  does  not  require  gaming 
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the  further  irregularly  distributed  random  number,  designated  further 
by  symbol  R,  since  it  is  formed/shaped  on  the  course  of  calculations. 

The  random  coordinates  of  the  start  of  molecules  have  density 
function: 

/(*!„*)  —  8(*j- JCw);  f(Xi+k)  =  {xu  —  *—1,2,3;  (3.18) 

f(Xls*-3)  —  *(x,  —  Xu);  f(xn.k-3)  =  (xu  —  x'u)-\  *-4,5,6;  (3.19) 


here  5(x)  -  Dirac's  delta-function. 

Since  the  drawing  of  trajectories  is  produced  relative  to 
coordinates  xxlt  then  to  components  &  of  particle  speed  after 
drawing  according  to  formula  (3.3)  for  cases  of  k=4,  5,  6  it  is 
necessary  to  confer  minus  sign. 

Page  79. 

4.  Drawing  of  random  particle  trajectories  under  different  laws  of 
reflection  of  particles  from  the  wall. 

Let  us  designate  xlt  and  xrl—  the  coordinate  of  initial  and  end 
point  of  one  of  the  elements/cells  of  break  -  particle  trajectory, 
which  tests/experiences  collisions  with  the  concave  surface,  with  the 
start  with  S*  we  have 
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•*«'<*=*>"■  Xu,  xll(ti.k)=zx]i  +  /?,  (x'u  —  x\i)\  *=1,2,3;  1  ^ 

X\  i  =  x j  i,  X\  t  *_3>  =  JCi  i  -f-  (jc2  i  —  Xi  /);  k  —  4,  5,  6.  j 

From  point  Xtl  the  molecule  moves  over  the  straight  line 

xi  —  xu  +  l-\ it;  «_u  =  £_ii(  £  jit  i)~hi.  (4.0) 

r=:l 

Index  -  "in  the  speeds  designates  the  emitted,  "  +  "  arriving 
flying  particle,  and  superscript"  0  "designates  the  unit  vector  of 
vector  of  velocity.  Solving  system  from  equations  (1.4)  and  (4.2),  we 
determine  values  of  ta  and  ta: 

fi . *  -  ( -  * ±  (»* - 4  “>' "'l  (2  a)—,  ^ 

where  &  =  £— -i  i  £— 1 1!  b  **  2  (alt  Xj  ^  Li  *  “f-  flj  *  Li  i);  r  =  ^  xt  4  4*  2  <z(  4  jc,  (-f-  au, 

for  which  straight  line  (4.2)  can  intersect  surface  (1.4). 

Due  to  the  final  accuracy  of  arithmetic  operations  by  ETsVM 
[  -  digital  computer]  the  coordinates  of  the  collision  of 

molecules  with  the  surface  can  be  computed  with  the  "short  round"  or 
the  "flight/passage".  In  order  to  eliminate  the  effects  of  these 
cases,  we  will  accept,  that  with 

(b*  —  4  or)  <  0;  *,<«„  f*  <  «i  (44) 

the  molecule  does  not  encounter  surface,  but  with 

/,  >•„  —  oo<f,<  00 ;  —  oc<f,<oo  (4.5) 
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is  possible  the  collision  of  particle  with  the  surface  with  tx  or  t, 
or  when  *i.2>  Calculations  show  that  it  is  possible  to  take  e^lO"5 

Thus,  are  found  one  or  two  points  of  intersection  of  particle 
with  the  surface 

*m  —  +  tli  itm,  m  -  1 ,  2.  (4.6) 

After  this  is  produced  checking  the  entry/incidence  of  points 
xlm  to  the  assigned  piece  of  surface,  determined  by  D  region,  for 
which  it  is  necessary  to  calculate  coordinates  in  that  system,  i 

which  is  assigned  D  region.  Values  with  which  the 

intersection  of  surface  did  not  take  place  or  the  corresponding 
points  do  not  fall  into  D  region,  are  substituted  by  number  LX>L, ;  L 
is  equal  to  maximum  size  Sk. 

The  cycle  of  this  checking  is  completed  for  all  surfaces  (1.4) 
of  body.  From  the  obtained  set  t  is  determined  in-  If  then 

particle  does  not  intersect  body,  otherwise  intersects  it  at  the 
point 

xrl  —  I  +  $lu  tmlm .  (4.7) 

The  reflection  of  particle  from  point  xH  is  realized  depending 
on  the  assigned  boundary  conditions  on  the  surface. 
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With  the  mirror  reflection  of  molecule  the  speed  of  reflections 
of  molecule  is  equal  to 

5—  ri  =  £+W  2  (S+rt  ®|)  (f -8) 

Here  ni—  -  internal  normal  to  the  surface. 

Page  80. 

In  the  case  of  diffuse  reflection  from  the  surface  is  introduced 
local  system  of  coordinates  with  beginning  W(xr'd  and  axes  w» 

such,  that  cosines  cnl  of  the  angles  between  and  are  equal  to 

ci  !“{£+'■< — «i)  rtflP  —  ($+rtni¥]-117’' 

(£„*,)<  0.  (4.9) 

If  (C.h«j)>0*  then  vectors  ctl  and  ctl  vary  directions  to  the 
opposite  ones;  if  (Cw«i)*»4;l  then,  accepting  cx,=0  and  taking  into 
account  that  knl~l,  (cuntisss0,  we  will  obtain  cll»-n,c1,/nw  cxl=- 
(l+(n,/n1 ) * ] *  y*; 

The  speed  of  the  diffuse  reflected  particle  is  developed 
according  to  the  formulas  (see  also  [11]) 

Z-rt  —  ;  £-h  —  cln =  (cos  psln  9, sin? sin 8, cos6}; 

?  =  2ir/?,;  sin  0  =  I  Ui  =  (-Ar'ln(/?4 /?,))'*; 

Kx  -  «.(Aw)-‘  +  ~a»)  . 


(4.10) 
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Here  a,—  accommodation  coefficient  of  energy  with  the 
mirror-diffused  reflection. 

With  the  mirror-diffused  reflection  each  molecule  with 
probability  p0  is  reflected  diffuse,  i.e.,  if  p*>R,,  then  reflection 
is  diffuse,  otherwise  it  is  mirror. 

Let  us  consider  law  of  reflection  of  molecules,  given  by 
Maxwell's  function 

A-MW’expf— VJ-  (4.11) 

L  »= i 

with  five  macroparameters  nnhr,Vm.  Designating  the  average/mean 
values  of  the  components  of  the  speed  of  the  molecules  reflected  in 
system  Wwn,  constructed  in  accordance  with  formulas  (4.9),  through 
—  {«n»  0,  T,,},  the  average/mean  value  of  energy  of  the  molecules  (per 
unit  of  mass)  reflected  -  through  >-  and  the  relative  average  speed 
-  through  Srn  —  {S,„  0,  5,,)= Vmh\n,  we  will  obtain 

*7,=  ^ sine,;  T„  =  vy;sr;  F*/2-^[2+  S2r  +  SriC?'-Sti)\;  ■ 
y  =  S„+0.5V^(l  +  erf5„)z-l(5,3);  cos6,; - (*,*6).  (4.12) 

Hence  we  obtain  the  transcendental  equation 
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'b  =  ?/,(2  +  5,,«')-1  -1?3(52  — I?,)-1*  (4.13) 

It  is  evident  that  this  model  does  not  limit  the  physical 

picture  of  phenomenon,  since  the  left  side  where  does  not 

contradict  right'.  3^ "yg, 

/*=  1 

Average/mean  values  ?r"  were  determined  with  14,^10*  cm/s  in 
work  [10]  in  the  form 

a,  =  5rl(tf-S?);  *»=sr‘ (C.J— 51»);  «t***r'(5u  -£i)-  (4.14) 

Components  *i«  and  *rH  are  expressed  in  system  Wwn- 

From  (4.13)  and  (4.14)  we  obtain  the  condition  of  the 
solvability 

(sin^  —  *,)*  +  (a„  —  cos«,)J<(l  —  a,),  cos  8t  =  (n,  C  u  ),  (4.15) 

i.e.,  the  terminuses  of  vector  (<*,,  *„)  must  be  placed  in  upper  half  of 
circle  with  the  center  at  point  (sin  6lf  cos  0X)  and  with  radius 
(1  — 
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The  sequence  of  calculations  is  the  following. 


DOC  =  83067303 


PAGE  Mi 7 


1.  Are  calculated  ®i,  <**,  E?,  fr9l  E,,  and  \p. 

2.  Is  solved  equation  (4.13).  With  0.04^^.0.9612  was  applied  the 
method  of  linear  interpolation.  With  i£>0.9612  it  is  obtained  with 
error  0.1% 

Sr ,  =  ( [1 ,5  +  (16 1  -  1 3,75 )«»]  (2(1-  «(.)]-» }  W,  (4. 16) 

while  with  ^<0.04 

-  -[(1+0  -  18 1)«*) (2t)“,F.  (4.17) 

These  relationships/ratios  are  found  from  asymptotic  expansions 
of  function  'H5r»)  when  S,t-+  oo  and  5,, -*•  —  oo . 

3.  It  is  calculated  <p',  then 

Sn^f'Ur,;  ^  =  Jr(Ar)-'«;  sj  —  5?  I  +  S*3-  (4.18) 

The  drawing  of  the  components  of  speed  Er„  of  particle  is 
produced  employing  the  procedure,  presented  in  Section  3.  Then  is 
computed  &-*  “*!/«&*»•  it  is  obvious,  in  the  general  case  with  the 
multiple  reflections  of  particle  from  the  concave  surface 
coefficients  a,  and  it  is  not  possible  to  consider  constants. 

Assuming  that  the  particle,  which  falls  at  a  high  speed  to  the 
surface,  reduces  its  energy  due  to  the  multiple  reflections  and  it 
reaches  equilibrium  with  the  surface,  with  which 
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we  will  describe  the  process  of  achieving  the  equilibrium,  by  the 
dependences,  which  give  a  change  in  the  values,  which  characterize 
energy  of  the  particle: 

52  ir2 = 2 h?  sr2  +  ?,[(i-  *,)  -  2A-1  it2}-,  lh  it2 = sr2 + 

+  ?n  IK  -  cos  8,)2  -  JL  A-'  IT2}-,  Vi  5r2  -  9,  (sin  8,  -  a,)2.  (4.20) 

Coefficients  <p«.  and  vary  from  one  (with  the  high 
energies)  to  zero  (with  low  energies).  Were  examined  the  following 
cases: 

1)  coefficients  <P„  9*  and  9-.  depend  on  the  number  of  multiple 

reflections  on  the  exponential  dependence 

—  <p„  =  —  exp  (ao,)r  a<0,  ax—p  —  p0,  (421) 

on  the  linear  dependence 

<P,  =  ?n  =  ?*  —  (pw  -  p)  (Pot  -  Po)~  Poo  -  Po  >  1 ;  (4-22) 

2)  coefficients  ?„  and  9*  are  assigned  by  the  arbitrary 
monotonic  functions  of  energy  of  incident  particles. 

In  formulas  (4.21),  (4.22)  p  -  reference  number  of  the  collision 
of  particle;  p,  -  number  of  reflections  (inclusively),  in  which 
during  calculation  and  coefficients  a,  and  «r  are  assumed 

to  be  constants;  a<0  -  arbitrary  coefficient.  Reflection  with  the 


DOC  =  83067303 


PAGE 


number  of  collision  p00  (inclusively)  in  (4.22)  becomes  diffuse  with 
the  complete  accommodation.  Let  us  note  that  assigned  dependences 
a„  <*„,  a,,  f>„  ?„  and  must  satisfy  the  condition  of  solvability  (4.15). 

5.  Calculation  of  the  total  particle  fluxes,  impulse/momentum/pulse 
and  energy.  Calculation  of  dispersions. 

Let  us  designate  the  flows  from  the  first  collision  by  index 
from  first  reflection  and  from  the  second  collision  and  the 
reflection  and  so  forth  to  the  escape  of  particle  into  infinity  from 
point  xrl —  by  index  r.  Particle  fluxes  relate  to  n^Va, 
impulses/momenta/pulses  -  to  mn*Vi/2,  energies  -  to  4mn®  Vl/2.  During 
the  calculation  of  aerodynamic  flux  coefficients  of  forces* relate  to 
moments/torques  +o  -nut* V*SudJ2.  when  — 0  instead  of 
V*  is  used  A51;* 

Page  82. 

Mathematical  expectations  of  the  aerodynamic  force  coefficients 
and  moments/torques  of  the  compound: 
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c±i 


V-=M„  . 

=  ±A  2  *±u;  C«-A  2  (5+a-5-«); 

iSi  i1-1 

j—  2  (-^1  i  X  »±u)'» 


*"  *=1 


4  ,l=,Y«  /  .  .  x.  a  _  2BPn  . 

2  K^iXS+a)  — (^iXs-rt)!;  A—  \z(nNasu 


mr 


A  e,  k=A 

B  =  -p §-  2 

/ w5ao  *Tl 


(5-1) 


Here  u  -  number  of  the  particle,  played  with  £*  and  that  fallen 
on  body.  Total  coefficients  of  complicated  body: 

ct  —  c+i  -4-c—  i  +  ml  —  m+i  -I-  m~i  +  m,i-  (5-2) 


Examining  the  sequence  of  the  random  values  of  the  aerodynamic 
coefficients,  obtained  according  to  the  results  of  the  drawing  of 
each  u-th  particle  trajectory,  we  will  obtain  for  calculating  the 
dispersions  the  following  recursion  relations: 


M[X]=a0M}>;  D[X\  —  Oo* D*  Af„  —  [G*  -  1) M*_i  +  ZJ/k 
D*  =  {[D*-.  +  G*  -  l)+(Z,  -  Mv)2}l\x\  a0  _  . 


(5.3) 


Here  M  and  D  designate  mathematical  expectation  and  dispersion; 
Z^,  *  (1+n  ^ — a ) ,  —  3^  npH  s  Cii 

Z(i.  —  (Xj  4  X  5+  If)  (^rt  X  ri)>  =>  SHdH  /ipH  X  ^  mt. 


Key:  (1).  with. 
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The  central  axis  of  the  system  of  the  forces,  applied  to  the 
body,  is  parallel  to  vector  ct  and  it  passes  through  the  point 

xi »  =  —  dH  [nl  X  ct\  [|  ml  jJ  |  ct  j*  -  (ct  mt)'\  [|  ct  | 2 1  ml  X  ct  (5.4) 

moreover  moment  with  respect  to  this  point 

mio  —  c<(C/'Ti/)|ri|-2,  /HjoJc,.  (5.5) 


Flux  coefficients  of  energy  to  the  body: 


*mNn  , 


£=£+ 4- £_  +  £,;  E±  =s+  At  2  4>;  Er=Ai  2 

|U>1  |ut 

A  Bpa 
1  ' 


(5.6) 


Flux  coefficient  of  the  energy  through  the  control  surface 

E»  —  Ax  2  1 5*,  |*. 

1*=*! 

When  value  E+  is  approximately  equal  to  the  area  of  the 

"shadow"  of  complicated  body  for  the  plane,  perpendicular 


Is  calculated  also  relative  value  of  further  collisions 
(A/«,  +  N>)A« ,  where  and  N, —  quantity  of  all  particles,  which 

fall  on  body  from  infinity  and  after  the  first  reflection. 


With  the  diffuse  reflection  of  particles  from  the  concave 
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surface  with  the  complete  accommodation  easily  is 

establ ished/installed  the  dependence  of  coefficients  on  temperature 
of  the  surface: 

cri~SZ',  E-t~SZ\  E,~SZ\  S.-Sa,(7W7’Jw*.  (5.7) 

Here  and  throughout  symbol  designates  proportionality. 

Page  83. 

6.  Calculation  of  the  aerodynamic  interaction  coefficients  of 
surfaces  and  bodies. 

We  will  examine  the  aerodynamic  force  coefficients  and 
moments/torques,  which  function  on  the  surface  of  concave  body  as  a 
result  of  multiple  collisions  of  particles  with  the  surface.  Indices 
and  r  will  designate  flows  from  the  first  collision,  from 
the  first  reflection  and  from  the  multiple  collisions  of  particles 
with  the  surface.  Aerodynamic  force  coefficients  relate  to 
fnthoVlaSJi  moments/torques  -  to  mnvty^tSytdytj2.  When  V*  — 0  instead  of 
Voe  is  used  . 

Mathematical  expectations  of  the  aerodynamic  force  coefficients 
and  the  moments/torques,  which  function  on  v-th  surface  of 
complicated  body  in  body  axes  Xxt: 
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•*-*0 

C±*l  —  A  2  S±l»i  I  Cnl  A  ^  (5+  2,i  £ — 2W  ~f~  ...  £— r«<  ): 

|i-l  |wl 


A  ^ 

fft±tl  =:  i  ~j  2  (^1  i  X  ^ilxi  )» 

““  |m] 


fflni 


~T~  2  It"**  I  X  2+J»/  )  —  (-*2  i  X  &— 2ii  )  “{“  •  •  ■  —  (xri  X  £— rti)]; 

““  l*=l 


4 


2S/>„  . 


Q.5  y 

/it  5.  A 


x(v;*)5,. 


(6.1) 


In  these  formulas  are  summarized  the  impulses/momenta/pulses 
during  the  incidence/drop  in  the  particle  on  v-th  surface  with  the 
drawing  of  all  ^==\-*-Na  particle  trajectories. 


Aerodynamic  coefficients  of  the  surface: 

c,|S=c+,(-t'C-^  +  cni;  m,<  —  m+w+  ffi_w  +  mr,i,  (6.2) 


If  complex  body  consists  of  several  separate  bodies,  then, 
summarizing  coefficients  (6.1)  for  the  surfaces  of  which  consists 
separate  body,  it  is  possible  to  compute  the  aerodynamic  interaction 
coefficients  of  the  bodies,  which  fly  on  certain  distance  from  each 
other. 


7.  Calculation  of  the  local  particle  fluxes,  impulse/momentum/pulse 
and  energy  through  the  surface  of  complex  body. 
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The  coefficients  of  local  particle  fluxes  («+  and  «,)  we  will 
carry  to  «« V9\ normal  (P*  and  Pr )  and  tangential  and  xr) 
impulses/momenta/pulses  -  to  mn^  l/®/ 2  and  energies  (e±  «rVk)  'nntovla/2.- 


For  calculating  the  local  flows  the  surface  is  divided  on  the 
pads  AS:  first  is  divided  the  rectangular  region  D  of  the  assignment 
to  surface  On  the  large  number  m*  of  equal  rectangles  with  the  help 
of  2(m-l)  straight  lines,  carried  out  through  each  side  of  rectangle 
D  in  parallel  to  its  sides.  Let  the  centers  of  all  areas/sites  AS 
have  integral  coordinates  on  the  external  variable/alternating  region 
D/3=l-m  and  in  terms  of  the  internal  variable/alternating  7=l-m  and 
the  general/common/total  numbering: 

•  -(P—  i)«  +  t»  •e(i+«‘).  (7.i) 


Page  84. 

Each  of  m1  rectangles  D  region  determines  on  the  surface 
area/site  of  the  arbitrary  form  whose  area  is  computed  from  the 

formula 

*»+V  r  t.+^t  i  as'  - 

AS*“  J  [  j  n^h)dzAdz%-  A7=  Jl_Il;(7>2) 

here  (/3lf  /32),  {ylt  y2)  -  limits  D  region.  Functions  h)  and  the 
variable/alternating  during  the  assignment  regions  D  take  the 
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form: 


in  Cartesian  coordinates:  ( y} ) 

* 

f(zf,  zT)=  — /fc-1;  z9^yx;  z,myt;  (7.3) 

in  the  polar  coordinates  (<p,  r,  y,): 

f(z9,  z1)  =  —  rn~u,  z?* i?;  Zfmr,  (7  4) 

in  the  cylindrical  coordinates  (<p,  r,  y ,) 

f(zt,  —  —  r  («,  cos  <p  4-  /»,  sin  <?)-*;  z9  —  ?;  r7  as  jr,;  (7.5) 

in  the  spherical  coordinates  (<p,  0,  r) 

f(z9,  2,)=  —  sin  M(«,  cos  «p  sin  6  4-  /*,  sin  ?  sin  6  +  nt  cos  8)—*; 

2T  =  9i  (7.6) 

/a  ,  \-i/a 

;  /y  =  -  <?/^(^  yj/fy  =  -  (^ -f b)  t).  (7.7) 


The.  components  of  the  unit  vector  of  internal  standard/normal 
rij,  value  ya  in  the  Cartesian  and  polar  coordinate  systems,  value 
in  the  cylindrical  and  spherical  coordinates  are  computed  with  the 
help  of  the  equation  of  the  surface 

F(yj,  yn)  =*  *y* j ')  yn  +  bi *yi  +  b** =°-  O'.  «=  i.  2»  3)»  (7-8> 

moreover  in  the  solutions  of  corresponding  square  equations  before 
the  radicals  is  taken  sign 

By  the  values  of  the  coordinates  of  the  impact  point  in  the 
particle  on  the  surface  are  determined  coordinates  (z*.  zt)  in  D 
region,  through  which  is  located  the  number 


DOC  =  83067303 


PAGE  (p 


«,  =  «(£— 1)  +  *f.  (7.9) 

where 

P  -  entier  [( z%  -  ft)  AjT1  +  11,  t  =  entier  [(2T  -  Tl)  b~l  +  1] , 

using  the  address  of  area/site  &sm.  Here  entier  (x)  indicates  the 
near  whole,  which  does  not  exceed  x. 

For  each  area/site  in  terms  of  the  value  of  address  u  occurs  the 
accumulation  of  the  corresponding  sums  from  which  are  computed  the 
flux  coefficients: 

the  normal  impulses/momenta/pulses 

psAfg  l*  •  /Vjj 

p±m1=iAi  2  |5±i«irt||;  Prm*mAi  2  (|£+sw«j|  +  |5-»*«j|  4-  •••  +  «il):  (7.10) 

|l»l  ll=sl 

the  tangential  impulses/momenta/pulses 


*±-=±At  2  d±i-e.i);  *. - A,  2  (7.11) 

t»«i  1,-1 


the  energies 


^.-^2  *  Ai  2  —  E— 2» -|-  ...  —  (7.12) 

t»S*I  (Ml 


the  particles 
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(7.13> 

(7.14) 

If  (vaoinl)  =  ±  l,  then  "  {s,n3o,  cosa,,  0}. 

Page  85. 


/!+.  —  Y  (*!«.;  rirm  =  Ai  y  0**.+  1*3-4*—  + 

Bpn 


H=1 

2Bp„ 


a  - m  .  A  —  _ _ 

i/  \r  a  o  *  s  i /2 


;  ^44  — 


Bpa  . 

ALAS.-’ 


V’aaJV.A5.  ’  3 

(x^ooi  «i)  «|]  [1  —  (t»aO/ ^QOi  =  V'oo  (  I 


-1 


In  formulas  (7.13)  fM-=I  during  the  incidence/drop  the 
particles  on  &Sm  from  infinity,  otherwise  im„  =  o;  w«.  (*a»»  —  .  Pr.  are 
equal  to  one  during  the  incidence/drop  in  the  particle  on  AS.  after 
the  first  reflection,  otherwise  they  are  equal  to  zero. 


In  the  case  of  the  symmetry  of  flow  can  be  of  Interest  medium 
stresses/voltages  along  the  bands,  parallel  to  the  sides  D  region. 
Designating  by  the  index  8  averaged  stresses  over  the  bands  in 
parallel  to  side  y  of  D  region,  and  by  index  y  the  averaged 
stresses/voltages  on  the  bands,  parallel  to  the  side  B  D  region,  we 
will  obtain 


1  T-"  1  »-m 

?+*  D,  S/+"  P+f—  —  2  P+-’  ••• 


fcT*l  __  T.  ?=1 

and  medium  stresses/voltages  on  over  the  entire  surface 


(7-15) 


«  •««* 

P+  =  : i  2  Pi¬ 
rn  m»l 


(7.16) 


( 


X 
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Varying  in  formulas  (7.15)  and  (7.16)  p  on  r,  e  and  n,  and  index 
"+"  on  and  r,  we  will  obtain  formulas  for  all  averaged  flows. 


With  the  diffuse  reflection  of  particles  from  the  concave 
surface  with  the  complete  accommodation  easily  is 

established/installed  the  dependence  of  coefficients  on  temperature 
**•  of  the  surface: 

p_m~sz'-,  pn»~szK,  'cr-~s»1;  «—  er»~<SE*0-i7) 

8.  Calculation  of  the  local  parameters  of  the  flow  of  the  strongly 
rarefied  gas  in  the  vicinity  of  compound. 


Let  us  consider  the  procedure  of  calculation  of  the  penetration 
probabilities  of  the  particles  through  the  specific  surfaces,  and 
also  densities,  the  average  speeds  and  the  temperatures  of  gas, 
stress  tensors  and  fluxes  of  momentum,  vectors  of  complete  and 
thermal  energy,  mass  flow  rates  and  static  pressures  at  arbitrary 
points  in  the  vicinity  of  compound. 


Let  us  introduce  formally  into  the  examination  the  flat  surface 
which  will  serve  for  fixing  of  the  parameters  of  gas  and  will  not 
affect  particle  trajectory.  Index  X  will  designate  the  number  of  this 
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"formal"  surface,  Sx  —  its  area. 


We  will  examine  the  flows  of  sign/criterion  <p,  transferred  by 
the  particles  through  the  "formal"  surface  from  the  outer  side  [index 
"  +  (E+y tij) >0J  and  from  inside  [index  The  parameters 

without  the  index  will  relate  to  entire  particle  flux.  Let  f  - 
function  of  particle  distribution  in  the  vicinity  of  this  surface. 

Let  us  write  formula  for  the  flow  of  sign/criterion  <p  through  surface 


Si: 


SI 


1  [/  /  G+'nj) dS*  d*J  +  j  f  nj)  dS\  ft,  ]  = 

Sl  (E+y  «y)X>  (E_y  Mj)  <0 

n^BVv,  1 


(8.1) 


Here  is  conducted  averaging  over  surface  which  must  be  sufficient 
small . 

The  penetration  probabilities  of  the  particles  through  the 
formal  surface  X  from  the  external  and  from  inside  are  equal  to 
lF±x-Xx  A/— ‘ .  Here  N±x—-  number  of  molecules  from  Np,  that  fell  on 
formal  surface  from  the  external  or  from  inside. 

Page  86. 


The  densities  (concentration)  of  particles,  averaged  over 
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surface  Sx,  through  which  they  pass,  and  those  referred  to  the 
particle  density  on  infinity  nmt  are  expressed  from  formula  (8.1) 
when  ?  — fey/iy)-': 


nix  -  «*x/i*o.  -±  Bt  SC1  *))-'■’  nx-nxln*,  -  + *-*. 


BX-BV*N?'. 


(8.2) 


The  average/mean  mass  flow  rates  of  molecules,  in  reference  to 
are  obtained  from  (8.1)  when  ?“s&/(S±/*y)'“1: 


\-l. 


^,1)^=  ±  S2(Svnl,rl  S 

.  !»■=  t 


(83) 


Taking  into  account  that  the  temperature  of  gas  is  expressed  as 
mean-kinetic  energy  of  the  particles 

_ i  r% 

(8.4> 


where  n,  m  -  density  and  the  mass  of  particles,  we  will  obtain  from 
(8.1)  when  =  fey'1/)-1 


Tlx-l^-ifT^  (n±x  Sx)~ 1  -  V4\  ; 

(T+x.  +  r_x.)  («;  Sx)-1  -  y  tf2;  r±x.=±  0,  "2 ' V  fey  «>)“ 

B,  =*fl2Aw  ;  V^x  =  |^'±/(x,|;  K-|VyW|.  ’  (8.5) 
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Determining  the  tensor  of  the  flux  of  momentum  by  the  expression 

=  Mtj  (mnoo  V 1 :  Mtj  —  mftiijfdi,,  (8.6) 

we  will  obtain  from  (8.1)  when  ?  —  5*Sy(S±/«/)“l 

M‘±iJ  (X,  =»+  B^sr1  "v'  Ij  (E±  J  n,rl ;  Af;/(X)-Af+*/(X)  +  ^-«/<X).  (8-7) 

H-=I 


Determining  stress  tensor  by  the  expression 


Pv  =  P«/  (««»  Ai ')- »;  Pij  =  m  f  fa  —  l/s)  (Sy  _s  l/y)/c«y  = 

—  m  S  M  &//<*£/  -  m/I  l/y,  (8.8) 

we  will  obtain  from  (8.1)  when  ®  —  SiSy^yrty)-*1 

P-*J  (X) =  (X)  (X)  l^y  (X) ,  psj  (t)  =  /Mjy  —  /tv  VTs  (i)  Vj  (\f  .  (8.9) 

Determining  the  vector  of  the  flow  of  total  energy  by  the 
expression 

E'/=E)^~n^ft^nyi ;  (8.10) 


we  will  obtain  from  (8.1)  when  <p  —  SMyfS+y/iy)-1 


fi/ft>-±54Sx-l^5 *5±y(?±/«y)-*;  E'/(^E'+/M  +  E'.m\  B<-B>hg.  (8.11) 

- 
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Determining  the  vector  of  thermal  energy  by  the  expression 

q, -  jtf-  ^*(5/-  Vj)fdt„  (8.12) 
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we  will  obtain  from  (3.1) 

P±t/M~  T±\»S\  ;  |  (8.13) 

q)  (x>  —  E  i  w  —  2  Vi  (n)  p [j  (i)  —  Vj  (\)  (T+x*  +  T -x«)  ST1.  I 

Static  pressure  P\ 

p\  =  pxpZ1  —  rt\T\ , 

but  average/mean  particle  flux  through  S* 

ul  =  u\  (mrtaa  hZwY~ 1  =  n>-  ; 

here  v'nX  projection  K  on  the  normal  to  Sx. 

Let  us  note  that  the  started  from  the  formal  area  of  the 
particle  can  for  a  second  time  clash  with  the  same  surface  due  to  art 
error  in  the  calculations.  So  that  this  could  not  occur,  the 
coordinates  of  particle  with  the  start  must  be  displaced  by  the  low 
value  e ,*10 " 4  in  the  direction  of  particle  motion. 

9.  Calculation  of  aerodynamic  coefficients  in  the  flow  of  Newton  and 
in  the  flow  of  light/world. 

Employing  the  procedure,  comprised  for  calculating  the  free 
molecular  flows,  it  is  possible  to  calculate  the  aerodynamic 
coefficients  of  the  bodies,  which  fly  in  the  flows  of  Newton  and 


(8.14) 

(8.15)  ‘ 
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light/world. 


The  total  aerodynamic  coefficients  of  bodies  and  the 

coefficients  of  the  local  flows  of  normal  impulse/momentum/pulse 
^•*in  Newton's  flow  are  equal  to 


Ct  h  —  (c+i  +  c-i)l 2;  mIH  =  (m+i  +  «_,)/ 2;  (9. 1 ) 

/>•«=*  (/>+•  +  />— )/ 2.  (9.2) 

Entering  the  right  sides  of  these  formulas  coefficients  c±i<  m±‘' 
and  P±»  [see  formulas  (5.1)  (7.10)  are  calculated  from  the  procedure 
of  calculation  of  free  molecular  flows  under  the  following 
conditions:  the  falling/incident  from  infinity  particle 

encounters  surface  only  one  time  and  is  reflected  mirror. 


Sometimes  appears  the  need  for  calculating  the  aerodynamic 
characteristics  of  bodies  in  the  luminous  flux,  which  proceeds  from 


any  point  source.  The  luminous  flux  consists  of- the  photons,  which 


have  energy  *♦  —  mass  m4,  =  —  ,  impulse/momentum/pulse  p*=-~  and 

cl  e 


numerical  density  n  _ here  h  -  Planck's  constant;  c0,  v  -  speed 


k*c% 


and  the  frequency  of  light/world;  light  energy,  which  falls  per 

unit  time  Der  unit  of  area. 


With  the  drawing  of  particle-photons  necessary  to  assign 
Soo^l;  the  speed  of  the  falling/incident  and  reflected  particles 


DOC  =  83067303 


PAGE 


are  equal  to  the  designed  aerodynamic  characteristics  will  be 

related  to  velocity  head 


0  _  __ 
2  —  2ct  ’ 


(9.3> 


Let  k2,  k2  and  k2  be  equal  to  the  portions  of  that  reflected, 
absorbed  and  passing  through  the  surface  of  the  luminous  flux, 
moreover  ki+k,+k,=l.  With  the  drawing  of  the  reflection  of 
particle-photon  from  the  surface  of  the  probability  of  its 
reflection,  absorption  and  passage  through  the  surface  are  equal  to 
respectively  k2,  k2  and  k2. 

We  will  examine  the  surface,  which  does  not  pass  light/world 
(k2=0)  and  isolated/insulated,  so  that  the  absorbed  photons  will  be 
again  scattered  according  to  "cosine  law"  (diffuse  reflection)  with 
the  function  of  the  distribution 

/«-  (9-4) 

Here  5(c-c0)  -  Dirac's  delta-function.  In  this  case  with 
probability  p0=k2(l-kl)  particle-photon  will  be  reflected  diffusively 
and  with  probability  k  -  it  is  mirror. 


Let  us  note  that  the  worked  out  methodology  makes  it  possible  to 


the  diffusion  of  light  in  the  arbitrary 
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Page  88. 

Chapter  II. 

DESCRIPTION  OF  UNIVERSAL  PROGRAM. 

The  program,  given  in  appendix  1,  is  comprised  in  accordance 
with  the  algorithm,  described  in  Chapter  I.  Turnings  to  IS-2,  to  the 
standard  programs  (SP)  for  calculating  the  functions  ln(x),  exp(x), 
sin(<),  arcsin(x}  arctg  00,  erfCA  to  SP  of  the  multiplication  of 
matrix/die  by  the  vector,  the  multiplication  of  matrices/dies  and  to 
SP  of  the  group  translation/conversion  into  the  binary  system  from 
the  decimal  (10-»2)  of  the  array  of  numbers  are  written  in  accordance 
with  [26].  Remaining  turnings  to  IS-2  are  given  in  Table  1  and  are 
isolated  in  the  program  with  feature  to  the  left. 
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*_  H--HrC°n“-  BTopaa  crpoica*  o<SpameHaa  k  Cfl 


ripHMeqaHHH 


0041  lleiaTb  Mac-  *,0x,  52 

CHaa  qacea 
or  a  Ao  b  c 
nepeaoAOM 
H3  2-»10 


CoAepucHMoe  aqeeic  a—b  ae  nop* 

TMTCft  _ 

To  . 


0150  ^HaeftHaa  ooo  °0  </»>  0150  <jcq>  n  —  mhcjo  y3aoa  RHiep- 

HHTepnoaa-  0_  7«  noaauaa  6ea  OAaoro; 

aaa  <t>ymc-  1  *•  5  </(*o)>  </(-«)>  ^  /(jc0)  —  apryMeat  a  ^jukuhs 

uafl  oahoB  a  nepaoM  y3Jie  aaiep- 

aepeMeHHofi  noaaaRH;  (£> 

s-  x  —  reaymee  3u<ieHHe  ap- 

_ (.  1)  _ ryaeHTa  (£) _ 


0122  BfeiiHCJxeHHe  000  00  <a> 

AaoflHoro  (jqo  QO  <+  > 
aarerpaaa 

MCTOAOMfei  000  00  <+«j> 


CHMncoHa^J;  Qoo  ^  <e(x)><f(x,y)>  <x> 


/=  f  F(x)dx,  F( x)~  f  / (x, y)dy\ 

Q2)a  c  (jt> 

apryMeHTU  x  w  y  Oepytca  aa 
<*>  a  <jc>+1;  a  weflaax 
7-t*7+m — 1  cieT  c( x)  a  d(x), 
npaqeM  e  (je)-><c(x)>,  a  d(x)~* 
<c(-f )>+•*;  »  aqefticax 
+  r—  1  cqet  f  (x,  y);  7  +  m  a  3  -f- 
+  r—  nycrue  aqeflKa;  norpemao- 
cth  ±  «,  a  +  *,  auqHcaeHHii  /  a 
F(x)  a0coanTHMe,corAa  3aax  (— ), 
a  OTHOCHTeJibHwe,  Koraa  3«ar  (*~); 
3aaqeaae  l  a  aqeflxe  0001 


Key:  (1).  Number.  (2).  Designation.  (3).  Second  row  1  of  turning  to 
SP. 

FOOTNOTE  1 .  The  first  row,  which  is  located  in  nucleus  x-1,  takes  the 
form  016  x  7501  7610J  <*>  -  cell  in  which  is  located  number  x. 
ENDFOOTNOTE. 

(4).  Notes.  (5).  Press/printing  array  of  numbers  from  a  to  b  with 
translation/conversion  from  2-*10.  (6).  Contents  of  nuclei  a-b  does 
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not  spoil.  (7).  Linear  interpolation  of  functions  of  one 
variable/alternating.  (8).  number  of  interpolation  points  without 
one;  argument  and  function  in  first  interpolation  point.  (9). 
instantaneous  value  of  argument.  (10).  Calculation  of  double  integral 
by  method.  (11).  Simpson.  (12).  argument  x  and  y  they  are  taken  from 
<x>  and  <x>+l;  in  nuclei  7-7+m-l  calculation  c(x)  and  d(x) ,  moreover 
c(x)— ->  <c(x)>,  and  GO — >  <c(x)>+l;  in  nuclei  a-cr+r-1 

calculation  f(x,  y);  7+m  and  0+ r  -  empty  nuclei;  error  ±et  andie, 
calculation  I  and  F(x)  absolute,  when  sign  (-),  and  relative,  when 
sign  (+);  value  I  in  nucleus  0001. 

Page  89. 

I.  Block  1(0)  for  calculating  of  total  aerodynamic  coefficients  and 
local  particle  fluxes,  impulse/momentum/pulse  and  energy. 

This  block,  which  occupies  nuclei  0001-6763,  is  basic  the  work 
of  all  remaining  blocks  II-IX  it  occurs  with  block  1(0). 

During  calculations  with  block  the  I  (0)  reflection  of  molecules 
is  mirror-diffused,  moreover  the  accommodation  coefficient  of  energy 
is  received  as  constant.  For  the  calculations  with  the  reflection  of 
molecules  according  to  the  law  (4.11)  are  used  blocks  II  and  III. 


DOC  -  83067304 


PAGE  /  7  f 


Into  the  initial  information  of  the  version  of  calculation 
enters  the  information  about  the  surfaces  of  combination  of  which 
consists  entire/all  surface  of  the  compound?  the  signs/criteria  of 
surfaces,  calculation  of  local  flows?  the  sign/criterion  of  the 
mirror-diffused  reflection  of  molecules  from  the  wall  (into  nucleus 
0332  are  sent  zero)?  number  M*.  more  than  which  there  cannot  be  the 
number  of  reflections  of  molecule  from  the  concave  surface  of 
compound  with  the  drawing  of  one  trajectory  of  molecule?  the  angles 
of  attack  a,  and  slip  /3,  (as  and  0,  are  assigned  in  degrees)? 
parameter  H?  the  accommodation  coefficient  of  energy  a,;  the  portion 
of  diffuse  reflecting  molecules  p„?  area  and  the  diameter  of 
midsection  \  and  dn;  the  coordinate  of  the  first  x'rt  and  the  second 
x'i  of  the  points  at  which  it  can  be  placed  the  center  of  gravity? 
the  coordinate  of  the  faces  of  parallelepiped  (control  surface) 
and  xu(xu<xti)\  lx  and  L,?  Tw;  number  **pi  n*p  and  Nip, 

controlling  the  readout  of  calculation.  Are  assigned  also  numbers 
n(P),  n(a,),  n(/3„),  n(H)  and  *(*«)»  the  designating  quantities  of 
surfaces,  angles  a,,  £,,  H  and  Addresses  for  the 
introduction/ input  of  the  parameters  enumerated  above  are  indicated 
into  block  I (O) . 

Information  about  each  surface  is  designed  in  the  form  of  the 
array  of  30  numbers  in  the  sequence,  indicated  in  Table  2. 
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In  this  table  v  -  the  number  of  surface;  a«l,  if  surface  is 
flat/plane,  otherwise  of  a*0,  b*0  in  the  absence  of  the  recording  of 
local  flows  over  the  surface  v,  otherwise  of  b«-l  during  the 
calculation  of  local  flows  according  to  inside  of  surface  and  b=+l 
over  the  external  surface.  The  recording  of  local  flows  can  be 
produced  only  for  one  surface.  With  the  work  of  program  with  the 
recording  of  local  flows  into  nucleus  0330  is  sent  a  number  one, 
otherwise  -  zero.  If  limits  /3X,  /3a,  7*  and  7,  D  region  designate 
angles,  then  they  are  assigned  in  the  radians. 


PAGE 


Key:  (1).  position.  (2).  Parameter, 


naQtMerp 


roi 


Page  90. 


The  information  comprised  in  accordance  with  Table  2  about  all 
surfaces  is  introduced  into  nuclei  0420-1415  (0420-0455  about  surface 
of  1,  0456-0513  about  surface  of  2,  etc.).  In  all  it  is  possible  to 
introduce  information  about  17  surfaces  (nucleus  0420-1415). 


For  the  designation  of  the  system  of  coordinates  of  the 
assignment  of  limits  /3lf  7,,  and  7,  the  rectangular  region  D  each 
surface  is  accompanied  by  the  indicative  code  (Table  3).  The  codes 
are  joined  into  the  array  by  way  of  the  location  of  surfaces  in 
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nuclei  (0420-1415)  and  are  introduced  into  nuclei  5201-5221. 

Let  us  note  that  during  the  assignment  to  D  region  in  the 
Cartesian  and  polar  coordinates  it  is  necessary,  in  order  to  for  all 
points  of  surface  y,>0.  In  the  case  of  cone  it  is  necessary  to  reject 
a  small  vicinity  about  the  apex/vertex.  According  to  calculation  data 
the  linear  dimensions  of  this  vicinity  must  be  not  less  than  10'* 

(for  machine  M-20),  otherwise  can  occur  the  stop  of  machine. 

Calculation  according  to  the  program  can  be  produced  for  certain 
limited  quantity  of  the  parameters:  a,  (12),  /3,(12),  H(3)  and  «,  (3). 
Are  here  in  the  brackets  indicated  maximum  quantities  of 
corresponding  parameters.  During  the  calculation  occurs  the 
consecutive  sorting/excess  of  all  combinations  of  the  values  of  these 
parameters:  first  is  sorted  out  every  a#,  then  /30,  H  and  a*.  Numbers 
n(P),  n(a0),  n(/3,),  n(H)  and  *(a*)#  equal  to  a  number  of  assigned 
surfaces  a#,  /3,  H  and  ■„  are  introduced  by  the  octal  codes  with  the 
second  address  into  the  nuclei  (Table  4). 

During  the  calculation  according  to  the  value  of  parameter  H 
with  the  help  of  linear  interpolation  are  computed  values  V x,  of  the 
average  speed  of  molecules  V  and  The  interpolation  points  H,  and 

value  tfx,  V  and  T„  in  them  are  assigned  in  the  nuclei,  indicated  in 
block  1(0)  (see  numerical  information).  Number  .V,*  is  introduced  by 
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the  octal  code  with  the  second  address  into  nucleus  6615. 

The  drawing  of  particle  trajectories  through  all  faces  is 
produced  by  "portions"  NXp.  A  maximum  quantity  of  developed  particle 
trajectories  is  equal  to  n2p.  Printout  of  total  aerodynamic 
coefficients  is  produced  periodically  after  the  drawing  of  next 
"portion"  of  trajectories  Nlp’=*ktNxp  (k,»l,  2,  3  ...),  and  the 
coefficients  of  local  flows  -  after  the  drawing  of  "portions"  of 
trajectories  N*p  =  ltiNXp  (k4*l,  2,  3,  ...).  Number  N2p  must  be 

integral  multiple  of  numbers  N3P  and  NXp.  Numbers  NXp-*-Nip  are  sent 
into  nuclei  Nx  p  -  0253,  M,  — 0254,  yvt,-0256  and  Nlp  —  0333. 

Let  us  note  that  the  values  of  some  parameters 
K  •»  V,  7*00,  Tw,  Lx,  1])  are  given  in  block  1(0) . 

The  described  initial  information  is  introduced  into  the  memory 
of  machine  with  block  1(0)  in  the  following  sequence:  is  introduced 
block  1(0)  with  the  previously  specific  and  pierced  on  the 
latter/last  punch  card  check  sum  (KS);  after  the  comparison  of  KS  of 
block  I  control  automatically  is  transmitted  to  the  first  nucleus?  on 
the  command/crew,  which  is  located  in  the  fifth  nucleus,  is  produced 
the  input  of  initial  information  and  begins  the  calculation  of 
problem  from  nucleus  0006. 
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Table  3. 


CacreMa 

Koa 

KOOpAIMT  ^ 

_ Ql _ 

JesapTOBaa 0) 

0  00  0000  0000  0000 

floiapaaa  (*i) 

0  00  0001  0000  0000 

UuHHAp  h  aecKan(% 

0  00  0000  0001  0000 

C^epaaeciaa  (it) 

0  00  0000  0000  0001 

Key:  (1).  Coordinate  system.  (2).  Code.  (3).  Cartesian.  ( 
(5).  Cylindrical.  (6).  Spherical. 


Table  4. 


AApec 

Koa 

Oo 

6616 

0 

00 

0000 

«(n) 

0000 

6617 

0 

00 

0000 

rt(«o) 

0000 

6620 

0 

00 

0000 

«(?o) 

oooo 

6621 

0 

00 

oooo 

«(//) 

0000 

6622 

0 

oo 

0000 

«(«•) 

oooo 

Key:  (1).  Address.  (2).  Code. 

Page  91. 


) .  Polar 


Calculations  according  to  the  program  are  produced  in  the 
following  sequence  (in  the  brackets  are  given  the  numbers  of  the 


DOC  »  83067304 


PAGE 


nuclei  of  program) : 

1)  the  exchange  between  MOZU  [core  storage]  MB-1  and  MB-2,  the 
press/printing  the  check  sum  of  version  (0006-0063); 

2)  the  expansion  of  operating  field  (RP)  IS-2  (1500-1501); 

3)  the  translation/conversion  of  the  array  of  numbers  of  initial 
information  from  the  decimal  system  into  the  binary  (1502-1503); 

N 

4)  printout  of  initial  information  (2507-2512); 

5)  the  recalculation  of  the  coefficients  of  the  equations  of  all 
surfaces  upon  transfer  from  one  axes  Yyt  to  the  next  Xx, 

(6015-6111); 

6)  calculation  and  printout  of  the  values  of  areas/sites  A>£.,  if 
contents  of  nuclei’s  0330  is  equal  to  1.0  (5327-5712); 

7)  the  calculation  of  values 

<S £,  a0,  Po,  H,  a*,  v®/,  V®,  V,  T ®,  S®,  ft®,  2x  (V^®*),  , 

/■»(&*)•  5if.  (1514-1530,  4010-4266,  4610-4653,  6507-6513,  1541); 

8)  the  drawing  of  the  random  speeds  and  coordinates  of  particles 
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with  the  start  from  faces  S*  (1541-1753,  2005-2022); 

9)  the  program  of  the  drawing  of  the  pseudorandom  evenly 
distributed  numbers  Rx-R10  in  interval  of  0-1  (6676-6704),  of 
constant  to  this  block  (6713-6750); 

10)  the  calculation  of  the  collision  of  particle  with  the 
surface  of  compound  (2030-2536); 

11)  the  calculation  of  the  components  of  internal  normal  rtt  to 
the  surface  at  the  impact,  point  in  the  particle  and  value  l\rint 
(2537-2577); 

12)  the  calculation  of  the  components  of  particle  speed  with 
diffuse  reflection  (2403-2633,  3425-3434,  3151-3166,  2634-2654, 

1505- 1513,  2655-2660,  3435-3440),  with  mirror  reflection  (2661-2673, 
6237-6241)  and  with  the  reflection  according  to  the  law  (4.11) 
(6320-6342,  2622-2634,  6343-6463,  6165-6207,  6465-6506,  6140-6160, 

1506- 1513,  2655-2660,  3435-3440); 

13)  the  calculation  of  total  flows  (2674,  3037,  6271-6317, 
2676-2715,  3066-3104); 

14)  the  calculation  of  local  flows  (2717-3036,  6161-6164); 


DOC  =  83067304 


PAGE  “7 


15)  calculations  before  the  printout  of  total  aerodynamic 
coefficients  (1754-1771,  3112,  6530-6535,  3113-3124,  3145-3150, 
6536-6542,  4272-4306,  4667-4777,  4314-4433,  5000-5036,  4434-4450, 
5037-5060,  4451,  4520-4545); 

16)  the  command/crew  of  printout  of  total  aerodynamic 
coefficients  (4546-4547); 

17)  calculations  before  the  printout  of  the  coefficients  of 
local  flows  (the  same  commands/crews  as  in  p.  15  with  the  addition  of 
commands/crews  in  nuclei  3124-3144,  3577-3631,  3057-3062,  6562-6607); 

18)  the  command/crew  of  printout  of  the  coefficients  of  local 
flows  (3615-3616,  3627-3630,  6573-6574); 

19)  the  sorting/excess  of  particles  N2a  (4550-4554,  4654-4655); 

20)  sorting/excess  a,  /3,  H  and  a,  (4655-4633,  4555-4567, 
4452-4475,  3340); 

21)  the  restoration/reduction  of  program  and  the  introduction  of 
the  new  version  of  calculation  (3340,  3145-3147,  4571-4577, 
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0066-0076,  0006). 


Let  us  note  that  during  the  transfer  of  control  to  nucleus  3340 
at  any  moment  of  calculation  occurs  the  restoration/reduction  of 
program  and  the  introduction  of  new  version. 

To  monitor  the  correctness  of  the  calculation  of  problem  is 
possible  on  the  printouts  of  initial  information  (2507-2512)  from 
nuclei  0260-0346  and  0420-0646,  and  also  by  tracking  after  nucleus 
3521,  in  which  occurs  the  accumulat ion  of  number  of  the  played 
trajectories,  and  by  nucleus  3550,  in  which  occurs  a  cyclic  increase 
in  the  number  of  multiple  reflections  of  particle  from  the  concave 
surface  in  the  developed  trajectory. 


Page  92. 


During  the  calculation  of  areas/sites  AS.  varies  the  contents  of 
nucleus  5672.  To  the  press/printing  areas/sites  AS.  are  put  out  by 
one  array  (5704-5705)  in  the  sequence,  which  corresponds  to  numbering 
<w«(/J-l)  m+7. 

In  the  program,  the  number  10,  *6(0 — 100),  06(0—10),  76(0—10). 

The  accuracy  of  calculation  AS.  is  equal  to  0.0001.  The 

check  of  the  correctness  of  the  operation  of  program  and  machine  can 

be  also  produced  with  calculation  for  the  standard  version. 

Printout  of  total  aerodynamic  coefficients  is  produced  by  one 
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array  (4546-4547);  the  order  of  the  location  of  the  parameters  in  the 
array  is  given  in  Table  5. 

In  the  printed  array  of  the  group  of  the  parameters  are  isolated 
with  the  help  of  conditional  number  n,  of  equal  to  — ' — (-00999999999, 
and  with  the  help  of  index  *  in  the  45th  digit  of  the  binary  code 
of  the  number  (numbers  of  these  numbers  in  Table  5  are  isolated  with 
prime).  Aerodynamic  coefficients  in  the  drag  axes  are  isolated  with 
index  and  moment  coefficients  relative  to  points  x'u  and  jct\  — 

by  index  ,Ti“  in  accordance  with  one  and  two  primes. 
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The  coefficients  of  local  flows  are  put  out  to  press/printing  34 
by  arrays  of  numbers  (Table  6).  Arrays  1-11  give  the  coefficients  of 
the  average/mean  local  flows  through  areas/sites  45..  Array  12 
consists  of  11  numbers  p\,  n'+t  p'_,  ^.With  the  help  of 

the  parameters,  which  are  contained  in  arrays  12-34,  easily  are 
expressed  the  coefficients  of  average/mean  flows  along  the  sides  /i 
and  7  of  the  rectangular  region  D 


P+t 
10"  : 


(1.1) 


and  average/mean  flows  in  the  entire  surface 


(1.2) 


Varying  in  formulas  (1.1)  and  (1.2)  p  on  r,  e  and  n,  and  index 
»+•'  to  •  and  r,  we  will  obtain  formulas  for  all  averaged  flows. 

In  the  arrays  1-11  numbers  are  divided  into  the  groups  (by  10  numbers 
in  the  group)  with  the  help  of  index  •  in  the  45th  digit  of  the 
binary  code  of  a  number. 


2.  Blocks  II (M)  and  I I I (  M)  for  the  drawing  of  the  reflection  of 
particles  from  the  surface  according  to  the  law  (4.11). 


Basic  part  of  the  program  for  the  simulation  of  reflection 
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Table  5. 
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Key:  (1).  No.  of  position.  (2).  Parameter. 
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Key:  (1).  Number  of  array.  (2).  Array.  (3).  Quantity  of  numbers. 

Page  94. 

Let  the  accommodation  coefficients  for  each  surface  be 
determined  and  dependences  on  the  angle  01  (angle  between  vectors  of 
the  speed  of  the  falling/incident  molecule  and  vector  of  internal 
normal  to  the  surface).  In  the  program  are  provided  the  nuclei  for 
dispatching  six  functions  «.u*“f  (d4 ) ,  flui.j—fCflx)  and 
The  values  of  these  functions  at  points  0!*O?  15° ?  30°;  45°;  60°; 

75°;  90°  are  sent  into  nuclei  0134-0205  [see  the  numerical 
information  of  block  1(0)].  For  each  surface  is  comprised  also  the 
special  code  whose  first  address  is  the  number  of  the  first  nucleus 
of  dependence  «,=/(•,)  (for  this  surface),  the  second  address  -  number 
of  the  first  nucleus  of  dependence  /&)»  the  third  address  -  number 
of  the  first  nucleus  of  dependences, = All  these  codes  are  joined 
into  the  array  and  are  introduced  into  nuclei  5222-5242  (into  5222  - 
for  surface  of  1,  into  5223  -  for  surface  of  2,  etc.).  In  the  given 
program  values  *•  [see  the  numerical  information  of  block  1(0)] 

are  undertaken  from  [10]. 

With  block  II  are  introduced  numbers  (p.-l)  (1420),  a<0  (0313) 
in  the  implementation  of  dependence  (4.21)  and  (p,-l)  (1420), 

(0313)  in  the  implementation  of  dependence  (4.22).  Are 
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here  in  the  brackets  indicated  the  numbers  of  nuclei. 

Block  II+III  makes  it  possible  to  simulate  arbitrary  monotonic 
dependences  9,,  9,  on  the  square  of  speed  of  incident  particle. 

With  block  II+III  are  sent  the  values  of  interpolation  points  g* 
into  nuclei  6764-6774  and  respectively  the  value  of  functions  in  them 
9,  (6775-7005),  <?«  (7006-7016) ,  fc  (7017-7027) . 

3.  Block  IV  (P)  for  calculating  the  coefficients  of  the  aerodynamic 
forces,  which  function  on  each  surface. 

Printout  of  the  aerodynamic  coefficients  of  surfaces  is  produced 
periodically  after  the  drawing  of  next  "portion"  of  trajectories  Nip 
(see  the  description  of  block  I).  The  coefficients  of  surfaces  are 
put  out  to  the  press/printing  by  the  arrays  each  of  which  relates  to 
the  specific  surface  with  respect  to  the  location  of  information 
about  the  surfaces  in  nuclei  0420-1415.  Table  7  gives  the  order  of 
the  location  of  the  aerodynamic  coefficients  of  surface  in  the  array. 

After  printout  of  the  coefficients  of  all  surfaces  is  produced 
the  printout  of  total  coefficients  just  as  in  block  I. 

4.  Block  V(V)for  the  simulation  of  free  molecular  flow  in  the  concave 


Key:  (1).  position.  (2).  Parameter. 

Page  95. 

Concave  surface  is  oriented  so  that  the  plane  of  entrance  would 
coincide  with  the  plane  of  the  fourth  face  of  control  surface  (k«4), 
moreover  the  center  of  the  area  of  the  figure  of  entrance  to  the 
center  of  the  fourth  face  of  control  surface  they  must  coincide. 
Depending  on  the  form  of  the  duct/contour  of  entrance  into  nucleus 
3344  is  introduced  the  special  code,  indicated  in  Table  8. 

In  the  case  of  the  rectangular  form  of  entrance  its  duct/contour 
must  coincide  with  the  duct/contour  of  the  fourth  face  of  control 
surface,  in  the  case  of  the  elliptic  form  of  the  entrance,  inscribed 
into  the  duct/contour  of  the  fourth  face,  into  nuclei  0243  and  0244 
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are  sent  the  sizes/dimensions  of  the  semi-axes  of  ellipse,  parallel 
with  respect  to  axes  xa  and  xa. 

With  the  printout  of  total  coefficients  are  printed  parameters 
1-120  (see  Table  5),  then  they  are  printed  (4551-4552)  two  numbers: 
number  of  particle  trajectories,  played  from  the  fourth  face  and 
which  did  not  fall  on  the  entrance  into  the  concave  cavity  and  which 
fell  on  entrance.  Printout  of  the  coefficients  of  local  flows  is 
produced  just  as  in  block  I. 

5.  Block  VI  (K)  for  the  calculations  of  the  local  parameters 
according  to  the  field  of  internal  and  external  flows. 

With  the  work  with  this  block. the  drawing  can  be  produced  over 
the  entire  surface  of  parallelepiped  -  control  surface  or  only  on  the 
fourth  face  of  control  surface  (entrance  into  the  internal  cavity, 
see  block  V).  In  the  first  case  into  nucleus  3344  are  sent  zero,  in 
the  second  case  into  nuclei  3344,  0243  and  0244  is  sent  the 
information  according  to  the  description  of  block  V. 

Information  about  the  "formal"  surfaces  whose  number  must  not 
exceed  16,  is  introduced  into  nuclei  0420-1415  just  as  for  the  usual 
surfaces,  the  order  of  the  location  of  "formal"  and  usual  surfaces 
not  playing  role.  In  the  sequence  *>«1,  2,  3  ...of  all  surfaces  there 

is  the  numeration  A *1» 2, 3  •••  of  the  formal  surfaces.  In  nuclei 
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3242-3261  are  introduced  the  sizes/dimensions  of  areas  &  of  all 
"formal"  surfaces  (into  the  nucleus  3242-Sx,  in  3243-S,  and  so 
forth).  In  nuclei  5222-5242  with  the  information  described  in  blocks 
II  and  11+  III,  is  introduced  further  information  about  all  surfaces: 
instead  of  the  code  of  operation  are  introduced  zero  in  the  case  of 
usual  surface  and  number  X  (by  octal  code)  in  the  case  of  "formal". 
For  example,  if  surface  of  1  and  3  usual,  surface  is  2  "is  formal", 
then  into  nuclei  5222-5224  are  introduced  the  codes 

0  00  0000  0000  0000 

0  01  0000  0000  0000 

0  00  0000  0000  0000 

• 

In  the  third  address  of  nucleus  6137  by  octal  code  is  introduced 
a  number  of  "formal"  surfaces.  • 

Printout  of  the  results  of  calculating  the  local  parameters  of 
gas  according  to  the  field  of  flow  is  produced  periodically  after  the 
drawing  of  a  number  of  trajectories  Nip. 


GE  /f? 


npllBIO  jrOJIbHHK  00 
3ui  DC  (*0 


0  00  0001  0000  0000 
0  00  0000  0001  0000 


Key:  (1).  Form  of  entrance.  (2).  Code.  (3).  Rectangle.  (4).  Ellipse. 
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Are  printed  5  arrays.  In  arrays  2-5  to  the  press/printing  are  put  out 
55  groups  of  numbers  (in  array  2  -  group  1-15,  in  array  3  -  group 
16-30,  in  array  4  -  group  31-45,  in  array  5  -  group  46-55)  in  the 
sequence,  indicated  inTable  9.  In  each  group  of  numbers  are 
contained  the  values  of  the  corresponding  parameter  for  all  "formal" 
surfaces,  in  the  sequence,  which  corresponds  to  numbering  X*l,  2,  3 
....  For  example,  in  the  beginning  of  array  2  are  printed  all  values 
W+x,  then  all  values  W.x  and  so  forth.  The  number,  arranged/located 
in  the  beginning  of  each  group,  i.e.,  the  value  of  the  corresponding 
parameter  for  A*l,  is  isolated  with  the  help  of  index  •  in  the  45th 
digit  of  the  binary  code  of  a  number. 
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Apparently,  in  certain  cases  it  is  useful  to  use  this 
calculation  procedure:  through  each  point  of  field  of  flow  whose 
vicinities  are  investigated,  to  carry  out  through  three  mutually 
perpendicular  "formal"  surfaces,  and  the  results  of  calculation  to 
average  according  to  the  formulas: 


<  X»3 k  .  i-3* 

s  ^  r‘4  I  & 


i=3*-2 

,  X=3»  «  X-3* 

Z  ^  =  t  Z 

3  X»3*-2  3  '  -  - 


0  X=3*_  2 

fl 

J  X=3*-2 


(5.1) 


In  array  1  are  printed  h„,  Tk,  p'kV *  for  all  k**l,  2,  3,  ...  the 
points  of  the  field  (numbers  with  numbers  1-4  designate  n\,  T\,  p\,  V\\ 
numbers  5-8  they  designate  th,  T,,  pi,  and  so  forth). 


Printout  of  the  coefficients  of  total  flows  is  produced  just  as 
in  block  V,  after  printout  of  the  coefficients  of  the  local 
parameters  in  the  field  of  flow. 
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Key:  (1).  the  group  of  numbers.  (2).  Parameters  of  array. 
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6.  Block  VII (p)  and  VII (7)  for  the  detailed  calculation  of 
coefficients  of  the  local  flows  through  the  surface. 
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In  block  I  during  the  calculation  of  local  flows  the  rectangular 
region  D  is  divided  on  100  equal  pads  with  the  help  of  18  straight 
lines,  carried  out  through  each  side  D  region  in  parallel  to  its 
sides.  In  certain  cases  can  be  required  the  division  of  region  on  the 
finer/smaller  pads. 

In  blocks  the  VII  and  VIII  rectangular  region  D  is  divided  into 
100  equal  bands  with  the  help  of  99  straight  lines,  carried  out  in 
parallel  to  the  side  7  of  region  D  (block  VII)  or  in  parallel  to  the 
side  /3  of  region  D(block  VIII).  Let  us  designate  the  areas  of  bands 
AS.,  moreover  the  number  of  area/site  »  =  p€(0— 100)  for  block  VII  and 
<#  =  y€(0 — 100)  for  block  VIII.  This  division  D  field  proves  to  be  useful 
during  the  detailed  study  of  local  flows,  constants  along  any  one 
side  D  region. 

With  the  printout  of  the  coefficients  of  local  flows,  which 
occurs  periodically  after  the  drawing  of  a  number  of  trajectories 
/V,„,  are  printed  23  arrays  of  numbers.  The  arrays  1-11,  which  contain 
according  to  100  numbers,  and  Massey  12,  who  contains  11  numbers, 
contain  information  about  the  average  coefficients  of  the  local  flows 
through  areas  AS,  and  through  entire  surface;  these  arrays  are 
analogous  described  in  block  I  and  given  in  Table  6. 
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Let  us  note  that  here  index  w«/S  during  the  use/application  of 
block  VII  and  u-y  during  the  use/application  of  block  VIII. 

In  view  of  the  smallness  of  areas/sites  A 5.  are  possible  the 
considerable  fluctuations  of  the  values  of  the  parameters  upon 
transfer  from  one  area/site  to  another.  Therefore  it  is  convenient  to 
also  have  the  averaged  values  of  the  parameters  on  several 
areas/sites  AS..  Arrays  13-23,  which  contain  according  to  20  numbers, 
give  the  averaged  values  of  the  parameters,  which  are  contained  in 
arrays  1-11,  on  five  areas/sites  AS„;  for  example,  the  parameters, 
which  contain  in  array  13,  are  obtained  according  to  the  formula 

P+»  —  ~*  V  />+.,  (6.1) 

5  mmSi-4 

where  —  parameters  of  array  1,  and  index  k  takes  values  of  k*l, 
2,  3,  ...,  20;  for  array  14  parameters  are  obtained  analogously  from 
array  2,  etc. 

Printout  of  the  coefficients  of  the  total  of  flows  is  produced 
just  as  in  block  I. 

7.  Block  IX  for  the  simplified  assignment  of  the  parameters  of 
surfaces. 


DOC  «  83067304 


PAGE 


Applying  block  IX,  it  is  possible  to  decrease  and  to  simplify 
information  about  the  following  surfaces,  most  widespread  in  practice 
(Fig.  1-6  and  Table  10);  parallelogram,  trapezoid  (triangle), 
elliptical  ring,  elliptical  cylinder  and  cone,  ellipsoid. 


Table  10  gives  the  order  of  the  location  of  the  parameters  of 
the  information,  which  contains  12  numbers,  for  each  surface. 
Information  about  all  surfaces  whose  number  is  not  more  than  17,  is 
introduced  into  nuclei  0420-0733  (into  nuclei  0420-0433  about  surface 
of  1,  0434-0447  about  surface  of  2,  etc.). 


II  0  0  !-fV3«  !—?/*■  l-t/J»  l-t/w 

I*  0  -1.0  0  0  a  !-•/* 


Page  100. 

In  Table  10  parameters  1-9  are  the  coordinates  of  points  A,  B 
and  C,  arranged/located  relative  to  surface  in  accordance  with  Fig 
1-6  and  figures,  given  in  Table  10.  Coordinates  x\,x\  and  x-  are 
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assigned  in  the  system  of  coordinates  Xxt.  The  angles  of 

longitude/length  ?*  and  <p  and  the  angle  of  latitude  9,  with  the  help 

of  which  is  selected  the  piece  of  surface  in  question,  are  changed  in 

1 6  <0-180°). 

the  ranges:  ?* 6 (0—360°),  ? 6  (0—360°),  A  Angles  <p  and  <?*  are  counted  off 
from  AB,  angle  9  is  counted  off  from  the  line,  perpendicular  to  plane 
BAC. 


If  necessary  for  the  recording  of  the  local  particle  fluxes, 
impulse/momentum/pulse  and  energy  into  nucleus  0330  is  sent  the 
number,  equal  to  the  reference  number  of  simple  surface  and  which  has 
plus  sign  during  the  recording  of  flows  over  the  external  surface  and 
minus  sign  during  the  recording  of  flows  over  the  internal  surface. 
Let  us  note  that  the  recording  of  local  flows  can  be  produced  only 
over  one  surface. 

Each  surface  is  accompanied  by  the  indicative  codes  (Table  11), 
which  are  joined  into  the  array  by  way  of  the  location  of  surfaces  in 
nuclei  (0420-0733)  and  are  introduced  into  nuclei  5201-5221. 

In  the  case  when  there  are  only  simple  surfaces  or  together  with 
the  simple  ones  there  are  complicated  surfaces  information  about 
which  can  be  been  given  only  on  Table  2,  the  information,  which 
composes  version,  is  introduced  in  the  following  order:  with  the 
address  code  (AK)  0420  information  only  about  all  simple  surfaces. 
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comprised  on  Table  10;  with  AK5201  and  AK5222  -  codes  for  all 
surfaces  (simple  and  complicated);  AK6616  -  number  of  all  surfaces 
(simple  and  complicated);  with  AK  1364  -  the  number  only  of  simple 
surfaces;  then  is  placed  all  remaining  information;  the  necessary 
blocks,  moreover  block  IX  is  arranged  always  for  information 
processing  about  the  simple  surfaces:  after  block  IX  are  introduced 
zero  into  nucleus  1345  in  the  case  of  the  presence  together  with  the 
simple  ones  of  complicated  surfaces,  otherwise  it  is  not  introduced; 
then  is  placed  punch  card  with  the  check  sum,  equal  to  zero;  without 
the  address  code  is  placed  the  information  about  the  complicated 
surfaces,  comprised  according  to  Table  2;  version  is  finished  with 
punch  card  with  the  check  sum,  equal  to  zero.  Information  for 
calculating  the  local  flows  according  to  the  complicated  surface  is 
introduced  according  to  the  description  Section  of  1  Chapter  II.  The 
examples  of  assignment  to  information  are  given  in  appendix  2. 

Let  us  make  some  observation  about  the  algorithm  of  the  program 
of  block  IX,  according  to  which  are  calculated  the  parameters  of 
surfaces,  described  in  block  I  (see  Table  2). 


Cosines  a},  are  determined  from  the  formula 

a]i  =  a'Jma’mi(l,  m=l,2, 3);  (7.1) 


aa  =  a3 iXa’u,  a. 
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Here  a*,—  cosines  of  the  angles  between  axes  xt  and  axes 
y'm  of  the  intermediate  coordinate  system;  ajm  —  cosines  of  the 
angles  between  axes  y,  and  y„  (see  Fig.  1-6). 
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Key:  (1).  Surface.  (2).  Code.  (3).  Parallelogram.  (4).  Trapezoid 
(5).  Elliptical  ring.  (6).  Elliptical  cylinder  and  cone.  (7). 
Ellipsoid. 
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The  sides  of  parallelogram  AB  and  AC  it  is  placed  in  the 

coordinate  planes  y,  y3  andV^V, .  After  conversions  we  find 

an  =■  \y l  —  cos  0,  /cos*}; 

„•  _  {  —cos I  I* cos 4* ( 1  —  cos  . 

°J 2  }/  1  -+•  cos  ’  V  \  +  COS  t  '  L  *  +  COS  ‘■J*  J  i 

1,2  [  1  —  COS  <!f 


P 

1  COS  <|*  L 


3 


.5’ 


DOC  »  83067304 


PAGE  £tf/ 


Limits  D  region: 

—  0;  —  cos*;  i,  -  0;  i,  -  AC  /I  -  cost-  (7-6) 

The  side  of  trapezoid  AB  it  is  directed  in  parallel  to  axis  -y x . 
-  y\>  -  angle  between  the  plane  of  trapezoid  and  the  plane  7i7 *• 

From  Fig.  2  we  note 

a’j\  —  {— 1,  0,  0};  a/j={0,  —  cos<b,  sin<|*};  o>3—  {0,  sin +,  cos'}*}.  (7.7) 

Equation  of  the  plane  of  the  trapezoid: 

sin  <j» +>»  cos  •(»  —  Af  A/cos  ^  =0.  (7.8) 

Coordinates  of  point  Y,  which  coincides  with  point  N, 

MA-- ACT- AB-ACI(AB-b),  j  g) 

MN=—  ati (air AC°)- M A,  NA=MA~MN,  > 

Limits  D  region: 

P,  -  ?,  -  arcsin  (Z  -f?A%  {AB°-  NA°)  <  0, 

Ps  =  =  arcsin  (Zr  NB3),  (AB°  •  NB°)  <  0; 

P,  =  =  it  —  arcsin  (aT,-  NA°),  ( AB°-NA° )  > 

&  =  <?,=*-  arcsin  (<£r NB°),  ( AB°  ■  NB°)  > 

71=^11-0.  ^-yn-ian^AC. 

r 
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Here  A?0,  NA°,  n!°,  aS°  -  unit  vectors. 


For  the  elliptical  ring  of  axis  y'm  it  is  directed  so  that 

aji  -  {cos-K  0,  -  sin*},  a]7  -  {0, 1,0},  a;a=(sln  0,  cos  +},  cos  <|»  -  ~  .  n  j 

AB 

Equation  of  the  plane  of  the  ring: 

yt  sin  ij»  -j-y,  cos  <j» — (Aflsinf-|-  l)cos<{»— 0.  (T.  14) 


Coordinates  Xo»  of  point  Y: 
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*o  i =•*/'  -  <h  t  +  1) 


(7.15) 


The  limits  D  region  are  calculated  from  the  formulas: 

f*  =  *.ir/2  +  A«p*,  A«p*<  — ,  *  =  0,  1,  2,3,  4; 

Pi  -?!=  0; 

*  *0 

k=**<p,  =  *  —  npw  | A?* | <0,001; 

*»  • 

P,-?a-afctg^,  A-0; 

P,=?3-*  +  arctg^,  *=i,2; 
k^-21  +  arctg^l  *-3; 


(7.16) 

(7.17) 


<o 

npH  |  A®*|>0,  001  (7.17') 


Key:  (1).  with. 


Tl  —  ACgjj,  IjmmAC. 


(7.18) 


For  the  elliptical  cylinder 


DOC  »  83067304 


PAGE 

y\  y\  03 

AB2  AC-  —  ^  (ypaBHeHHe  noBepxHocTH); 

a'jm  -  *Jml 
X0  i  **  xi » 

=0,  p,  =  f,  Tl  =  0,  T,  =  Z.. 

Key:  (1).  the  equation  of  surface. 


(7.19) 


(7.20) 


For  the  elliptical  cone 

y\  y\  y*  CD 

— ~H2~^  (TPaBHeHHe  noBepxnocTH)  (7.21) 

Key:  (1).  the  equation  of  surface. 


where  H»AC*L* (AC-a) ; 


ajm  i« 


xoi  —  xi  <hi  | 

Pt—  0,  &  =  <?,  7,  —  -f  10~\  T|  =  //.  I 


1 


(7.22) 


For  the  ellipsoid 


y\  ,  y\  ,  y\  . 

"45^  +  +  "C*"  “  *  rypaBHeHHe  noBepxHoctH); 


a/m 


-<0  <  =  x'{, 

—  o,  &  =>  ?,  Tl  =  0.  T*  = ( 


(7.23) 

(7.24) 


Key:  (1).  the  equation  of  surface. 


8.  Calculations  in  the  flow  of  Newton  and  in  the  flow  of  light/world. 


DOC  -  83067304 


PAGE 


For  the  calculations  in  Newton's  flow  it  is  necessary  to 
introduces® ;§>  1,  A^=l,p,*0.  In  nucleus  0251  to  send  number  0.01.  The 
corresponding  coefficients  are  calculated  from  formulas  (9.1)  and 
(9.2)  Chapter  I . 

For  the  calculations  in  the  flow  of  light/world  it  is  necessary 
to  be  given  S®^  1,  p0»ka  (portion  of  the  light/world  absorbed  by 
surface) ,  ^,->1  and  to  introduce  into  nucleus  3637  the  command/crew 

0  00  0051  0000  6520 

and  into  the  nucleus  0251  number  0.01. 

Page  103. 

The  designed  coefficients  in  the  flow  of  light/world  are 
referred  to  "velocity  head"  . 

* C9 

During  calculations  in  the  flow  of  Newton  and  in  the  flow  of 
light/world  into  nucleus  0332  are  sent  zero  -  sign/criterion  of  the 
mirror-diffused  reflection. 

9.  On  the  sequence  of  the  introduction/ input  of  blocks. 


During  conducting  of  calculations  with  the  use/application  of 
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several  blocks  whose  possible  combinations  are  given  in  Table  12, 
blocks  are  introduced  after  block  I  and  version  of  calculation  in  the 
order,  which  corresponds  to  an  increase  in  the  number  of  block.  In 
Table  12  blocks  le  and  I*  designate  block  I  for  calculating  the 
coefficients  of  total  (without  the  local  ones)  or  local  (and  total) 
f lows . 

10.  On  the  sensors  of  0-1  pseudorandom  numbers  R^Rn  evenly 
distributed  in  the  interval. 

For  the  work  of  program  it  is  necessary  to  develop  ten  numbers 
Ri“Ri,.  All  ten  programs  for  generating  these  numbers  are  obtained  of 
one  program  (No  1)  of  obtaining  pseudorandom  numbers,  given  in  [26]. 
Were  preliminarily  calculated  the  constants  of  sensor  d+1,  d+2  and 
c+1  after  obtaining  i*2*10*  of  pseudorandom  numbers.  Here  i*l,  2,  3, 
...,  10.  The  obtained  constants  are  given  in  the  nuclei  of  6713-6750 
blocks  I.  Turning  to  the  program  of  obtaining  random  numbers 
6676-6704  is  produced  at  the  value  of  index  register  (RA) ,  equal  to 
the  number  of  address,  in  which  is  placed  constant  d+1.  For  example, 
during  the  turning  to  the  program  of  sensor  6676-6704  at  the  value  of 
RA,  equal  to  6713,  constants  d+1  and  d+2  for  the  work  of  sensor  are 
taken  from  nuclei  6713  and  6714,  and  next  pseudorandom  number  is 
obtained  in  nucleus  6715. 
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For  the  work  with  another  sensor  it  is  necessary  to  introduce 
the  new  program  of  sensor  into  nuclei  6676-6703  and  new  constants 
into  nuclei  6713-6750. 

In  appendix  2  in  example  3  are  given  the  program  and  the 
constants  of  another  version  of  the  sensor  of  pseudorandom  numbers 
Initial  sensor  with  period  2,J-1  was  undertaken  from  work  [26] 
(program  No  5).  Constants  for  all  ten  sensors  in  example  3  were 
obtained  after  drawing  i«6*10*  (i»lf2,3,  ...»  10)  pseudorandom 
numbers. 
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Table  12. 


Ic 

Ij 

II 

II  +  III 

IV 

V 

VI 

VII,  VIII 

IX 

lc 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

I. 

+ 

— 

+ 

+ 

— 

+ 

— 

+ 

4 

II 

+ 

+ 

— 

+ 

+ 

+ 

+ 

+ 

+ 

II  +  HI 

+ 

+ 

+ 

— 

+ 

+ 

+ 

+ 

+ 

IV  | 

+ 

— 

+ 

+ 

+ 

— 

— 

+ 

V 

+ 

4* 

4- 

+ 

-4- 

— 

— 

+ 

+ 

VI 

+  • 

— 

+ 

+ 

— 

— 

— 

— 

+ 

VII,  VIII 

+ 

+ 

+ 

+ 

— 

+ 

— 

— 

+ 

IX 

+ 

+ 

+ 

+ 

+ 

+ 

4* 

+ 

— 

Page  104. 


11.  Examples  of  calculation. 


Let  us  consider  the  examples  of  assignment  to  initial 
information  for  the  calculation  according  to  the  universal  program  of 
aerodynamic  coefficients,  local  flows  and  parameters  of  gas  for  some 
bodies  of  the  complex  form  whose  diagrams  are  given  in  Fig.  7-9.  The 
assignment  of  numbers  and  commands/crews  for  the  calculation  in  these 
examples,  and  also  the  results  of  calculation  are  given  in  appendix 


For  the  compound  (see  Fig.  7.)  consisting  of  simple  surfaces  of 
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1-6  (1  -  parallelogram,  2  -  cone,  3  -  cylinder,  4  -  trapezoid,  5  - 
ring,  6  -  sphere)  and  complicated  surface  of  7  (paraboloid  of 
revolution),  it  is  necessary  to  calculate  aerodynamic  coefficients 
for  each  surface  of  1-7  with  the  diffuse  reflection  of  molecules  with 
the  complete  accommodation  and  with  the  reflection  of  molecules  in 
accordance  with  law  (4.11).  The  order  of  the  location  of  information 
in  the  version  and  their  numerical  values  are  given  in  example  1  of 
appendix  2. 

In  nuclei  0420-0527  is  placed  the  information  about  the  simple 
surfaces,  comprised  according  to  Tables  10?  then  into  the  appropriate 
nuclei  are  placed  the  angles  of  attack  a,,  of  slip  0,,  parameter  H  in 
terms  of  value  of  which  with  the  help  of  linear  interpolation  are 
calculated  values  V  and  Too,  accommodation  coefficient  <*„, 
portion  p,  of  the  diffuse  reflecting  molecules,  Su,  d„,  coordinate 
and  *2,  the  faces  of  parallelepiped  (control  surface),  coordinate 
jc'4  and  x’rl  of  two  points  X  and  0  relative  to  which  must  be 
designed  the  moment  coefficients,  the  temperature  of  surface  Tw,  of 
number  /V,p,  iVl9,  Nif,  Nif,  the  signs/criteria  (nuclei  5201-5207)  of  all 
surfaces  of  1-7  (simple  and  complicated),  number  AU,  a  number  of  all 
simple  and  complicated  surfaces  (nucleus  6616)  and  only  simple 
surfaces  (nucleus  1364)  and  a  number  of  remaining  parameters  (nucleus 
6617-6631).  Then  are  placed  blocks  II,  IV  and  by  IX. 
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In  nucleus  1345  are  sent  zero  in  the  case  of  the  presence 
together  with  the  simple  ones  of  complicated  surfaces,  otherwise  it 
is  not  sent.  After  punch  card  with  the  check  sum,  equal  to  zero,  is 
placed  the  information  about  complicated  surface  of  7,  comprised 
according  to  Table  2. 


Fig.  7. 


a  5 


v- 


Fig.  8. 


Fig.  9- 
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Page  105. 

The  introduction/ input  of  version  {example  1)  occurs  into  two 
stages:  first  is  introduced  the  first  part  of  the  version 
(command/crew  from  nucleus  0006),  then  the  second  (loading  order  from 
the  nucleus  of  1350  blocks  IX).  In  appendix  2  is  given  the  array  of 
the  aerodynamic  coefficients  of  surface  6,  obtained  first  during 
the  calculation  without  block  II  (diffuse  reflection  of  molecules), 
and  then  with  block  II  (reflection  in  accordance  with  the  law  (4.11). 

In  example  2  is  examined  the  problem  of  calculating  the  local 
flows  along  the  height/altitude  7  according  to  inside  of  the  cylinder 
(see  Fig.  8)  with  those  determined  of  value  H,  v JO,  V  and  T sent 
in  the  appropriate  nuclei  (see  also  the  numerical  information  of 
block  I).  Were  used  the  following  blocks:  V  -  for  the  drawing  of  the 
start  of  the  molecules  through  the  round  entrance  and  block  VIII  - 
for  the  detailed  recording  of  local  flows  on  the  areas/sites,  which 
are  obtained  during  the  division  of  height/altitude  7  into  one 
hundred  parts.  Is  given  an  example  of  the  calculation  of  the  flows  of 
part  n,k  (k-1,  2,  3,  ...,  20)  of  the  multiple  reflections  on  the 
areas/sites,  which  are  obtained  during  the  division  7  into  twenty 
parts,  moreover  report  k  is  conducted  from  the  basis/base. 
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In  example  3  are  given  the  results  of  calculating  some 
parameters  of  gas  at  the  flow  through  the  air  intake,  which  is  the 
element/cell  of  the  hyperboloid  (see  Fig.  9). 

The  results  of  some  interesting  calculations  with  the  help  of 
the  universal  program  are  given  also  in  works  [27-29]. 

Was  thoroughly  checked  the  method  of  drawing  of  the  random 
parameters  of  particles.  From  given  for  an  example  Fig.  10  it  is 
evident  that  the  agreement  of  histograms,  obtained  with  the  drawing 
according  to  the  method  of  Baird  of  random  values  [see  density 

function  (3.3)],  also,  with  drawing  of  speed  with  the  diffuse 
reflection  [see  formula  (4.10)],  with  the  theoretical  data  /(?*,*)  and 
fr (K) “ 2 ;'3  exp(—  V?)  is  good  even  with  comparatively  small  numbers  of 
played  trajectories.  In  Fig.  10  N,  -  a  number  of  diffuse  reflected 
particles;  by  lines  they  are  designated  the  results  of  calculations 
according  to  precise  formulas,  by  points  -  according  to  the  Monte 
Carlo  method. 


4  2  3  4ft 


4  2  3  4ft 

Fig.  11. 


4 


a  a  4^d* 
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We  will  designate  aerodynamic  coefficients  in  the  drag  axes  by 
index  "0"  on  top.  For  the  positive  reference  direction  of  drag  let  us 
take  direction  V**,  and  for  the  projection  of  forces  on  body  axis  Xxx 
-  opposite  direction  of  axis  Xxw  i.e.,  in  the  manner  that  this  is 
accepted  in  aerodynamics. 

The  coordinates  of  the  point,  relative  to  which  is  designed  the 
moment/torque,  we  will  designate  <.  and  moment  coefficients 
relatively  X,*—  by  subscript  ,t*. 

Fig.  11  gives  the  results  of  the  calculations  of  total 

3 

aerodynamic  coefficients  for  ellipsoid  2-y*ara“’T  at  t*le  values  of 
parameters  a^lO,  a,*3,  a, -2,  x0_,»0,  S.265,  7®  —  1404°K, 

7’»  —  320oK,  **  =p,»l,  xt  i=  {lQ;  0?  0  !  =  rf,,  =  2a3.Here  (j, 

i*l,  2,  3)  -  unit  matrix.  Fig.  12  compares  the  coefficients,  obtained 
by  the  Monte  Carlo  method,  with  the  calculated  ones,  obtained  by 
calculating  the  integrals  over  the  surface.  From  the  given  results  it 
is  possible  to  conclude  that  an  error  in  the  method  of  Comte-Porta 
for  total  coefficients  c]  and  mlT  already  when  jV,—  10*  is  within 
the  limits  of  3-5%,  for  obtaining  the  coefficient  c\  with  the  same 
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error  it  is  required  ^=(3-5)  by  10s  and  it  is  more. 

Were  produced  also  calculations  for  the  compound,  which  consists 
of  two  ellipsoids,  distant  behind  each  other  at  certain  distance. 

3 

First  ellipsoid  ai*50J/*'  a,=10,  x0im*o,  the  second 

/»! 

ellipsoid  ax^lO,  a,=3,  a, *2,  jcfli  =  {0;  12.5;  -5  \.  For  both 

ellipsoids  — 8,,.  Remaining  parameters:  Soo  =  8,265,  1404°  K, 

r.  -  320°  K,  *.-p0-  1,  2  *(3  +  5 1/90),  4,-10,  jct  ,  —  { 10;  0,  0). 

•  •  •  «  _  _  _  _ _ _ 

In  Fig.  13  results  of  this  calculation,  in  which  are  considered 
multiple  collisions  of  particles,  they  are  equal  with  the  results  of 
the  approximate  computation  of  coefficients  with  the  help  of  the 
integrals  in  the  surfaces.  During  the  calculation  of  integrals  in  the 
surface  the  local  coefficients  of  tangential  and  normal 
impulses/momenta/pulses  are  considered  equal  to  zero,  if  vector  is 
opposite  Ko  and  carried  out  from  the  point  in  question,  it 
intersects  concave  surface,  otherwise  they  are  computed  from  the 
formulas  for  the  single  area/site.  An  error  in  this  calculation  will 
be,  obviously,  small  with.S«I>l,  the  diffuse  reflection  of  particles 
and  when  Tw  order  7**.  Fig.  13  shows  a  good  agreement  of  results  of 
calculation  in  limits  of 
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Pig.  12. 


Monte  Carlo  method.  (2).  Calculation  of  integrals 


dmi MmS  Mooti-JUms  ■,  CD 

om—  Mover  Mnnugi «  mim  im— «  (.SLJ 


Fig.  13. 


Key:  (1).  Monte  Carlo  method.  (2).  Calculat 
taking  darkening  into  account. 
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The  results  of  calculations  for  the  concave  part  of  the 
hemisphere  confirm  dependences  (5.7)  (Fig.  14)  and  with  the  error 
less  than  1%  with  JV,>  1000  coincide  with  the  calculation  according 
to  the  formula 

f,  -2  +  1,06  /5/S..  (41.3) 

obtained  in  work  [13]  for  the  hypersonic  approximation/approach  (Fig. 
15).  Fig.  14  and  15  results  of  calculations  by  the  method  of  Monte 
Carlo  depict  as  small  circles,  and  according  to  formula  (11.3)  -  by 
solid  lines. 

Fig.  16  gives  the  results  of  calculating  the  dispersion  of  the 
drag  coefficient  c[  and  lift  c't  during  the  flow  around  single  plate 
of  the  flow,  normal  to  its  surface.  Here  when  =  0  (cold  surface)  is 
valid  the  formula 

D  [dl  -  8  S£*  x"‘  (S* )  ( exp  (-  si)  {( 1  +  Si)/2  +  M?/2]  + 

(1  +  erf S« )  [/5s»  (3/2  +  sL)i2  -  (1/2  +  Si)  +  V*S»  Af?/2]}, 

D  |ftj «  D [cl]  =  2/Si,  K — *«  Mp = -{- V^(l  +  erf S.)  x~'  (S.)/ 2  (11.4) 

or  when  S*^  1 

D  [rtl  =  D  (*;i  -  D  (<*]=?  2  S£*. 


(11.5) 
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From  Fig.  16  it  is  evident  that  when  r.—O^the  calculations  by 
the  Monte  Carlo  method  with  an  error  in  less  than  1%  will  be 
coordinated  with  the  calculations  according  to  formula  (11.5). 

The  author  is  sincerely  grateful  to  M.  N.  Kogan  for  the  theme 
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APPENDIX  1 


UNIVERSAL  PROGRAM  FOR  STUDYING  INTERNAL  AND  EXTERNAL 
FREE-MOLECULAR  FLOWS  NEAR  AN  ARBITRARY  GROUP  OF  COMPLEX  BODIES 
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oafi  ooo  mi  oooo  i*ao 
•a9l  004  02*9  0000  1*21 
•192  09*42**  9109  l*aa 

•29*  094  0090  14*0  0001 
•299  09}  0001  7722  0001 

•297  01*  9170  9449  9*99 

•241  019  9990  0000  0000 

•243  001  7741  991*  991* 
•2««  001  1*74  9907  9907 
•249  094  0000  2949  0000 
•20«  001  7741  9914  9914 
*247  001  1*7*  9912  9912 
•270  094  0000  2949  0000 
4271  019  9990  0000  0000 
•272  074  0000  *277  0000 
•279  002  9911  4921  9911 

4274  092  0000  000*  0000 

•279  014  4249  904*  4902 

427ft  000  0000  0000  0000 

4277  002  9912  4921  9912 

•900  09ft  0000  9040  0000 

•901  902  9927  1*22  9927 

4909  019  9990  0000  0000 

•90*  074  0000  4910  0000 
•909  001  1*7*  999*  999* 
*904  092  0000  0000  0000 
•907  014  3170  2474  2713 

4910  019  3990  7722  0000 

•311  074  0000  77l7  0000 
4312  092  0000  0017  0000 
•313  ol*  2717  2474  2719 
*91*  *19  3990  7722  OOOO 
•319  034  OOOO  4290  0004 
U}«  *«I4  0417  0400 

•917  014  9109  3044  4902 

•320  0>2  OOOO  0009  OOOO 

4321  909  00*7  l*l7  1*23 

4322  112  OOO*  4321  0001 
4923  702  0030  1*20  00*7 
•92*  1*2  0010  4323  0001 
•329  009  1*17  1*1?  1*31 
•924  002  0292  1*91  1*32 
•327  019  7741  1*31  OOOO 
4330  034  OOOO  4333  OOOO 
•931  0**  1*92  OOOO  1*32 
•332  09*  0004  2*l7  OOOO 
•393  44*  0099  0092  1*20 


433*  0*9  1*20  1*20  0099 

4339  001  7741  0099  0099 
433*  o«*  0039  OOOO  0099 
4997  00*  7741  0099  0094 
43*0  009  1*20  0034  0099 
43*1  002  OOOO  0099  0099 
49*2  094  0004  2422  0097 

43**  o 9*  00**  1**0  1*97 
43*9  072  OOOO  1*97  00*0 
49*4  *00  OOOO  OOOO  1*20 
43*7  *00  0001  OOOO  1*21 
*390  *00  0002  OOOO  1*22 
*391  **00  0003  OOOO  1*23 
4992  *94  OOO*  *210  1*2* 
•393  014  *33*  7901  7010 

4999  079  1*92  1*20  1*99 
4994  004  OOOO  OOOO  OOOO 
4997  019  4399  7721  4399 

4340  072  0040  1**0  OOOO 
*341  01*  4342  *3*4  *394 

4342  013  *333  7721  *393 

4343  09*  Oil*  1**0  1*37 
434*  01*  *349  43*9  4394 
*343  033  4333  7721  *333 
434*  033  4333  7721  *393 
494?  009  1*17  OOOO  1*20 
49?0  019  1**3  7741  0000 
•371  074  OOOO  4732  439* 
•372  094  9991  4212  4*27 

4373  009  39*7  1*99  1*22 
437*  002  77*1  1*39  1*21 
4379  009  1*74  1*21  1*21 

4374  009  774*  1*21  l*2l 
*377  001  l‘tl  1**2  1*21 
4*00  0**  1**1  OOOO  4*27 

4*02  009  2401  0313  1*22 
4*03  014  4*0*  7501  7*iO 
4*0*  079  1*22  0003  1*23 
4*09  009  7742  33*7  2*92 
4*0*  00*  7794  7742  2*93 
4 * 0 7  009  2*33  39*7  2*93 
4*10  0  9*  2*39  2**0  **09 
4*11  OOO  l 3*4  9000  1*29 
f Vt 2  OOO  1973  OOOO  1*1* 
6*19  900  1*04  OOOO  1*29 
4*1*  OOO  1*22  OOOO  1*2* 
4*19  OOO  1*24  0009  1*27 
4*14  OOO  39*2  OOOO  1*30 
4*17  OOO  397o  OOOO  1*31 
4*20  032  OOOO  9000  OOOO 
4*21  900  3724  OOOO  1*32 
4*22  112  0009  4*21  OOfll 
4*23  300  3712  OOOO  00*2 
4*2*  900  4107  OOOO  37i2 
4*29  112  0010  4*23  0001 
4*2*  013  l94o  *120  19*0 
4*27  01*  1 9 7*  4{20  1979 
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44  So 

013 

1611 

6129 

1*11 

43?* 

300 

ooco 

00  30 

3  000 

4431 

000 

6124 

9000 

1  96* 

*323 

300 

9000 

00  03 

3  3  0  C 

6432 

000 

*129 

0000 

1973 

*926 

300 

OOCO 

0000 

0  0  0  0 

6433 

000 

6126 

3000 

160* 

4327 

300 

0000 

0000 

0  0  0  0 

6  4  3  4 

oon 

*127 

3000 

1*22 

«33o 

036 

0000 

6333 

0900 

6439 

300 

*130 

3000 

3562 

*331 

330 

3012 

00  0  1 

0  323 

6436 

000 

*131 

3000 

3  3  ?o 

4932 

3  70 

000  1 

6331 

OOOC 

643? 

3  00 

0000 

0000 

1*22 

*933 

01* 

693* 

00  06 

0922 

6440 

099 

1*20 

7744 

oo*i- 

*934 

000 

0000 

0030 

3322 

6441 

009 

30*1 

30*1 

00*2 

*939 

31* 

6936 

3113 

3130 

64*2 

001 

7764 

30*2 

0  0*2 

633* 

300 

3337 

0030 

3150 

64*3 

0** 

3042 

3000 

00*2 

*337 

0  3  0 

0*11 

2131 

2*03 

64*4 

001 

0041 

00*2 

0  3  77 

*9*0 

0  ?0 

2101 

6337 

0  300 

6**9 

002 

007? 

1*20 

03*1 

65*1 

316 

*3*2 

131* 

4*35 

64*6 

009 

0041 

30*1 

00*1 

49*2 

03* 

0000 

*272 

4*33 

64*7 

002 

3003 

30*1 

00*1 

69*3 

o3a 

0012 

*030 

*11* 

6*30 

016 

*431 

733  1 

7410 

*9** 

o?o 

*300 

63*3 

ooco 

6  4  3  1 

0  79 

00*1 

300  3 

00*2 

49*9 

000 

0900 

0030 

0009 

6432 

000 

0000 

3000 

0000 

*34* 

000 

9000 

0000 

0000 

6*93 

003 

0042 

00  77 

0  0*3 

69*7 

000 

3720 

0000 

004  1 

6*3* 

000 

1*20 

3000 

0041 

6990 

000 

3721 

0000 

3728 

6*39 

016 

6496 

*446 

6492 

6391 

000 

3722 

0000 

3  72  1 

6*36 

000 

OOOO 

3000 

6492 

6332 

336 

ro*i 

13*3 

3722 

6*37 

01* 

6460 

7331 

7613 

6333 

309 

330* 

330* 

1*32 

6*60 

073 

1*20 

03H 

00*1 

6934 

372 

cooo 

0000 

OOCO 

6*61 

001 

77*1 

30*1 

00*1 

6939 

109 

coco 

0000 

9060 

6*62 

003 

3041 

3*1* 

00*1 

6996 

112 

C  0 1  2 

6999 

oool 

6*63 

396 

3041 

61*3 

0077 

6537 

33* 

0003 

6331 

OOOC 

6*64 

000 

0000 

0000 

0000 

6560 

300 

COOO 

0000 

0000 

6*69 

002 

1*21 

0231 

3727 

6  561 

300 

G  0  0  3 

0000 

0000 

6*66 

302 

0000 

0231 

3730 

6962 

300 

696* 

0000 

*  *  4  * 

6*67 

000 

0000 

0030 

0000 

6963 

372 

0000 

cooo 

0300 

6  4  ?o 

002 

00*1 

372* 

3791 

*5*4 

*52 

cooo 

9300 

*977 

'6*71 

001 

3291 

0291 

3732 

*9*9 

*32 

cooo 

0000 

*371 

6*72 

301 

3291 

0291 

3733 

*94* 

*01 

*00  1 

A«2* 

*62* 

6*73 

316 

*47* 

13*8 

1*2* 

6967 

112 

oon 

69*6 

oool 

6*7* 

000 

1*23 

0030 

1366 

6570 

013 

49*6 

3394 

6366 

6*73 

300 

1*2* 

0000 

1973 

6571 

090 

0000 

00  00 

0000 

6*76 

300 

1*23 

0030 

160* 

6372 

112 

coll 

*9*9 

3  0  C  1 

6*77 

000 

1*26 

0030 

1*22 

6373 

014 

6374 

7331 

7  *  l  0 

6900 

000 

1*27 

0030 

1*2* 

63?4 

032 

6*2* 

00*1 

6  6  ?  3 

*301 

000 

1*30 

ooon 

35*2 

6373 

•3 1  3 

**** 

J*7l 

6*4* 

*302 

000 

1*31 

0000 

3970 

69  ?6 

0  54 

66** 

2023 

636* 

*303 

032 

3000 

0030 

0  000 

6977 

000 

0000 

0  000 

0900 

*30* 

3  0  0 

1*32 

0000 

372* 

*4  0  0 

112 

3  012 

*3*4 

9001 

6309 

112 

0003 

633  4 

0001 

6601 

300 

*4*3 

0000 

656* 

6906 

Q5« 

0000 

6140 

1420 

6602 

009 

0000 

0008 

0900 

6907 

432 

0300 

0030 

*912 

6603 

009 

32*9 

0000 

6567 

6910 

190 

3000 

3033 

*900 

6604 

300 

32*6 

0003 

6570 

6911 

112 

211* 

*  3  i  0 

00C1 

*603 

01* 

*60* 

63*2 

*  6  c  2 

6912 

000 

0000 

30  30 

0  SCO 

6606 

000 

32*7 

0000 

63*7 

6913 

000 

0000 

3  0  3  0 

0  0  C  9 

6607 

336 

32’0 

31*9 

*370 

691* 

09* 

3000 

131* 

0000 

6610 

300 

0000 

OOCO 

0000 

6919 

900 

3000 

3000 

0003 

6611 

900 

0000 

0000 

0000 

6916 

000 

0030 

3000 

0000 

6612 

300 

0000 

0000 

0000 

6317 

000 

0000 

0000 

0  OC  9 

6613 

300 

0000 

0000 

0000 

6970 

900 

0000 

3000 

0300 

661* 

000 

0000 

0000 

oooc 

6321 

000 

QQ30 

•3000 

0000 

6613 

000 

0900 

1000 

ooooN*# 

69*2 

390 

0000 

0000 

0000 

6616 

300 

0300 

0090 

ooooit(n) 

6923 

300 

OCJf 

oooc 

tooo 

**1  ? 

390 

coco 

0  0.0  0 

330Cft(4j 
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*62* 

000 

0000 

oooo 

oooofUp.) 

6?02 

713 

0001 

0090 

0002 

i«2t 

000 

0000 

0090 

ooooalH) 

67)] 

921 

0002 

0009 

6002 

6*22 

000 

0000 

0360 

oooo  It (<(.*}[ 

670* 

000 

OOOO 

oooo 

00006 

6*23 

000 

0000 

0600 

oooo 

6703 

ooo 

OOOO 

oooo 

OOOO 

6*2* 

000 

0000 

0900 

oooo 

67o* 

000 

OOOO 

oooo 

OOOO 

662} 

000 

0000 

oooo 

oooo 

.  * 70 7 

002 

2*73 

2602 

2*01 

*626 

000 

0000 

oooo 

9000 

*710 

03* 

OOOO 

6*02 

OOOO 

6627 

000 

0000 

oooo 

oooo 

*711 

003 

1*23 

1*37 

1*37 

6*30 

000 

0000 

oooo 

oooo 

6712 

036 

OOOO 

3*77 

OOOO 

*63i 

000 

oooo 

oooo 

oooo 

6713 

3*3 

0*13 

263* 

20  7  3 

6*32 

000 

0000 

oooo 

oooo 

671* 

100 

OOOO 

OOOO 

OOOO 

•  633 

000 

0000 

oooo 

oooo 

*713 

073 

*132 

*3*2 

0730 

663* 

000 

0000 

oooo 

oooo 

*716 

130 

3360 

3661 

0 7o  7 

f!nn 

6*33 

000 

0000 

oooo 

oooo 

6717 

100 

OOOO 

OOOO 

OOOO 

6636 

000 

0000 

oooo 

oooo 

*720 

077 

*7*o 

73*2 

1*1*1 

for  draw- 

6637 

000 

0000 

oooo 

oooo 

6721 

620 

*737 

*003 

7336 

66*0 

000 

oooo 

oooo 

oooo 

6722 

100 

OOOO 

OOOO 

OOOO 

6**1 

000 

0000 

OOOO- 

oooo 

6723 

100 

*237 

*003 

7336 

6**2 

000 

oooo 

oooo 

oooo 

*72* 

771 

3636 

037* 

702* 

L 

6**3 

000 

oooo 

oooo 

oooo 

6723 

100 

OOOO 

OOOO 

oQqo 

66** 

000 

oooo 

oooo 

oooo 

672* 

0  77 

7*7* 

1373 

*030 

i 

66*3 

000 

oooo 

oooo 

oooo 

6727 

367 

7  30  3 

337* 

*7*2 

66*6 

003 

67*3 

*730 

6631 

*730 

100 

OOOO 

OOOO 

OOOO 

lR* 

66*7 

016 

6630 

7301 

7610 

6731 

loo 

7303 

337* 

*7*2, 

6*30 

000 

6631 

000* 

*631 

*732 

377 

7313 

0019 

73ft*1 

»D 

6*31 

002 

oooo 

6*31 

6631 

*733 

100 

OOOO 

OOOO 

OOOO 

6632 

000 

oooo 

oooo 

oooo 

*73* 

loo 

7313 

0019 

7366, 

*633 

000 

oooo 

oooo 

oooo 

*733 

0*0 

6133 

0313 

2322 

*0 

6*3* 

000 

oooo 

oooo 

oooo 

673* 

109 

OOOO 

OOOO 

oooo 

*633 

000 

oooo 

oooo 

oooo 

*737 

100 

*133 

9913 

23  2  2, 

6*36 

000 

oooo 

oooo 

oooo 

67*o 

012 

*71* 

0669 

6*73 

>D 

6*37 

000 

oooo 

oooo 

oooo 

67*1 

100 

0000 

OOOO 

oooo 

*6 

66*0 

000 

oooo 

oooo 

oooo 

67*2 

100 

*71* 

0*69 

6*73, 

6661 

000 

oooo 

oooo 

oooo 

67*3 

*33 

*1*9 

6907 

0167 

n 

6*62 

000 

oooo 

oooo 

oooo 

*7** 

100 

OOOO 

OOOO 

OOOO 

6663 

000 

oooo 

oooo 

oooo 

67*3 

100 

*169 

*307 

01*7, 

Is 

666* 

000 

oooo 

oooo 

oooo 

*7*6 

07* 

3332 

772* 

1*09 

6663 

000 

oooo 

oooo 

oooo 

67*7 

loo 

OOOO 

OOOO 

OOOO 

,R« 

666* 

000 

oooo 

oooo 

oooo 

*730 

0  »7 

7323 

7*3* 

30  0* 

*6*7 

000 

oooo 

0  0-0  0 

9000 

*731 

032 

oooo 

*733 

OOOO 

6*70 

090 

oooo 

oooo 

oooo 

6732 

01* 

*733 

6*7* 

6  7g  * 

6*71 

000 

oooo 

oooo 

9000 

*753 

000 

*737 

eooe 

*632 

6672 

000 

oooo 

oooo 

oooo 

*73* 

092 

OOOO 

*7*0 

OOOO 

6673 

000 

oooo 

oooo 

oooo 

*733 

01* 

*796 

6*76 

*70* 

6*7* 

000 

oooo 

oooo 

9000 

6736 

000 

*7*2 

oooo 

*633 

6673 

000 

oooo 

oooo 

oooo 

*737 

092 

OOOO 

*7*3 

oooo 

*67* 

39* 

010* 

oooo 

0002  6 

*7*0 

01* 

*7*1 

*6  7* 

*7o* 

■ 

6*77 

n5 

OOOO 

9002 

000  3  Transducer  6  76  1 

032 

OOOO 

*7*4 

0060 

6700 

33* 

00*9 

oooo 

0  00 2  for  uni- 

*7*2 

016 

66*6 

6*7* 

*7j  * 

»701 

?g? 

0002 

oooo 

0  0  0  0 

wywg  formly 

67*3 

000 

0009 

0009 

OOQO 

i 

, 
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(0) 


•  100 

•  ♦03 

100000000*1 

0  172 

♦  ♦00 

0101 

•  ♦0  J 

toioooooo 

0  173 

♦*oo 

0102 

♦  ♦03 

102000000 

>Ho 

0  174 

♦*o< 

0103 

♦  ♦03 

160000000 

0175 

♦♦Of 

0104 

♦  ♦0  3 

200000000 

0  176 

♦♦Of 

0105 

•  ♦03 

300000000 

0  177 

♦♦0# 

0106 

♦  ♦03 

500000000 

0200 

•  f  00 

0107 

♦  ♦  G  4 

isoooooo  cn 

020  1 

♦♦Of 

ono 

♦  *0  4 

700000000 

0202 

♦♦II 

0  111 

♦  ♦06 

700000000 

A/  ,  ,  oa°3 

♦♦•1 

0112 

♦  ♦04 

70  1  QQQQOQ 

^Yoo[m/s]  0  20* 

♦♦Oo 

0113 

♦  *Q  4 

770600000 

0205 

♦  400 

0  114 

•  •o* 

77200300C 

0206 

♦  ♦•0 

0115 

+  ♦0* 

761500000 

02Q7 

♦  4#0 

0116 

•  ♦04 

1 120379001 

0210 

♦♦Of 

0117 

112037900 

0211 

♦*Oo 

0120 

♦  ♦•4 

1 12037900 

.7  [  m/  s  ] 

0212 

♦♦Oo 

C121 

♦  *0  3 

954700000 

0213 

♦♦Oo 

C  122 

♦  ♦0  4 

106300000 

0214 

♦  ♦0  1 

C  123 

117200000 

0215 

♦-•  l 

C  124 

♦  •04 

135100000, 

0216 

♦  *0  1 

0  125 

♦  ♦04 

lOOQOOOOO 

0217 

♦  *00 

9126 

•  •04 

100003000 

0220 

♦  ♦Oo 

0  127 

•  •04 

100000000 

T«fK] 

0221 

*  *0  0 

C  1 3o 

♦  ♦0* 

120700000 

0222 

♦  *9, 

0131 

•  •O  4 

1404000C0 

0223 

•♦It 

0132 

♦  ♦04 

1423030C0 

0224 

••00 

0133 

♦  ♦04 

1576030001 

0225 

•♦00 

013* 

♦  ♦00 

6004030017 

022« 

♦♦00 

0;35 

♦  ♦.3  0 

575903000 

0227 

•♦00 

0136 

♦  ♦30 

540000003 

0230 

-♦0  1 

3137 

♦  ♦30 

500000003 

oLd 

0251 

-♦01 

3  140 

•  ♦00 

445000000 

0233 

-•01 

01*1 

•  ♦00 

362000000 

C  233 

-♦01 

01*2 

•  ♦00 

24000300a 

023* 

•♦00 

0143 

•  •0  0 

900003000 

0235 

♦  ♦ot 

0144 

♦  ♦00 

890003000 

023  c 

•♦0  1 

3145 

♦  •00 

870003000 

I  „ 

0237 

♦♦0  1 

51*6 

♦  ♦00 

040009000 

02*0 

♦  ♦0  1 

31*7 

♦  •90 

795000000 

02*  1 

••0  1 

0V50 

♦  ♦00 

730003000 

02*2 

•  ♦01 

0151 

♦  ♦00 

630000000 

02*3 

•♦0  1 

0152 

♦  ♦01 

1*0000600 

02*4 

•  ♦ci 

0153 

♦  •01 

146000000 

02*5 

•-03 

0154 

♦  ♦01 

139000000 

'aim 

02*9 

♦♦01 

0155 

♦•01 

124000000 

02*7 

•  ♦01 

0  156 

♦  ♦01 

110000900 

0230 

•  ♦03 

0157 

♦♦00 

9*0000000 

3251 

•  •01 

0160 

♦*0o 

780000001; 

0252 

♦  ♦01 

0161 

♦♦01 

12000000c 

3253 

•  •00 

0162 

4*0, 

117500000 

9254 

♦  •00 

0163 

♦♦•l 

1  1  1000000 

Jju 

82  55 

♦•06 

0164 

♦  ♦Of 

ooooeoooo 

C2$6 

*•00 

0165 

♦  ♦»* 

OOOOOOOOQ 

0257 

•♦02 

0166 

♦*0q 

505000000 

0260 

♦♦00 

0167 

♦♦do 

300000000 

02*  t 

•♦00 

0170 

♦  ♦«, 

ooooeoootf 

0262 

♦♦00 

0171 

♦  40* 

875000000 

Ln 

0263 

♦  ♦00 

5*00000001 


JLTi 


775000001 
100000000 
121000000 
utoooooo 
ooooooood 

145000000  I 
200000000 
♦*5ooooooUt> 

575000006  *■ 

645000000 

700000000 

ooooooood 
2S8000000 
500000000  ,  D 

707100I10C’  jULva 
•66000000 
065000000 
1000C0000 
400000000 
821000000 
129*00000 
222*70000 
542700000  W 
6040  ]  0 (i  0  C  I 

7715s*too 

071480000 
919760000 
9*5810000 
901220001 
45500000 j 
300000000 
200000080 
100000000 
000080000 
lOOOOOOOOV^ 
200000000  «« 
300000000 
400000000 
500000900 
50000009S 
t  ooeooooo 
i  ooeooooo 
ioooooooo€, 

*55000000 

152100900 

40ooooeoo7wl*Kj 

325000000 
100000001 
000  000900  Nip 
ooeoooooo Nip 
lOCOOOOOOCi 
OOCOOOOOO 
100000900  r 

ooeoooooi 
ooooooooq 
ooooooooq  . 
ooooooooiflL*  [ de 
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C2*4 

*♦00 

030100000 

1  r  i 

.1342 

4+00 

088000000  Xh. 

32*5 

•*0# 

0U0O0O5Q0 

tc«[  deg] 

03*3 

4*80 

000000903  Xa 

C2«6 

•♦00 

000003300 

13*4- 

+♦3  3 

ooooooooo 2^3 

0267 

••00 

OU04033QO 

034b 

4+15 

l 00 000800  Lt 

0270 

•♦00 

000000300 

1346 

441+ 

♦Oooooooo  Li 

32  M 

*♦00 

0OOO033QO 

3347 

4*80 

ooooooooo 

0272 

*♦00 

000800300 

3360 

4*80 

OOUOOOOOO 

3273 

♦♦00 

OU00003QO 

03b  1 

♦*00 

000003000 

0274 

4*0  0 

00300030 J 

3  3b2 

••80 

ooooooooo 

3275 

4*00 

ooooonooo 

0363 

*•88 

ooooooooo 

027* 

443  0 

ooooooooo 

0364 

••00 

OOOOC3QOO 

0277 

4*40 

ooooooooe 

0355 

♦  *eo 

000033000 

4300 

3*00 

OQOOOUOOU 

C  356 

♦•00 

OOOOOOOOO 

0  Jo  1 

4*10 

ooooooooo 

^[deg] 

0357 

*•80 

OOOOOOOOO 

4302 

♦  ♦30 

ooooooooo 

0  36  C 

**00 

800000000 

43C3 

♦  •00 

ooooooooo 

336  1 

*•00 

OOUOOOOOO 

4304 

♦  ♦30 

OOOOOOOOO 

0362 

♦  •ou 

ooooooooo 

13Cb 

♦  ♦00 

ooooooooo 

3363 

♦•CO 

000003000 

4306 

4*00 

OOOOOQOOU 

8364 

♦  •oo 

ooooooooo 

4307 

♦  •00 

ooooooooo 

0365 

**00 

000009000 

33  10 

+  •00 

OOOOOOOOO1 

0366 

♦•00 

ooooooooo 

43  1  1 

♦■•oo 

ooooooooo 

•H 

0367 

♦•00 

ooooooooo 

4312 

4*0  0 

OOOOOOOOU 

0  J70 

**00 

QOQOCOOOO 

43  13 

4*00 

UOOOOOUOO1 

1 

037  1 

**oo 

OOOOOOOOO 

1314 

4*00 

OOOOUUOUU 

0372 

♦*00 

OOOOOOOOO 

131b 

+  •00 

ocoooOooc 

3473 

•*00 

OOUOOOOOO 

1316 

+♦00 

ooooooooo  It* 

0374 

*•00 

ooooooooo 

1317 

4*0  0 

OOOOOOOOO 

Xft 

5375 

•*00 

OOUOOOOOO 

1420 

♦•00 

ooooooooo  X-rz 

3376 

♦*00 

OOOOOOOOO 

132  1 

4*0  0 

OOOOOOOOO  JCfi 

0377 

♦  *co 

OOOOOOOOO 

1322 

4*00 

ooooooooo  Xrs 

0«QQ 

♦♦CO 

OOOOOOOOO 

1323 

4*0  0 

OOOOOOOOO  <It3 

3*0  1 

*♦00 

OOOOCOQOO 

■>324 

4*0  0 

OOOOOOOOO 

3402 

♦*00 

000003000 

•32b 

4+00 

5*« 1805804WT 

3403 

♦•00 

080000000 

C  326 

4-09 

loooooooa 

0*04 

♦•00 

800080000 

0327 

•♦00 

ooooooooo 

043  5 

•  •oo 

ooooooooo 

0330 

+  ♦00 

Q0030OU00 

0  4  3  6 

•  •oo 

000005000 

0331 

-♦♦0  1 

loooooooa 

3407 

•  •oo 

OOOOOOOOO 

0332 

4*10 

OOOOOOOOU 

a  a 

0*  10 

•  •oo 

000003000 

0333 

4*00 

ooooooooo-N+p 

C*  1  1 

•  •ol 

350000030 

0334 

4*0  0 

C‘ 0000000  O-i© 

3*12 

••81 

473238898'3% 

033b 

4*0  0 

ooooooooo 

3*13 

♦  *  8  l 

14142135*^2" 

0336 

4*00 

ooooooooo 

it1 

3*1* 

••01 

177245384^y 

0337 

*♦0  0 

ocooooooo 

xf 

14  15 

♦  •oo 

80*22  7  0  0  0 

CJ40 

4*00 

UOOOQQOOO 

52, 

3*16 

♦•8  1 

*  2  83  1 8  53  1  2.% 

«  J*t 

♦+•0 

UO  00  OOOOOjHU 

3  4  1  7 

♦  ••1 

31+1592*5  % 
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BLOCK  II  (M) 


2001 

013 

3930 

64  J  5 

0000 

2002 

03* 

0000 

306* 

0000 

2313 

03* 

1420 

3632 

2402 

2033 

030 

0000 

2640 

0000 

2030 

003 

1*17 

1*3? 

1*33 

2037 

000 

0000 

0000 

0000 

2040 

004 

2*32 

1476 

1*23 

2041 

010 

2042 

6046 

2637 

2*42 

002 

1*33 

1*29 

1*22 

2043 

002 

1*3* 

1*2? 

1**4 

2001 

002 

1*33 

1439 

1*22 

2002 

004 

1*3* 

1*22 

1*24 

2003 

002 

1*24 

0213 

0000 

20*4 

03* 

0230 

6796 

1*20 

2**9 

002 

0227 

1424 

0000 

3049 

001 

77*1 

2673 

2673 

3041 

013 

3990 

7722 

3990 

*373 

002 

0313 

0000 

OOOO 

0574 

*3* 

0000 

0402 

0000 

«37» 

092 

0313 

2001 

142* 

0*40 

004 

2493 

1470 

1*27 

0*47 

002 

7701 

1*33 

1*33 

**5e 

002 

1  4  3  * 

1*2(1 

1*36 

C710 

•56 

0008 

6373 

0008 

*732 

002 

2*02 

2073 

oooo 

•  733 

030 

7701 

*707 

1423 

474* 

005 

142* 

031* 

1*23 

*747 

001 

1423 

033* 

1429 

•  730 

044 

1*23 

0000 

1*29 

*73l 

001 

0313 

1423 

1423 

*?32 

002 

7761 

1*2* 

1426 

0793 

.00  9 

7762 

1*20 

1426 

^0794 

004 

1423 

1*26 

1*20 

3314 

002 

0940 

0000 

oaooNumber 

« 

0332 

001 

04C0 

0000 

oooo  Number 
1,0 

Note:  introduce 

data  into 

2003  OS*  0000  3*0*  0000 

2004  000  0000  0000  0000 


2***  00*  14*3  14**  14** 
2*45  00*  1423  14*4  1424 
2*40  00 *  1*23  1439  l4|S 
2049  ooi  i4**  14*9  1*33 
2* So  001  1*24  142?  1434 
2091  03*  *1*?  2**1 


2***  03*  00*0  *?4*  0000 
200?  010  2*7*  7501  7*10 
2070  000  0012  0190  o2l$ 
2*71  0?»  14*4  0230  1420 
2*72  030  0000  «?30  00*0 

3042  03*  0000  2001  0000 


037*  03*  0000  8572  00*0 
*377  004  1422  0313  1423 
*400  09*  0000  0409  0000 

•*3l  002  143*  1433  14*7 
0*52  003  143*  143*  1454 

0033  Q30  0000  2*30  oO|0 

*734  010  2*31  0*47  **3? 


*733  044  1420  OQOS  14*0 

•73*  044  14*3  0000  14*3 

073?  000  1**0  0000  0041 

07*0  010  07*1  *439  01*7 

•7*1  044  1434  0000  145* 

*7*2  QOS  770*  0041  0041 

•  70S  03*  141*  07ii  *1  o  7 

0313  ooi  o9*o  oo*o  on o o  Number 

1.5 

0334  xo*  04*7  2400  oooo  Number 

•n,75 

cells  0313,1420,5222-5242 
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BLOCK  III  01) 


2313 

036 

0000 

3*32 

0300 

6401 

05* 

2039 

2640 

6371 

64*2 

092 

OOOO 

OOOO- 

OOOO 

2*01 

030 

3000 

0011 

OOOO 

6403 

0  1  6 

6404 

7301 

761» 

2602 

l  7  J 

1476 

6773 

6173 

6404 

000 

0010 

0130 

67*4* 

• 

• 

*403 

1  *3 

147* 

*773 

*373 

2633 

0  7* 

0000 

2640 

6336 

6406 

013 

640  9 

2001 

640  3 

6407 

ll  2 

0002 

6-0  3 

0001 

2641 

016 

66c2 

6646 

2637 

6410 

03* 

2602 

2642 

*403 

2644 

003 

1422 

6373 

1622 

6707 

003 

2633 

334> 

2693 

2643 

003 

1624 

6374 

162  4 

6710 

03* 

OOOO 

6373 

0900 

2646 

003 

1435 

6  3  73 

1433 

*711 

903 

6373 

1437 

1*37 

6373 

092 

7032 

0300 

*732 

003 

77*2 

35*7 

2692 

6374 

000 

7354 

3000 

7  54ijrn 

*733 

034 

7736 

77*2 

2633 

6373 

016 

6376 

7  301 

7613 

*734 

096 

OOOO 

*7*7 

OOOO 

637* 

352 

*764 

0042 

7031 

6377 

J92 

7032 

3  5  3  0 

673* 

0*4 

6373 

oooo 

6373 

6400 

030 

7304 

0030 

7541JP'1 

BLOCK  IV 

on 

2217 

036 

0000 

3037 

oooo 

362g 

QC3 

OOOO 

3000 

0503 

(3*21 

01* 

3622 

7301 

7610 

3064 

096 

0000 

6230 

oooo 

l>622 

0  >2 

4463 

00*1 

*312 

3623 

112 

oooo 

3140 

0001 

3122 

000 

0000 

0000 

oooo 

3624 

096 

3137 

3143 

3623 

' 

>*23 

032 

OOOO 

oooo 

OOOO 

3130 

304 

4001 

1633 

4001 

3626 

034 

Oil* 

6*1* 

0001 

>131 

306 

4232 

1634 

4232 

3627 

013 

3*23 

OOQl 

3623 

>132 

112 

0230 

3130 

0001 

3630 

033 

3423 

772* 

3623 

3133 

092 

0000 

oooo 

oooo 

3631 

036 

OOOO 

3140 

OOOO 

>134 

303 

4001 

0331 

4301 

3133 

112 

0210 

3134 

0021 

53*4 

0  3* 

OOOO 

3706 

0300 

>136 

0>6 

0000 

3623 

0300 

>137 

H2 

0000 

3140 

0001 

3713 

034 

Oil* 

772* 

1437 

3140 

600 

4001 

0000 

4463 

371# 

096 

1437 

6022 

1436 

>141 

400 

4022 

0000 

4464 

3717 

033 

>020 

143* 

9920 

>142 

600 

4o*3 

0000 

4663 

3720 

053 

1*37 

1*36 

1434 

3143 

400 

4232 

0000 

4466 

3721 

03* 

OOOO 

616* 

0000 

>144 

436 

4233 

3377 

44*7 

.  *043 

'*36 

3222 

*01  7 

9029 

3377 

400 

4274 

0000 

447Q 

3600 

600 

4064 

0000 

4471 

6234 

034 

0090 

1440 

0901 

3601 

400 

6103 

0000 

4472 

6233 

033 

0001 

7722 

0001 

>602 

600 

4126 

3000 

*473 

6236 

072 

oooo 

0001 

OOOO 

3603 

6Q  0 

4313 

oooo 

4476 

*237 

01* 

31>0 

3645 

3633 

3*06 

600 

4336 

0000 

4*73 

3603 

600 

4337 

oooo 

4476 

6273 

01* 

3*67 

3046 

0302 

3*06 

400 

6147 

oooo 

4*77 

1*07 

600 

4170 

oooo 

*300 

4300 

032 

OOOO 

6317 

0017 

'610 

400 

4211 

0030 

4301 

36U 

400 

4400 

oooo 

4302 

6306 

07* 

oooo 

6312 

OOOO 

3612 

600 

4421 

3000 

4  lQ  3 

3613 

4oo 

4442 

3000 

430* 

6307 

01* 

0234 

247* 

2713 

3614 

432 

0000 

0030 

3620 

»*13 

701 

4463 

6471 

6303 

6313 

01* 

3*43 

2*76 

2713 

3616 

701 

43:3 

4477 

k  Jg  3 

JM  ’ 

1 1  2 

c : :  5 

J  6  l  3 

0?01 

4317 

ol« 

3073 

1046 

0302 
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BLOCK  V  (B) 


9 

4914 

0  9* 

0000 

4019 

411 

0 

9 

4019 

0** 

491* 

4000 

002 

9 

4194 

000 

002* 

0000 

431 

2 

1 

4l»3 

39* 

002* 

42*7 

431 

9 

4299 

000 

0000 

0000 

002 

*290 

000 

0000 

oeoo 

002 

9 

9 

4291 

0*« 

4*14 

*117 

410 

9 

4497 

Q?« 

9011 

490* 

*39 

9 

4992 

0** 

4309 

430* 

*99 

2 

49g* 

0000 

3344 

000 

1 

4999 

on 

9114 

4911 

111 

9 

4910 

09* 

9012 

4090 

1*3 

9 

*311 

0** 

*012 

3230 

403 

9 

4912 

000 

*01* 

0000 

411 

1 

4313 

0*« 

4014 

4**1 

412 

0 

4914 

001 

0049 

0044 

004 

2 

4919 

000 

0000 

0000 

004 

9 

491* 

000 

0000 

0000 

000 

*919 

0*0* 

*000 

000 

1 

2 

*929 

i'S 

eooo 

4337 

000 

1 

0 

4949 

092 

9901 

0041 

9*7 

i 

1*991 

01* 

4392 

7901 

7*1 

« 

3 

1*992 

092 

3271 

0041 

327 

« 

*9*0 

00* 

0000 

0000 

000 

9 

49*1 

•  00 

0000 

*000 

000 

9 

4*9* 

000 

000* 

0000 

327 

0 

9 

4*97 

*00 

0000 

0000 

327 

9 

9*12 

*01 

0043 

0343 

009 

0 

*10* 

*09 

0249 

0243 

000 

*109 

909 

*244 

0244 

000 

9 

*10* 

00* 

77*1 

0001 

32* 

9 

*197 

004 

77*1 

000* 

32* 

2 

*110 

009 

0*49 

0244 

000 

9 

’111 

009 

0*1* 

0001 

32* 

2 

*112 

004 

32*2 

0024 

32* 

9 

*119 

09* 

*11* 

311* 

*11 

9 

*114 

004 

3317 

3921 

142 

9 

*119 

09* 

000* 

3237 

000 

*11* 

019 

»72* 

3344 

*00 

3 

*119 

03» 

00*0 

4912 

000 

4 

*120 

019 

77*4 

3344 

000 

1 

*121 

09* 

77*1 

491* 

32* 

Note:  introduce  data  into  cells  0243.0244.3344 
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BLOCK  VI  (K) 


09*  0000  0603  0000 

1992  001  0003  0303  0090 
1933  090  0000  3033  0000 
193*  001  00*3  03*3  0099 
1939  q9*  0000  322*  0000 

292T  §9*  0000  2720  2924 

2717  o9o  0000  3037  0003 

2720  094  003*  1**0  0001 

2721  ols  0001  0000  0003 

2722  03*  0000  2930  0000 

2723  000  0000  0000  0000 
272*  000  0000  0000  0000 

2729  000  294*  QOOO  0301 

2726  000  2«00  0000  0002 

2727  ooo  3346  qqoo  294* 

2730  016  2731  3l?l  2*09 

2731  ooo  QOOl  0000  2944 

2732  000  0 0 0 2  0000  2600 

2733  009  1*17  0o91  9ooa 

2734  ol9  9000  0090  0000 
2733  j3*  7761.  3091  4777 

2736  016  2737  7501  7*10 

2737  q92  1490  0033  0069 
27*0  0 92  1420  0000  0003 
27*1  0Q4  7761  9qqo  9000 

2742  Q09  1*74  9 o o o  900* 

2745  43*  ooQO  0000  279o 
2744  309  1420  9000  9001 
2749  7o9  9001  1*20  9009 

2746  3o9  1420  300*  9013 

2747  112  0002  274*  OOOl 

2790  OOO  0000  0000  0003 

2791  ooo  9oal  1421  3010 
2732  009  9ao2  1422  9311 
2 »33 .00  9  9QQ1  1422  9012 
2794  q9«  00>0  1440  oqo! 
2799  031  0001  7722  0001 
2796  072  0000  OOOl  0000 
279?  002  1417  0000  0009 
27*o  03*  0000  2766  0000 
27*1  301  4777  *001  4901 
27*2  013  27*1  27?J  27*1 

2743  019  2741  27?*  (TOOO 
27#*  07*  3000  27*1  OOOO 
27*3  096  277*  333"  27*1 
27**  3oi  4777  4021  4321 
27*7  OH  2»6*  17  7  3  27*6 

*04*  09*  38Q0  301*  0900 

3106  OQO  0000  3000  OOCC 

3107  09*  0000  3111  3990 


1336  003  0062  0002  0*02 

1337  ooft  0091  32*3  0601 

1940  0 ®*  QOOO  3026  0600 


2770  013  2706  27»j  ,*aa 

2771  070  0600  270*  0600 

2772  036  2777  3000  27*6 

2773  OOO  3001  00*0  0640 
277*  3q1  3010  4741  4741 
2773  301  3016  4701  *7*1 
2770  301  *77?  4001  4001 
277?  301  4777  4021  *621 
3006  066  6666  0000  0600 
3001  060  *329  0660  0846 
3662  OOO  *920  0090  004? 
3003  OOO  *927  0000  0090 
300*  OOO  0000  0000  OQOQ 
3qq3  09*  0000  7430  0000 
3006  OOO  *523  COOO  3011 
3*0 7  OOO  692*  0000  3012 

3010  036  6927  627}  J0}3 

30 1 1  OOO  0000  0000  0000 

3012  ooo  0000  0000  OOoo 
3al3  OOO  0000  0003  C 0 0 9 
>01*  302  3322  *922  3922 
3019  OOO  3011  0000  *929 
301*  OOO  3012  0000  692* 
3317  096  3013  3Q67  *927 
3:20  OOO  03*0  00*0  0020 
9q21  000  0000  0000  0600 
3022  OOO  0020  0000  0000 
3029  OOO  OOOO  0020  0000 
302*  OOO  OOOO  OOOO  0020 

3029  OOO  00*0  OOOO  0060 
3026  009  0002  3264  e002 
>027  SOI  OOOl  0002  0*01 

3030  002  77*1  0001  0  9  0  9 
3331  03*  0900  3033  OOOO 
3932  0  9*  OOOO  3039  OOOO 
3033  031  77*}  4404  6404 
303*  03*  OOQO  1)94  OOOO 
3339  091  77*i  *607  64*7 
33  36  096  OOOO  1*3 J  oOOO 
3037  69*  QOOO  3006  OOOO 


3110  034  OOOO  3lii  OOOO 


>122  OOO  OOOO  3930  33?' 
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*12*  003  330*  773’  0300 
*130  07*  3321  3132  1*21 
*131  03*  331*  2101  1*20 

*170  03«  0000  271?  0000 
31*1  000  0370  0000  00(70 

*22*  002  00**  3233  0001 
*223  002  0030  323*  0302 
*22*  003  0001  0001  0301 
*22*  03*  oooo  1336  0000 
*230  001  03*1  03*2  *233 
*2*1  001  0***  0***  3236 
*232  00*  3233  7762  3233 

*000  03*  017  *12*  002* 
*001  013  33**  772*  0003 
*002  0*6  *00*  *0 o 3  1632 
*003  013  33**  7722  0000 
*00*  0*6  *00*  *003  1*32 
*00*  03*  0003  *030  Q300 
*00*  0*6  0000  1332  0000 

*17*  000  0023  oooo  *3l« 

*267  000  OOOO  OOOO  0023 
*270  000  OOOO  OOOO  0026 

**74  003  330*  773’  0300 
**77  o’*  3*11  *30*  *333 
*300  000  *303  OOOO  *13<J 
**01  000  *30*  OOOO  311* 
**02  03*  *3q  3  *3o«  *333 
**0S  000  77*1  OOOO  00*1 
**0 *  00*  3310  77*2  1*17 
*3o3  003  3310  33iC  1*32 
*9o*  013  33**  OOOO  OOOO 

*320  03*  OOOO  *33’  OOOO 

*3*7  0  32  3301  00*1  3671 

*3*0  000  0  0  Q  3  OOOO  OOOO 

**3*  000  OOOO  OOOO  (ll( 

*7i0  03*  331*  5013  3317 

3012  001  0 C  *3  03*3  0330 
“•13  013  331*  OOOO  0300 
901*  07*  OOOO  *7 1 1  oooo 

910*  003  02*3  02*3  0001 
9103  003  02**  02**  0002 
9106  QO*  77*1  QQol  32*3 
”3 0 ’  00*  ’7*1  q 0 0 2  32** 
9110  003  02*3  02**  0001 
*111  003  0*1’  0001  32*2 
9112  00*  32*2  002*  32*2 


3132  003  3317  3513  1*20 

3133  0  3*  QOOO  2101  6000 


3l’2  001  02*3  00’9  0070 
Sl’3  0 3*  oooo  2331  OOOO 

3233  00*  373*  77*2  3236 
323*  03*  OOOO  310*  OOQO 
3233  000  QOOO  OOOO  OOoC 
323*  000  QOOO  OOOO  OOQO 
3237  00*  J51?  6*07  i*2o 
32*0  003  l*2o  3262  1*20 
32*1  03*  0000  31  1*  0 0 0 0 

*00’  09*  oooo  1334  0### 

*010  072  OOOO  6*1*  qooo 
*011  *3*  o?13  *01*  3303 
*012  056  *313  *0®6  *10* 
*013  036  OOOO  *17*  *071 
*01*  0 3*  3000  *013  *111 
*013  q36  *316  *000  0023 

*173  036  0026  *267  *3l7 

*27i  go*  4^4  *2  3  7  4l64 


*®0  7  07*  Ml*  *311  3113 
*310  03*  5012  *030  1632 
*311  036  *012  3230  *032 
*312  000  *013  OOOO  *11* 
*313' 096  *91*  *001  *122 
*31*  ooi  00*3  00**  00*6 
*313  000  OOQO  OOOO  00*6 
*316  000  OOQO  OOOO  OOQO 
*317  900  0000  OOOO  OOOO 


*361  QOO  OOOO  0033  OOQO 
*«17  990  9909  OOOO  ««07 


3013  QOO  oeoo  OOOO  *7 j 1 
9016  036  30  17  *7u  *3*3 
3017  100  OOOO  OOOO  1*1’ 

9U3  03*  3113  311*  3ii  3 
911*  00*  3917  3321  1*23 
3113  036  OOOO  3237  OOQO 
3116  019  7722  93**  oflo9 
»117  036  OOOO  *312  gooo 
9120  013  772*  33**  0000 
3121  «30  77*1  *312  3262 


DOC  «  83067305 


page  g&S 


Page  140 


>330 

•  34 

OOCO 

3*04 

OOOO 

9331 

013 

33*4 

0000 

OOOO 

>332 

•  >• 

0000 

2109 

OOOO 

9333 

003 

1*20 

3242 

1*20 

>33* 

030 

6*3* 

0000 

1*21 

9335 

•  04 

1*20 

1421 

1*22 

533* 

003 

3310 

i*a? 

1*23 

5337 

004 

1*22 

3910 

1424 

33*0 

003 

7762 

*749 

1*2* 

93*1 

009 

1*2* 

142* 

1*23 

3342 

004 

1*2* 

3910 

1*26 

5345 

932 

0000 

233* 

OOOO 

>>44 

>04 

*C0  1 

1421 

4001 

>345 

303 

*0*1 

1  *23 

4041 

534* 

903 

*2*1 

l*2» 

4241, 

3347 

112 

0  C  3  * 

211* 

0001 

9330 

0>2 

COOO 

2299 

OOOO 

>351 

303 

*131 

1422 

4101 

3332 

503 

*301 

1*24 

43«1 

5333 

303 

4**1 

1424 

4**1 

9334 

303 

*4ol 

1*24 

*601 

3333 

112 

013’ 

2121 

0001 

333* 

*32 

OOCO 

0000 

2132 

3337 

*01 

*c*l 

*041 

4*41 

5340 

112 

001* 

212’ 

0901 

3361 

013 

212* 

3020 

2i  a* 

>342 

ooc 

COOO 

0000 

OOOO 

9343 

112 

C  1  3  3 

212* 

0301 

>36* 

OOO 

0  0  '.0 

0000 

0  0  0  c 

3363 

0>2 

0000 

0000 

OOOC 

>364 

*13 

*7*1 

OOOO 

0300 

3347 

*36 

*  *,*  1 

21  ’2 

OOO* 

3370 

70* 

*  761 

4**1 

4*6; 

9371 

70  * 

3031 

4  ’4 1 

3001 

9372 

70  * 

3c  21 

4*41 

3021 

93*3 

70* 

3041 

*7*1 

33*1 

3374 

0  2  0 

0000 

OOOO 

0009 

9375 

OOO 

0003 

OOOO 

0300 

337* 

033 

*  *61 

4761 

0001 

3377 

633 

3  0  C  l 

9001 

0002 

3*00 

403 

9321 

3021 

0003' 

>401 

031 

0  0  Q  1 

0002 

COO* 

>402 

101 

OOC* 

0003 

3321 

5*03 

203 

1*2* 

3921 

OOO* 

>404 

302 

30*1 

OOO* 

30*1 

3*o3 

603 

*761 

9001 

030* 

3*0* 

403 

303  1 

3021 

0  303 

5*07 

60  3 

4*61 

3021 

0  006 

>*1C 

*32 

0306 

0  00  c 

2143 

3*11 

303 

OOC  1 

8  0  9  7 

OOCl 

3*12 

1 1  2 

3  C  0  ’ 

2161 

0  30  1 

3413 

303 

OOCO 

3009 

0003 

341* 

302 

3061 

9001 

3  301 

9*15 

30  <c 

31u  1 

9092 

3  3  21 

3*1* 

302 

3121 

9  0  0.9 

30*1 

*417 

902 

31*1 

OOO* 

3  561 

542C 

902 

3164 

0099 

3*01 

»  4*1 

302 

3  2  o  1 

OOO* 

3*21 

»  *22 

U2 

001* 

2136 

0  30  1 

9*23 

903 

3009 

3000 

0900  , 

9*24  III  0000  0000  0000 

9429  092  0090  0090  |000 
9*2*  *03  5301  4?*i  09»1 
942’  90S  lUj  5001  0>02 

9430  *03  3*21  5q21  00o3 

9431  001  0001  0002  00ft4 

9432  001  0004  0003  000* 

9433  003  7742  000*  000* 
943*  *01  *241  *2*1  0003 
9435  004  0009  1427  0003 
943*  203  0C  0  9  *7*1  O009 
9437  *02  3221  0004  0002 
94*0  102  0002  0009  34*1 
9441  H2  001’  2l’4  0001 
9**2  013  220’  3*22  22oT 

9*43  013  220*  3q23  22|* 
94*4  013  2210  Sol*  2210 
9443  000  0000  0000  0000 
944*  013  2 1  ’ 6  3  o  2  3  217* 
9*4’  013  2 1 ’4  3022  217* 
3430  033  2177  3023  2177 
9*31  033  2200  *022  2200 
9432  014  2223  2l73  22l3 
3*33  013  21  76  3-0  2  3  21*6 
3434  *13  2177  3023  2177 
3433  013  2177  3022  2177 
3*3*  033  2200  3023  2200 
949’  01*  2230  2 1  7 3  *2*3 
94*0  09*  0  C  3  0  2236  01)00 

9461  OOO  3000  0000  0*00 

9462  000  OCOO  0000  0"00 
3*63  000  0000  0000  00*0 
>46*  OOO  QOOO  0009  0000 
3*63  000  GOOD  0000  0000 
3*66  000  OCOO  0000  0000 
940’  000  OCOO  0000  0000 
3*70  000  3609  0000  0009 
3471  000  0000  0000  0000 
3*72  000  OCOO  0000  0000 
9473  q>2  OCOO  OOOo  0000 
9* ’*  *13  *0*1  0030  0000 
3* ’>  0  3*  0000  22>3  OOftO 
3*’*  7o*  *101  *0*1  *101 
94”  70*  41*1  *0*1  *1*1 
>300  70*  *201  *0*1  4201 
3301  304  *2*1  1*27  *2*1 
9302  OOO  3000  0000  OOflO 
3303  1 1 2  003*  22**  B00l 
590*  032  OCOO  226*  oOOO 
>903  *32  CCOO  3000  22*0 
9306  302  CCOO  *0*1  *061 
330’  112  031*  225*  o  ,1  n  l 
331C  OOO  0000  0000  QOOO 
9311  013  225*  22*5  2256 

3512  it*  3315  2255  0301 
5313  05*  226*  2266  2256 
351*  302  0000  *0*1  *061 

3513  OOO  3300  00*3  03*9 
3516  a  > a  3-C0  2121  j0<!3 
5917  *33  *33t  am  fc#ol 
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9520  112  1920  226?  0020 

9921  0*4  0190  613?  1*31 

9922  013  2>91  1*31  239l 

9923  013  3022  613?  1432 
9324  092  0000  0000  0000 

9329  9*4  3921*0000  9321 
9326  703  4741  9o*l  3941 
992?  112  001?  2273  0001 

9330  000  0000  0000  1430 
3331  032  0000  0000  0000 
9332  100  0000  0000  2900 
3933  112  0023  2302  0001 
9934  o  92  0000  QOOO  0000 

9939  492  0000  0000  2313 
9336  401  4741  2900  2900 
9937  401  9Q41  2901  2901 

9940  401  9941  2302  2902 

9941  4qi  3921  2903  2503 

9942  112  0002  2306  0001 

9943  000  0000  3000  0000 

9944  492  0000  000G  2317 

9349  913  2306  2322  2306 
9946  1 1 2  0 00 3  23i9  0001 
3947  000  0000  0000  0000 
9990  H2  0004  2303  0001 
9391  032  0000  2327  0000  * 

9332  000  0003  00Q4  0004 
9993  401  4741  2900  2900 
393*  401  3041  2501  2501 
9393  401  9341  2902  2902 
9356  40i  9921  2903  2903 
9997  900  2323  QOOO  2306 
9960  112  0 0 0 3  232?  0001 
3361  304  2474  7763  2474 

9962  112  002?  2331  0001 

9963  032  0000  2374  0000 
9364  094  0130  613?  1431 
3969  013  23*1  1*31  23*1 
9966  013  3022  302*  000* 
996?  492  0000  0000  2347 
9370  ?04  4041  32*2  4041 

9706  09o  0019  0001  6030 
3?0?  0?0  9331  0000  0000 


3971  112  777?  2340  0001 

9372  000  0000  0000  0000 

9373  o 1 3  23*o  000*  23*0 

9974  112  0033  2337  0001 

9373  0 32  0000  2114  OOOO 

9976  432  0000  0000  2364 

9977  4)2  oooo  0000  2332 

9600  300  *001  0000  4001 

9601  112  777?  2390  OOQl 

9602  000  0000  0000  0000 

9603  013  2361  613?  2361 

560*  Q13  2390  1*32  2350 

9609  013  1430  1*31  1*30 
5606  QOO  0000  0000  0000 
960?  al2  0016  2347  OOQl 

‘9610  016  2361  7301  ?6lO 
.3611  092  4001  00*1  4000 

9612  013  2361  1430  2>61 

9613  000  0000  0000  1430 

9614  000  0000  0000  0000 

9619  112  0002  234(  0001 

9616  000  0909  0000  0900 

9617  o 33  2397  2373  2337 

9620  013  2363  772*  2363 

9621  01*  31*9  23*6  2366 

9622  000  0000  0000  0000 
3623  000  0009  0000  0000 
9624  013  2*01  6137  2*01 

9629  092  0090  0000  0000 

9626  709  4 ?* l  9921  2101 

9627  1 12  0017  2376  OOBl 

9630  000  0000  0000  0000 

9631  000  OOOO  0000  0000 

9632  000  6606  OOOO  292* 
9*33  000  6607  OOOO  2929 

19634  016  2*09  7901  ?*10 
9639  092  2500  00*1  2529 
9636  000  0000  OOOO  OOOO 
363?  000  OOOO  OOOO  OOOO 
3640  032  OOOO  23*6  OOOO 
* 

3710  09*  4411  0001  6030 
9711  070  9331  9706  OOOO 


6*03  01*  6*0*  7301  ?*10  «*0*  000  OOOO  OOOO  OOOO 
*60*  032  32*2  00*2  3261  *607  000  OOOO  OOOO  0000 
06*9  096  OOOO  290?,  OOOO 

Note:  introduce  data  into  cells  0243,0244,3242-3261, 
3344,5222-5242,6137 
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BLOCK  VII  if)  ' 

BLOCK  VIII  (p 

*744  094  0000  *771  0000 

2744  |I(  00?*  9771  *074 

9334  000  14*1  0000  14*4 

3339  000  14*0  0000  14*3 

3474  11*.  0143  5394  0001 

\ 

9703  U*  0143  9391  0001 

■  ■ 

PROGRAM  A 

20*4  100  0000  0000  2101 

• 

*0*9  HI  00*3  2024  0001 

3040  100  0000  0000  2101 

3041  li*  00*3  3040  0001 

3*4?  112  0004  4944  0001 

Note:  program  A 
must  be  added  to 

blocks  VII  and 

4344  401  4001  *1*1  2101 
4347  ll*  0004  4344  0001 

VILI 

t)T4  #5f  210t  0001  21** 

4*03  054  3*4?  3X43  494? 

4404  001  ?? 41  7?43  0044 
4403  904  Slot  0044  2101 
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;i** 

603 

1*34 

3132 

*00  * 

12>0 

001 

1**7 

7741 

•  OOl 

ll*» 

JOV  6oo3 

6102 

*  10  2 

1221 

003 

6301 

5133 

i**7 

11?* 

112 

3032 

432* 

3001 

1222 

002 

0008 

1**7 

1**7 

112’ 

09* 

*132 

93*2 

*002 

1223 

0*2 

C  008 

0000 

0000 

ll’C 

33* 

1*33 

6302 

*008 

122* 

*0* 

1*3* 

6001 

l*<l 

ll’l 

335 

31*2 

50*2 

*003 

1223 

702 

3828 

1*41 

1  *6 1 

11-52 

’32 

*309 

6  03.3 

600* 

122* 

112 

0082 

*62* 

890  1 

il” 

202 

*91* 

600* 

1*5* 

1227 

092 

7741 

9032 

301* 

11’* 

’93 

*000 

6030 

*011 

1230 

093 

001* 

7742 

*0Q  1 

ll?3 

*91 

*311 

3193 

3  1  0  2 

1231 

003 

4801 

7762 

*  3  0  1 

11’* 

7a2 

*326 

1*3* 

*015 

1232 

<21 

6001 

7735 

4  0  o  2 

11” 

283 

*01* 

631* 

*020 

1233 

802 

*801 

4002 

*0  o  3 

11*0 

*01 

*02C 

3102 

3102 

123* 

002 

6003 

7737 

*00* 

11*1 

203 

*03* 

*00* 

*02* 

1213 

a  74 

0800 

*4*0 

1*4* 

11*2 

*01 

*92* 

3101 

3101 

1236 

00* 

7736 

6002 

1*«9 

lt*3 

112 

3035 

*3*0 

0001 

1237 

0  96 

*8*7 

*636 

1**7 

11** 

0** 

>101 

0300 

3101 

12*0 

003 

ooi* 

0*16 

*00  3' 

11*3 

0** 

3  l  o  2 

0000 

3102 

12*1 

016 

*6*2 

7301 

7610 

11** 

0** 

3103 

0300 

3103 

12*2 

073 

6003 

0011 

*0  o  * 

U*’ 

30* 

*30* 

3101 

•  02* 

12*3 

00* 

*886 

9135 

6  0  0  7 

ll’O 

30* 

*012 

9132 

<027 

12** 

31* 

46*9 

7301 

76 1  0 

ll’l 

2lJ  3 

1**3 

<02* 

<»32 

1  2  *3 

073 

4007 

0012 

6010 

1 1  32 

2  8  3 

1**3 

*022 

603* 

12** 

013 

6of  2 

0000 

0000 

ll’3 

2  3  3 

*022 

<02* 

60*0 

12*7 

374 

3o*7 

*691 

1*<7 

U’* 

2Q3 

*022 

<022 

60*3 

1  2  5fl 

3*4 

*010 

*696 

1*63 

il»3 

‘91 

*032 

310* 

310* 

12*1 

015 

7763 

6002 

0  Q  00 

11** 

*91 

*013 

3103 

3 1 0  3 

1232 

0  76 

cooo 

*633 

0000 

ll” 

‘01 

*0*0 

313* 

313* 

1233 

901 

9*14 

6010 

1*43 

:i*o 

*01 

*0*3 

3132 

3107 

123* 

0  36 

0090 

*636 

0000 

11*1 

112 

6010 

*3*2 

3001 

1233 

001 

0*17 

6010 

V*<3 

:  t*2 

332 

COOO 

i  130 

0009 

123* 

003 

96*1 

50*7 

1*66 

•  1*3 

31* 

636* 

7301 

7*10 

1237 

0*6' 

0000 

*723 

**63 

-1** 

373 

»10* 

903* 

*001 

1 2  •  0 

096 

0080 

*777 

0000 

;  1*3 

1 1  2 

0991 

*363 

8001 

12*1 

500 

5033 

0000 

1*90 

:  i*6 

332 

*10* 

1330 

0030 

12*2 

112 

0019 

*661 

OOQl 

'.1*2 

37* 

3142 

*321 

1**3 

12*3 

500 

5007 

0000 

1*30 

1 1  ’a 

03* 

*001 

*322 

1*6* 

1 2*  * 

112 

ooH 

*663 

0001 

!  1  2  1 

1  0  2 

0  *  l  2 

3  0  3  l 

l*** 

12*3 

003 

3031 

3  3  3 1 

3130 

l  l  ’2 

1C  £ 

3  1  0  7 

1030 

0000 

12*4 

.30* 

7  7|i  1 

9132 

1*23 

l  l  2  3 

J’* 

3142 

*373 

1**4 

1  2*7 

112 

OOl  * 

*666 

oool 

112* 

1  SO 

*032 

*376 

1**3 

1270 

002 

7761 

5032 

600  1 

1 1  29 

307 

3417 

3  3  J  2 

1*43 

1271 

009 

0*1  6 

6001 

1  *49 

11’* 

393 

*  1  0  2 

33*7 

1*67 

12’2 

013 

3  0  3* 

7722 

0000 

1 1  22 

39* 

3103 

31*2 

1**7 

1273 

976 

ooco 

*717 

1*70 

i?:a 

002 

3300 

1**2 

1**7 

127* 

002 

5031 

0000 

0000 

1201 

33* 

0300 

*223 

0030 

1273 

076 

0  0  0  0 

*709 

1*43 

1292 

30* 

5  0  6  7 

30*6 

3133 

1276 

002 

3000 

3031 

i*<7 

1293 

33* 

3133 

3135 

*001 

1277 

036 

573* 

*725 

1*47 

120* 

332 

2261 

son 

6001 

1300 

09* 

0102 

3033 

6010 

1203 

392 

9909 

3330 

0000 

1331 

03* 

0  0  74 

*010 

*00  1 

1206 

ICO 

3303 

3393 

313* 

1332 

009 

50*7 

90  31 

*002 

12  0  2 

112 

3  0  0  2 

66  3  6 

0001 

1303 

002 

50*7 

*001 

*  0  0  3 

»  2  1  C 

J** 

6901 

3  398 

3137 

130* 

00* 

6002 

*003 

1*67 

12'.  1 

392 

3903 

3132 

31*3 

1 30  3 

30* 

l  *<7 

1**7 

*00* 

1212 

3  9  C 

*  1  33 

3038 

31*3 

130* 

00* 

77*1 

*00* 

1**3 

1213 

3  3  C 

27*1 

3098 

31*1 

1397 

302 

oooo 

1**3 

1**3 

12'* 

31* 

*6  l* 

6*62 

44*3 

1310 

032 

0080 

0009 

0000 

1213 

OOC 

3132 

3000 

1*4* 

1311 

*0* 

1456 

1*4  7 

*00l 

1  2  1  • 

COC 

3163 

3800 

1  **• 

1312 

302 

3  0  2  3 

*00  1 

1  *6  1 

12” 

ca  s 

3q  46 

5137 

1**7 

1313 

112 

C  002 

*711 

oool 
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1314 

001 

1467 

0034 

1466 

1340 

013 

47J6 

3444 

4736 

1313 

002 

1466 

3031 

1466 

1341 

000 

0000 

0000 

0000 

1316 

0  96 

0000 

4723 

1447 

1342 

112 

7777 

4333 

0001 

1317 

002 

7761 

3033 

6001 

1343 

0  32 

6764 

0  00  0 

7613 

1320 

009 

6001 

0417 

1467 

1344 

000 

7304 

0  0  0  0 

7341 

1321 

0  56 

3734 

4723 

1447 

1343 

0  96 

0000 

00  30 

0000 

1322 

000 

0(H)  0 

0000 

0002 

1346 

034 

0130 

4736 

6000 

1323 

000 

0000 

0000 

0001 

1347 

013 

4730 

6000 

4750 

1324 

000 

0014 

0000 

0000 

1330 

030 

0000 

4731 

0077 

1323 

013 

5733 

7722 

5733 

1331 

033 

6616 

4741 

6001 

1326 

013 

5732 

5733 

0000 

1332 

034 

0114 

6001 

6002 

1327 

076 

0000 

4734 

1471 

1333 

0  34 

0030 

6000 

6003 

133o 

001 

0330 

0000 

0000 

1334 

013 

4762 

6000 

4762 

1331 

076 

7761 

4734 

1471 

1335 

013 

4762 

6003 

4?02 

1332 

000 

3734 

0000 

1471 

1336 

013 

4760 

6002 

4760 

1333 

000 

0000 

0000 

0000 

1337 

013 

4762 

3444 

4702 

1334 

013 

4334 

4724 

4334 

i960 

112 

7776 

4737 

0001 

1333 

0  3  2 

0000 

0000 

0  00  0 

1361 

016 

4762 

7301 

7<10 

1336 

300 

1440 

0000 

0420 

1362 

032 

0000 

0042 

0000 

1337 

112 

0033 

4736 

0001 

1363 

036 

0000 

0  0  3  0 

0000 

1363 

032 

0000 

0000 

0000 

1375 

000 

0000 

0000 

5733 

1366 

100 

0000 

0000 

9100 

1376 

032 

0000 

0745 

0000 

1367 

112 

0010 

4766 

0001 

1377 

002 

3031 

0000 

oooo 

1370 

g9« 

0000 

4523 

0000 

1400 

076 

0500 

3002 

0000 

1371 

002 

0000 

7761 

3734 

i4n 

032 

0000 

46«1 

0000 

1372 

02l 

0330 

7753 

9732 

1402 

002 

5033 

0000 

0000 

1373 

06l 

7731 

3732 

37*32 

1403 

076 

0000 

3001 

0000 

1374 

033 

3732 

7732 

5732 

1404 

032 

0000 

4304 

0000 

130  2 

000 

0733 

0000 

0  0  90 

150  3 

036 

0736 

0»40 

0051 

Note:  introduce  data  into  cells  1364,  66 16,  and  1345 
(except  for  block  IX) 
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APPENDIX  2 

CALCULATION  EXAMPLES 


Example  1 


Assignment  of  numbers  and  instructions  for  calculating  the  aerodynamic 
coefficients  of  a  complex  body  {see  Pig.  7) 

1  .  Numbers 


it  I 

Addressl  Number  IjAddress  Number  I  Address  Number  Parameter 
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2.  Instructions 


Address 

Instructions 

Address 

Instructions 

5201 

000  0000  0000  0000 

6615 

000  0000  0100  0000 

5202 

000  0000  0001  0000 

6616 

000  0000  0007  0000 

5203 

000  0000  0001  0000 

6617 

000  0000  0001  0000 

5204 

000  0000  0001  0000 

6620 

000  0000  0001  0000 

5205 

000  0001  0000  0000 

6621 

000  000  0001  0000 

5206 

000  0000  00(0  0001 

6622 

000  0000  0001  0000 

5207 

000  0000  0001  0000 

1364 

000  0000  0006  0000 

3.  Block  II;  numbers  and  instructions  to  block  II 


Address 

Numbers  and  instructions 

• 

Address 

Numbers  and  instructions 

0313 

-1 

5225 

000  0134  0152  0170 

1420 

10» 

5226 

000  0134  0152  0170 

5222 

000  0134  0152  0170 

5227 

000  0134  0152  0170 

5223 

000  0134  0152  0170 

5230 

000  0134  0152  0170 

5224 

000  0134  0152  0170 

4.  Blocks  IV,  IX,  instruction  to  block  IX 


Address 


Instruction 


1345 


000  0000  0000  0000 


5.  Punched  card  with  control  sum  equal  to  zero 

6.  Numbers  (without  address  code) 


n- 


c 
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7.  Punched  card  with  control  sum  equal  to  zero 

Results  of  calculation  of  the  aerodynamic  coefficients  of  surface  6 


with  Np  *  2000 


Parameter 

•  Calculated  values 

of  the  parameters  - 

without  the  use  of  block  II . 

using  block  II 

1 

-j - 00 

25333314 

c+«i  j 

H - 02 

38246653 

[ 

+  +—  03 

43034721 

-H —  02 

16433219 

'"•+«  t  | 

A - 02 

29514657 

+  +-  01 

14184753 

-L  09  S94H7S9I* 

J.  m  Qfil  1R55I 

C-6i 

-f +—  03  13349649 

+ -  02  63304634 

1 

_L  09  IfiQQfilfifi 

_L.  09  38108933 

1 

+  -i —  03  59884535 

+H —  02  60868448 

m-8i 

+■ -  02  47080428 

H -  01  12316268 

-H —  06  10030515 

+4 —  01  16308264 

+ -  03  23085248 

-H —  00  00000000 

'roi  | 

+• - 05  40420692 

-H —  00  oooooooo 

-H —  03  18622387 

+  H —  00  oooooooo 

+H —  04  71240455 

++—  oo  oooooooo 

mr*l  j 

-H —  03  42093157 

+  J —  oo  oooooooo 

+  +—  04  97449755 

+-|-  oo  oooooooo 

[ 

-i -  00  25881075 

4 -  00  34845169 

c*i  | 

1 

-1-  09  .VW5910Q 

-L  At  10155198 

O-  09  10830457 

-1-  09  .I'VUUTfin 

( 

-H —  02  23134077 

4H —  02  77301667 

n«l 

1 

q -  02  72386769 

4 -  01  15267734 

1 

-H —  01  14282303 

4-4—  01  30493017 

Example  2 


Assignment  of  numbers  and  instructions  for  calculating  local  flows  to 
the  inner  surface  of  a  cylinder  (see  Pig.  8) 

1 .  Numbers 


Address 

Number 

Address 

Number 

Parameter 

Address 

Number 

Parameter 

0420 

0 

0440 

1.0 

— 

0336 

2.0 

du 

0421 

0 

0441 

0 

— 

0337 

-0,01 

Ax 

0422 

0 

0442 

0 

0340 

3,0 

•*21 

0423 

0 

0443 

0 

— 

0341 

-1,0 

J> 

•*12 

0424 

1.0 

0444 

1.0 

— 

0342 

1,0 

*72 

0425 

0 

0445 

-3,0 

— 

0343 

—  1,0 

•*?*. 
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Continued 


Values 

of  index 

k 

Calculated  | 

values  of  i] 

Values 

of  index 

k 

Calculated 
values  of  nrt 

1 

• 

-++  01  1624I04S0 

11 

— H-  01  110155920 

2 

+++  01  166203118 

12 

+++  01  102991993 

3 

— H-  01  152718078 

13 

— h+  01  100126422 

4 

+  ++  01  151453855' 

14 

+  +  +  00  967551620 

5 

— ++  01  140244416 

15 

'  — 1-+  00  855457225 

6 

+++  01  147071217  j 

16 

+  ++  00  808259585 

7 

,  ‘  —++  01  129624947 

17'  j 

r++  00  751790980 

8 

+++  01  126759376 

•  18 

'++■(-  00  644753475 

9 

— (-+  01  1 19342604 

19 

•  — H-  00  589970500 

10 

+++  01  116561314 

20 

+++  00  473662030 
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Example  3 


Assignment  of  numbers  and  instructions  for  calculating  flow  parameters 
In  the  outlet  section  of  a  complex  channel  (see  Pig.  9) 

1 .  Numbers 


Address 

Number 

fi 

•i 

Address 

Number 

Address 

Number 

Parameter 

0420 

i  j 

0453 

0  i 

0506 

1 

_ 

0421 

i  ; 

0454 

0 

0507 

0 

— 

0422 

i 

0455 

0 

0510 

0 

— 

0423 

-0,12  ! 

0456 

2 

0511 

0 

— 

0424 

o 

0457 

0  i 

0512 

0 

— 

0425 

0 

0460 

0 

0513 

0 

0426 

0 

0461 

o  ; 

0310 

100 

H 

0427 

—i  ! 

0462 

0 

0107 

8000 

0430 

o  ! 

0463 

0  { 

0313 

1 

0431 

i  ! 

0464 

1 

0332 

0 

— 

0432 

0 

0465 

—i  ! 

0334 

1 

Po 

0433 

0 

0466 

0 

0335 

3,14159265 

s* 

0434 

0  j 

0467 

i  ! 

0336 

2 

0433 

1  i 

0470 

o 

0337 

0 

•*ii 

0436 

I 

0471 

0 

0340 

5 

xn 

0437 

0 

0472 

0 

0341 

-2.0 

-tl 2 

0440 

0 

0473 

1 

0342 

2,0 

*22 

0441 

0 

0474 

1 

1 

0343 

-2,0 

*?3 

0442 

0 

0475 

0 

0344 

2.0 

•*23 

0+43 

0 

0476 

0 

0243 

2.0 

— 

0444 

0 

0477 

—1 

0244 

2,0 

0+45 

:  6,25313531 

0500 

0 

3242 

3,14159265 

! 

0446 

0 

0501 

0 

0253 

1000 

<Vj  p 

0447 

1  5 

0502 

;  o 

0254 

3000 

A/j  p 

0450 

0 

0503 

i  6,28318531 

0256 

1000 

"ip 

0451 

0 

0504 

0 

|  0333 

1000 

p 

0452 

!  0 

0505 

1 

|  0250 

1000 

T m 

6616 

000 

0000 

0002 

0000 

5202  ; 

000 

0001 

0000 

0000 

6617 

000 

ocoo 

0001 

0000 

5222 

000 

0000 

0000 

0000 

6620 

000 

0000 

0001 

0000 

5223 

001 

0000 

oooo 

0000 

6621 

000 

0000 

0001 

0000 

3344 

000 

0000 

0001 

0000 

6622 

000 

0000 

0001 

0000 

6137  1 

ooo 

0000 

0000 

0001 

5201 

000 

0000 

0001 

0000 

j 

4. 

5. 


Block  VI 

Punched  card  with  control  sum  equal  to 


Results  of  calculation  of  the  parameters 
of  a  gas  in  the  outlet  section  of  a 
complex  channel  with  Np  -  3000 


Parameter 


Calculated  values  of  the 
parameters 


